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Zusammenfassung. Es werden die vier Tensorflachen — lineare, reziproke, quadratische 
und reziprok-quadratische Tensorfliche — samt ihren Eigenschaften kurz besprochen, und er- 
klart, warum im Zweidimensionalen der Unterschied zwischen der linearen und der reziproken 
Tensorellipse sich verwischt. Ausgehend von einer kinematischen Erzeugung der quadratischen 
Tensorellipse wird gezeigt, wie man durch Abanderung des Mechanismus zu den beiden bekannten 
Kreisdarstellungen kommt, dem Spannungskreis von Culmann-Mohr und dem Tragheits- 
kreis von Mohr-Land. 


Summary. The four tensor faces — linear, reciprocal, square, and reciprocal-square — are shortly 
discussed together with their properties, and an explanation is adduced why in the two-dimensional 
domain the difference vanishes between the linear and the reciprocal tensor ellipses. Starting 
from the cinetic production of the square tensor ellipse it is demonstrated how, by altering the 
mechanism, the two well-known circular presentations are arrived at, viz. the tension-circle of 
Culmann-Mohr and the inertia-circle of Mohr-Land. 


Résumé. L’auteur discute briévement les quatre surfaces-tenseur, & savoir les surfaces 
linéaire, réciproque, carée et réciproque-carrée, ainsi que leurs qualités et attributs, et il fournit 
explication pourquoi dans le domain des deux dimensions la différence entre l’ellipse-tenseur 
linéaire et l’ellipse-tenseur réciproque disparait. En partant de la production cinématique de 
Vellipse-tenseur carrée l’auteur démontre comment- en modifiant le mécanisme- on arrive aux 
deux représentations circulaires connues, le cercle de tension selon Culmann-Mohr et le cercle 
@inertie selon Mohr-Land. 

In einem kurzen Aufsatz in dieser Zeitschrift hat Jung! seinerzeit die enge Ver- 
wandtschaft zwischen den Tensorkreisen von Culmann und von Mohr (bzw. Land; 
vgl. dazu den zitierten Aufsatz) nachgewiesen. Das Folgende soll ein Versuch sein, 
diesen Zusammenhang noch unmittelbarer anschaulich zu zeigen. Dabei wird zunachst 
etwas weiter ausgeholt, um den Unterschied zwischen den verschiedenen Tensorflachen 
und -ellipsen zu diskutieren, weil diese Unterscheidung oft nicht klar genug heraus- 
gearbeitet wird. Eine riihmliche Ausnahme in dieser Hinsicht bildet Lagally’, auf 
den sich auch die folgende Arbeit zum Teil stiitzt. 

Es sei 2 ein symmetrischer, zweistufiger, dreidimensionaler Tensor mit den Kigen- 
werten /,, 4,, 2;. Legt man das Koordinatensystem in die Hauptachsenrichtungen, 
dann ist 

Die OS 4 Os 
Rhee Or AD Or (1) 


Ferner sei e ein beliebiger Einheitsvektor; dann ist 
3= Ue (2) 


1 Jung: Der Culmannsche und der Mohrsche Kreis. Osterr. Ingenieur-Arch. 1, 408—410 
(1947). ig 
2 Lagally: Vektorrechnung, 8. 197ff. Leipzig: Akad. Verlagsges. 1928. 
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ebenfalls ein Vektor. (Der Punkt bedeute die skalare Multiplikation.) Denken wir 
uns die zu simtlichen Lagen von e zugehérigen Vektoren 3 als Ortsvektoren aufge- 
tragen, so erfiillen ihre Endpunkte eine Flache. Es ist nun umgekehrt 


C= 2S 
und damit ergibt sich fiir die Flachengleichung 
Veena Us 


oder kurz 


a Dn! gic Ls (3) 


YA 


Yay 


TA, 


& reziprokes Maels- Elipsod 


C reziproke Tensortlache 


@ linesre lensortliche 


Abb. 1. Die vier Tensorflachen. 


Die Flache ist von der zweiten Ordnung, die Koeffizientenmatrix lautet nach (1): 
az 0 0 
12=f 0 AAP oO 
0 Os aa 


Die Koeffizienten sind simtlich positiv, es handelt sich somit in jedem Fall um ein 
Ellipsoid, mit den Hauptachsenlangen |A,|, |A,|, |A,|, das unter Umstanden, wenn ein 
Eigenwert Null wird, zu einer ebenen Ellipse abflacht. Entsprechend der Potenz 
von % sprechen wir hier vom ,,reziprok quadratischen Tensorellipsoid‘‘ oder auch 
vom ,,Mafellipsoid’ (Abb. la). Ein Beispiel dafiir ist das Lamésche Spannungs- 
ellipsoid. Kine unmittelbar anschaulich-geometrische Beziehung zwischen e und 4 
ist an diesem Ellipsoid nicht vorhanden. 
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Der Vollstandigkeit halber sei hier auch noch die Fliche 
joes toy ea 
erwaihnt, das ,quadratische Tensorellipsoid“ oder ,,reziproke Mafellipsoid‘‘ (Abb. 1b); 
dazu gehért z. B. das von Cauchy verwendete Spannungsellipsoid. 


Fiir manche Anwendungen benétigt man nicht den ganzen Vektor 3, sondern nur 
seine Komponente in der Richtung von e, das hei8t man interessiert sich nur fiir 
den Skalar 

ER SM Sp (4) 
Die Bezeichnung i? wurde hier in Anlehnung an den Trigheitsradius gewahlt; es 
ist aber festzuhalten, da im allgemeinen Fall dieser Skalar durchaus nicht immer 
positiv zu sein braucht, das hei&t 7 nicht immer reell ist. Trigt man nun den Vektor 


i 
Carr ne (5) 
als Ortsvektor auf, so ist e =7zy und daher nach (4) 
psy = ly (6) 


Wieder haben wir die Gleichung einer Flache 2. Ordnung vor uns, aber diesmal enthalt 
sie den Tensor in der plus ersten Potenz, wir sprechen daher von der ,,linearen Tensor- 
flache“ (Abb. 1d). Dazu gehért z. B. das bekannte Triagheitsellipsoid nach Cauchy 
und Poinsot. Die Hauptachsen haben die Langen 1/J/A, = 1/%,  1/|/A, = 1/%o, 
1///A; = 1/t; und nur, wenn alle drei Eigenwerte positiv sind, liegt ein Ellipsoid vor. 
Beim Tragheitstensor ist das bestimmt der Fall, wiirde man aber z. B. einen Spannungs- 
tensor in dieser Art auftragen, kénnte sich auch ein ein- oder zweischaliges Hyper- 
boloid ergeben; falls ein Eigenwert Null ist, ein Zylinder, und wenn zwei Kigenwerte 
verschwinden, zwei parallele Ebenen. 

Hier ist es ziemlich leicht, an Hand der Flache den Zusammenhang zwischen e 
und 3 zu rekonstruieren. Legen wir im Punkt xy die Tangentialebene an die Flache, 
und bezeichnen wir den Ortsvektor in der Ebene mit ), so lautet nach (6) die Gleichung 
der Ebene 

(ex) =1, 
oder mit (5) und (2): 
Veg =o: (7) 

Man ersieht daraus, da8 der Vektor 3 auf der Tangentialebene senkrecht steht. Wenn 
wir nun unter ) weiterhin speziell das Lot auf die Tangentialebene verstehen, so ist y 
parallel zu 3 und 

ec Mate 8) 

Isl = So = Ta 


Damit sind Richtung und Lange von 3 festgelegt. 
Eine weitere Darstellung ergibt sich, wenn man zu jedem Einheitsvektore eine 
dazu senkrechte Ebene legt im Abstand i vom Ursprung. Ihre Gleichung lautet: 
nies. (9) 
Alle diese Ebenen umhiillen eine Flache; den Berithrungspunkt erhalten wir in be- 
kannter Weise, indem wir (9) quadrieren: 
(Die)? a7 =e Lee, (9a) 
und diese Gleichung nach dem , scharparameter‘‘, das heift hier: nach e formal 
differenzieren : 


2(y-e)y = 2% -e. 
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Das ergibt mit (9) fiir den Ortsvektor y des Berithrungspunktes 
y= >We (10a) 


bzw. 

4=A-e=—1y. (10b) 
Nun kann man aber (9a) auch schreiben: 

Pee: WU t- Wae 

und mit (10b): ° 

Pei A a, 
das heiBt 

yAt-y=1. (11) 
Da hier der reziproke Tensor von % auftritt, sprechen wir von der ,,reziproken Tensor- 
fliche“’ (Abb. 1c). Ein Beispiel dafiir ist das Drallellipsoid nach Mac Cullagh, 
das reziproke Ellipsoid des Tragheitstensors. Die Hauptachsen der Flache haben 
die Langen //A, = i,, V/A, = ta, Y/43 = 13. Auch hier kann an Stelle eines Ellipsoides 
wieder eine andere Mittelpunktfliche 2. Ordnung auftreten. 


Fiir den Zusammenhang zwischen den Vektoren e und 3 gilt hier, da der 
Vektor 3 =7t zum Beriihrungspunkt der Tangentialebene zeigt, wahrend ander- 
seits e auf der Tangentialebene senkrecht steht: genau umgekehrt gegeniiber der 
linearen Tensorfliche, wie es der Reziprozitét der beiden dargestellten Tensoren 
entspricht. 

Die scharfe begriffliche Trennung dieser verschiedenen Tensorflachen wird dadurch 
etwas erschwert, daB im Zweidimensionalen der Unterschied zwischen der linearen 
und der reziproken Tensorkurve sich ziemlich verwischt. Die lineare oder ,,erste 
Tensorellipse“ hat die Achsen 1/)/A, und 1/\/2,, die reziproke dagegen |/A, und |/2,. 
Da die beiden Achsenverhaltnisse reziprok sind, sind auch die beiden Ellipsen 
geometrisch ahnlich (was bei den Ellipsoiden nicht der Fall ist). Allerdings sind die 
Achsen gegeneinander vertauscht, aber auch dieser Unterschied verschwindet, wenn 
man bei der reziproken Ellipse die Abstaénde 7 nicht in der Richtung von e, sondern 
senkrecht dazu auftragt (ein Verfahren, das man im Raum nicht nachahmen kann). 
Dadurch wird auch die Ellipse um 90° gedreht und kommt damit in dieselbe Lage 
wie die lineare Ellipse: sie wird jetzt als ,,zweite Tensorellipse*‘ bezeichnet. Da iiber- 
dies die GréBe der Ellipse eine reine Mafstabfrage ist (entweder |/2 oder 1//// ist 
keine Lange, man benotigt daher mindestens einen willkiirlich gewahlten Mafstab), 
so werden diese beiden Ellipsen, abgesehen von ihrer begrifflichen Entstehung, geradezu 
identisch. So spricht man beim Tragheitstensor von der Triagheitsellipse schlechthin, 
wobei die Achsen tiblicherweise gleich den Tragheitsradien 7, und 7, gewablt werden, 
das heiBt es wird eigentlich der Tensor 2 = 3/M dargestellt, wobei M die Gesamt- 
masse bzw. -flache ist. Dementsprechend ist das Ergebnis 3 immer noch mit M zu 
multiplizieren, damit man die gesuchten Tragheits- und Deviationsmomente erhalt. — 
Wollte man den Spannungstensor in analoger Weise darstellen, so ergibe sich je 
nachdem eine Ellipse, eine Hyperbel oder zwei parallele Gerade. 


Entsprechend ihrer doppelten Entstehung kann man nun die Tensorellipse auch 
in doppelter Art verwenden (Abb. 2; in der Abb. sind an den Achsen die Bezeich- 
nungen t, und 2, zu vertauschen; an den Vektoren der linken Halfte ist der Faktor 
7,7,/t zu ersetzen durch 1/%): 

a) Man zieht einen Radius in Richtung e, legt im DurchstoBpunkt die Tangente 
und fallt darauf das Lot; dann gibt dieses die Richtung von 3 = %-e, und das 
Produkt aus Radius und Lot, gebrochen durch 7,?* 7,2 ist der Reziprokbetrag von 3. 
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b) Man legt parallel zur Richtung e eine Tangente an die Ellipse; dann ergibt der 
Ortsvektor des Beriihrungspunktes die Richtung senkrecht auf 3, und das Produkt aus 
Radius und Lot den Betrag von 3- In diesem zweiten Fall erhalt man 3 gleich in Kom- 
ponenten zerlegt, namlich in Richtung e und senkrecht dazu: die Lange des Tangenten- 
abschnittes mal dem Lot ist gleich der Senkrechtkomponente von 3, z. B. beim 
Tragheitstensor gleich dem Deviationsmoment gebrochen durch die Gesamtmasse. 


Abb. 2. Die Tensor-Ellipse. 


Ebenso wie die lineare und die reziproke, werden in der Ebene auch die quadratische 
und die quadratisch-reziproke Tensorellipse identisch, abgesehen von ihrer Orientierung, 
die wir aber nicht bentitzen werden. Aus dieser sonst nicht beniitzten quadratischen 
Tensorellipse werden wir nun die beiden Tensorkreise herleiten. 

Ks seien A und yu die beiden Eigenwerte des Tensors, x und y die Komponenten 
von e, bezogen auf die Hauptachsen, dann ist 


fA O\ (x ‘he 

EP alo u(y) = (ao) 
das laBt sich aber umformen in 
A ay A— lM i A+ A— yp ~ 
Ta 2 (_3)= ae Payne (12) 
wobei e den an der x-Achse gespiegelten Vektor e 
bedeutet. Damit kommen wir nun zu der bekann- 
ten kinematischen Konstruktion der Ellipse (hier 
handelt es sich immer um eine Ellipse!) durch ein 
Kurbelpaar, bei dem die zweite Kurbel mit der- 
selben Winkelgeschwindigkeit, aber in entgegen- 
gesetzter Richtung dreht wie die erste (Abb. 3). 
Nach (12) hat dabei die erste Kurbel die Richtung 
von e, wir kénnen somit 3 leicht in seine beiden 
Komponenten parallel und senkrecht zue zerlegen. Abb. 3. Kinematische Erzeugung der 
Sie seien mit o und 7 bezeichnet, in Anlehnung an Tensorellipse. 
den Spannungstensor. 

Von dieser Ellipsenkonstruktion ausgehend, erhalten wir die beiden Kreise durch 
eine kleine Abainderung der Zeichenvorrichtung: bisher bewegten sich beide Kurbeln 
tiber das Zeichenblatt, der Punkt A war fest und der Punkt C zeichnete die Ellipse. 
Befestigen wir aber das Blatt am Hebel AB, so zeichnet der Stift in C einen Kreis 
um B, und wir erhalten die bekannte Darstellung von Culmann bzw. Mohr fiir den 
Spannungstensor (Abb. 4a). Befestigen wir hingegen das Blatt am Hebel BC und 
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bringen den Stift in A an, so zeichnet er einen Kreis um B, und wir haben damit 
die ebenfalls bekannte Darstellung von Mohr fiir den Tragheitstensor (Abb. 4b), 
haufig, aber zu Unrecht, auch nach Land bezeichnet. Selbstverstindlich kann man 
jeden der beiden Tensoren auch in der anderen Weise darstellen, was Klotter? bereits 
fiir den Spannungstensor gezeigt hat. 


Abb. 4. Die beiden Tensorkreise. 


Aus dieser Herleitung ergibt sich auch, warum in den Kreisdarstellungen der 
Radius gegeniiber e mit der doppelten Winkelgeschwindigkeit dreht, und zwar im 
Spannungskreis gegenliufig, im Tragheitskreis gleichliufig. Will man auch im 


Abb. 5. Positive Spannungsrichtungen. 


Spannungskreis eine gleichlaufige Darstellung, was fiir die Verfolzung der Zusammen- 
hange oft vorzuziehen ist, muB man den Kreis spiegeln, das hei8t man trigt z. B. o 
und 7+ als Rechtssystem auf, wenn man sie im Schnittbild als Linkssystem definiert 
hat (Abb. 5); es ware zweckmafbig, diese Festlegung allgemein einzufiihren*. In diesem 
Fall entspricht das Paar o, und t,, der allgemeinen Definition (Linkssystem), nicht 
aber das Paar o, und t,,, die ein Rechtssystem bilden. Man soll daher bei der 
Konstruktion des Mohrschen Spannungskreises ein positives t,, tiber o, nach oben 
auftragen, dagegen tiber o, nach unten: dann drehen Schnittflache und Kreisradius 
im gleichen Sinn. Beim Tragheitskreis wird die Gleichsinnigkeit bereits durch die 
Konstruktionsvorschrift (Auftragung des Deviationsmomentes) gewahrleistet. 


® Klotter: Graphische Darstellung zugeordneter Spannungs- und Verzerrungszustiinde. 
Ingenieur-Arch. 4, 354—364 (1933). 


* Vgl. auch Kappus, Polkonstruktion zum Mohrschen Spannungs- und Verzerrungskreis; 
Stahlbau 22, 138—140 (1958). 


(Hingegangen am 22. Dezember 1953.) 
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Bergman’s Linear Integral Operator Method in the Theory 
of Compressible Fluid Flow. 


C. Axially Symmetric Flow and Singularities. 


By M. Z. E. Krzywoblocki, University of Illinois. 


Summary. In this third part of his work on Bergman’s linear integral operator method, 
the author at first discusses axially symmetric flow where the method of infinite approximations 
by a sequence of linear equations is applied. A discussion of the properties of Bergman’s operator, 
singularities, fundamental solutions and generalization of Blasius’ formulae close this part. 


Part Vv. 
Axially Symmetric Flow 


1. General equations, explanation of the method of infinite 
approximation by a sequence of linear equations. 


Perhaps the most basic, and consequently the most important, particular case 
of a spatial flow is an axially symmetric spatial flow, as for example flow past bodies 
of revolution. Let x, 7, y denote the cylindrical coordinates with the z-axis coincident 
with the axis of symmetry of the flow. The velocity vector then no longer depends 
on the rotation angle zy, but rather entirely on # and r, and always lies in a meridian 
plane (plane through the x-axis). As a result, the flow in all meridian planes is the 
same, and needs to be studied only in a given meridian plane. Below we shall apply 
the hodograph transformation. Thus, x, r denote the cylindrical coordinates in the 
physical plane in which the flow actually takes place, whereas the corresponding 
velocity components, wu, v, determine the cylindrical coordinates in the hodograph 
plane. The velocity potential equation of an axially symmetric flow in the physical 
plane has the form: 

Pax (1 ay Px |") =f Prr (1 ay Pra?) wit 2 Par afr =i Q,/T rr 0, (5. 1 .0) 
where ¢ is the velocity potential and the subscripts denote partial differentiation. 
As is known, the following relations hold: wu = 9,, v = y,. When transforming this 
equation into the hodograph plane, we introduce a new function ® connected with @ 
by means of the formula: 

D=xr9,+1O, —/-. (5.1.1) 


Then in the hodograph plane (u,v), the equation for the function ® is of the form 


L(@) —N (©) —0,. or L(®) = (®), (5.1.2) 
where 
L(®)=(1— B)®ut(1 “)®,, +22 %,,, (5.1. 2a) 
a(D,, D, ; 
N (6) = —F (Pip — Put) =— Fe Hee (5.1.2b) 
a=» (g—w— v’), y= (k —). (O12) 


The quantities and g = dmax are constants. We note that in polar coordinates 
we have g = (wv + vy), @ = arctan (v/w) and equations (5.1.2a) and (5.1.2b) 
assume the form 
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L(®)=|¢9,,+4(1 1 \0,+ (1-4) @oolg, (8.1.24) 


N (®) = G+ {q-> sin 6 [®,, (¢ ®, + a6) — (H,9—q*Pe)]}, (5.1.2e) 
with. 
®, = ©, sin 6 + gq Ge cos 6. 


In the solution of the problem of an axially symmetric flow, we shall restrict our 
attention at first to the case of analytic solutions (for the definition of analytic 
functions, see Part VI), and then we shall consider these solutions for complex values 
of the arguments. The method of solution of equation (5.1.2) by means of successive 
approximations will be explained. 

If we should have a plane flow, then, as one may easily verify, the equation in 
the hodograph plane would have the form 


L(G) = 0. (S223) 
Denote a solution of (5.1.3) by ®. Then: 
BO = 4G + 7 —, — —, (Si: by 


where y is the velocity potential. Consider ® as the first approximation to the 
solution of equation (5.1.2). The method to be discussed is based on the fact that 
for those g (2, 7) which yield flow patterns past profiles of interest in applications, 
comparatively simple formulas can be obtained, using, for instance, the operator 
method. Hence, considering ®® as the initial approximation in that part of the 
plane where r (the distance of the flow under consideration from the rotation axis) 
is substantially larger than zero, the higher approximations can be obtained by solving 
the set of equations 

LOM) ce NO ye ee as (5.1.5) 


If the functions &” converge to a limit ®, then in many instances @ will represent 
a solution of equation (5.1.2) with slightly distorted boundary conditions. The 
corresponding 


p= qe, —®, (5.1.6) 


will represent a potential in the hodograph plane of an axially symmetric flow past 
a profile which is obtained by a slight distortion from the boundary profile corresponding 
to gp. By some additional considerations, these distortions can be diminished 
considerably. 

In the most general case, the equation L (©) = N (®) is a non-linear one. Thus, 
applying the method of successive approximations (5.1.5), we replace the non-linear 
equation by a sequence of linear non-homogeneous equations. An important task 
which arises is the development of a procedure suitable for solving these equations, 
for practical purposes, at least for a small number of approximations, say one or 
two. As is known in the study of linear equations, say in the case L (WY), the fact 
that if Y,,» = 1,2,..., are solutions of L = 0, the linear combination 2, c, ¥,, 
where ¢, are arbitrary constants, is also a solution of L, is of great advantage for 
numerical applications. Similarly, this fact permits us to prepare, once and for all, 
tables of a set of particular solutions Y, which can then be used conveniently in 
determining a solution satisfying prescribed initial or boundary conditions. Below, 
we shall show which tables can be used in the initial approximations, in order to 
enable one to determine a complete set of particular solutions. Moreover, such tables 
can be computed for those values of the arguments in which we are particularly 
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interested; in addition, this procedure gives us a certain insight into the variation 
of the solution if some parameters, for instance the boundary data, vary. 


In our particular case, a considerable simplification may be achieved if the ex- 
pression N (®"~) in (5.1.5) is replaced by 


(n—1) (n—1) 
a(® oo”) 
2 a 
NV. (Gi) — : . 
Ny (®”) (n—1) (n—1) a(u, v) , 
@D + @ 3 


vmax » min 


(On Laa) 
where ocala and oe” are the maximum and minimum values oe’ in the 
domain in which the function is considered. Below, we shall discuss the application 
of such a procedure in our case. 

As mentioned above, we shall limit ourselves to the study of analytic solutions. 
The advantage in making this restriction and in considering the equation in the 
complex domain is that we can simultaneously consider the sub- and the supersonic 
regions and obtain solutions for the mixed case. 


2. The determination of @, 


In this section, we shall determine the solution ®® of the equation L (®) = 0, 
corresponding to flows past profiles of interest to us. As it has been shown in previous 


sections, the stream funktion Y of a twodimensional flow pattern can be represented 
in the form 


Gael 
0 = tn | BZ, 0i(5 20 —8)) 0 — eva (5.2.1) 
t= —=1 
Z=i1(g +10, Z=A(q)—i8, (5.2.1a) 


where # is the generating function of the integral operator, and f is an arbitrary 
function of a complex variable. Here we may use alternatively Z = 4 — 10, 


J —=j,+70, where q 

A(q) = \ [ao — (y + al [ae? — y gg? dq. (5.2.2) 
If, instead of considering the “‘exact’’ compressibility equation, we consider the so- 
called “‘simplified’’ compressibility equation, then the generating function H becomes 


a hypergeometric function. This was indicated in Part IV. As is known, in this 
particular case, the quantities 2 and H are connected by the relation: 


A e=0(— A) (5.2.3) 
where c¢ is constant. 
Considering plane flows of compressible flows, the generating function FH is of 


the form: E =c(—2A)"F (1/6, 5/6, 1/2, X), (5.2.4) 
X =P[cH'(q) +40 — Al 2H), (5.2. 4a) 


where F' is the hypergeometric function and 4, is a conveniently chosen constant. 
Let us keep in mind that in various parts of the plane, the hypergeometric function 
has to be represented by different power series developments. Concerning the 
function f, it may be chosen arbitrarily. For example, in order to obtain flow patterns 
past an elliptical shape, the function / should be of the form [#, pp. 217 and 218]: 
1 “Ie SA SN 

p= 2 Te a | + & expZ— R?U | §’ 

where R = const. An analogous pattern will be preserved in the present case. 


Lote Bergman and B. Epstein: Determination of a compressible fluid flow past an oval- 
shaped obstacle. J. Math. Physics 26, 195—222 (1948). 


(5.2.5) 
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Equation L (6) = 0 may be solved by means of several methods. Applying 
the method of separating the variables, one obtains solutions of L (6) = 0 in the 
form of an expression which represents an approximation 


N bd 
DA, wilexp (iv ) F(%, 8.784), (5.2.6) 
y=—N 
” ie es = 
a= ZrP+(4y7+4, B=ZrE+(4y)4-A, Y=rH1, (6.2.68) 
4A=—vgtwgter—rvytt+(4y)4. (5.2.6b) 
The symbol F (x, B,...) denotes the hypergeometric function. It remains for us 


to calculate the coefficients A, which correspond to the stream function (5.2.1) with f 
from (5.2.5), say. At first, let us note that equation (5.1.4) may be transformed 
into the form [analogous to equation (5.1.6)]: 


gp =4 @) — G0), (5.2.7) 


Our next step is to find a relation between yp and y. From the elementary relations 
between velocity potential and stream function, we have: 


P= ero gy, (5.2.7a) 
go = — oo e7(1 — M*) q+ yh. (5.2.7b) 
From the relation (5.2.7), we have: 
go) = ¢ G3 — 6, (5.2.7¢) 
gp = q BY + BO — OY = qo. (5. 2a 7) 
Combining (5.2.7a) through (5.2.7d), we get: 
q D) = — ogo (1 — M*)q> ph, (5.2.8a) 
(0) (0) 
q Bij — BE = ego qy, (5.2.8b) 
where 
0 = 0 (9) = ell — ye 27 Q?]Ve, (5.2.8c) 


and cy = constant, denotes the velocity of sound at rest. Suppose that we have 
found the values of the function y by means of the integral operator method, for 
instance. Then, using (5.2.6), we express © by means of a polynomial or a series 
dD Byqlexp (ir 6) Fy, (5.2.9) 

y=—Co 
with N finite or N —> co. Assume also that B, Set. i.e., that B, and B_, are 
conjugate. The coefficients B, can be found from (5.2.8b) with the use of (5.2.1): 


+1 
| 2 (Z,Z,0)-f(52(—e)) (ie) eat (5.2.10) 


I 


( é 
q co, ~ Dy) = CoQ orga 


where H and f are functions (5.2.4) and (5.2.5), say. The convergence of the 
series (5.2.9) is assured if the solutions represented by (5.2.9) are supposed to be 
analytic functions of q and 6 in the interval — 2 < 6 <2. In particular, it is possible 
to insert for q a specific constant value, say g = q). This may be done since in the 
following, the expression (5.2.6) will always be considered in a comparatively small 
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neighborhood, in which q varies only slightly. If, then, (5.2.9) is approximated, 
for practical purposes, by a polynomial in the sense that 


t ee N 2 
, \ & B, qo” exp (iv 6) F),, — >’ B, qo?lexp (é 6) Wii d0—-tmin., (5.2211) 
Sage | = OO y=-—N 
then we may show that 
Bye) — By= A, (B,), BO) = B= A,@)[B, 014, > + 0; (5.2. 12a) 
+1 
a 7 Z | ¥! 
A, = —He(1 — acing) ph) | B(Z,Z,t)-{(>4(1—2)] are ey at , (8.2.12b) 
—1 ; C10) 


By = 47 0 (G) Wol2 40? Fo” (407/2) +P Fo (G0?/4?) + 40? Fo (907/4?)], (5.2. 120) 


25 uy i) ‘Uke pe 
A,=¢ ong |i] | E(Z,Z,t)- 
—7 al 
1(320—0)} 00) exp(— veya : (5.2.12d) 
Y= % 


By = 22190 (9) G0" [P Fy (Go°/9?) + 2 G0 Fy (G0/P)]- (5.2. 12e) 


In previous sections, we discussed widely the application of certain tables to the 
calculation of y. These same tables can be used to calculate py (5.2.1), and, in turn, 
for the computation of further tables, for the evaluation of the right-hand side 
of (5.2.12). 


3. Determination of the solution of the equation in the transonic case. 


As indicated above, it is convenient to approximate the set of equations (5.1.5) 
by the equations (5.1.7). As it will be explained in detail in the following, this permits, 
for certain values of n, the use of tables which can be computed once and for all. The 
set (5.1.5) can be approximated by (5.1.7) only in a domain where the variation 
of ®,”—» is comparatively small. This item complicates the procedure a little. 
In the following, the determination of ©“ will be discussed. The same procedure 
can then be repeated for larger n. 


The domain B, in which ® is considered, is divided into a number of sufficiently 
small domains N,, so that in every neighborhood N, the variation of ®, is 
sufficiently small. In order to determine @, it will be necessary to use different 
procedures, depending upon whether WN, lies in the sub-, super-, or transonic region. 
In this section, we shall consider an NV, which lies in the transonic region. 


In order to solve the equation 
L (®™) = Np), (Ga) 


in a neighborhood N,, of the sonic line, it is convenient to replace g by a new variable 
Q=78-W +E (5.3.2) 


Notice, that Q@ = 0 corresponds to points on the sonic line, i. e., we consider Q to 
be small in the transonic case. The expression L (®) assumes then the form 


L(®) =(y 4) v2—@ (+9]7 (47 VF — QP @ + Q) Poe — 27 —Q): 
IvV(P@+Q) +O (y +1] Ge +9 (vy + 1) Doo}- (5.3.3) 


242 M. Z. E. Krzywoblocki: 


We also write 
co Quen 


yytYAGCG+WDYEG-MIM=LY DS 6, rQtexp(s46), (5.3.4) 


s=0 A=—(2N +1) 


and now apply the method of power series to find ®. The result is 


oo 2N+1 
BY = Ay (0) +4, (0)Q +H, + |(Nnys)? 7 Ma+s,cexp (¢46)| Q"*5, (5.3.5) 
n=—1 A=—(2N +1) 


where ®, is an arbitrary solution of (5.1.3), and A, (6) and A, (6) are arbitrary 
analytic functions of 0; 


2N+1 
Gp Uy Nae oyotd -AR(O) a Gy x exp AD), (for mi 2, (5.3.6) 
== (20N fl) 
Here 
m 
NG en ES (S325) 
pei 
and 
Mns3a= 
Catia Te Uy Ve 0 ) Oe a: 000 
Cn, Spiga ch eta Chea ael, a eeamnae (us eee 000 
Cu 0 Sohail gl piel ial Viger Od Foe 000 
ee ae ene Ter CO Ee en ES Lhe, RY A ws oC os c ‘ea 
C34 (0) 0 0 oe ee ee ee ee 0) 0 S, ihe Us, Voy 
(Co, ie Via ay ,4) 6) O Ou ef a 0), 4) Macleay et bye 0) 0 (0) Si Te U; 
(Cup oa Conn ola ted, Giada ree eae Oh OO Ol oS aeaas 
(Co, — T_14,,,) ORO HARE SHER Ore, OF Fe Ot Oe, ESS 
for n =2 
Here 


Sy = 4(m +38) (n + 2) yg, 

T, =[4(n + 1l)yee — (8m + 10) G4] (n + 2), (5.3.8a) 
U, = [(4n — 2) @ — (8m + 2) @7] (n + J), 

Vin (Ae = 6) = 4)yR ye EA. 


Similar formulas hold for M,,, M;,, and M,,. In the computation of these 
constants, tables which can be prepared once and for all greatly facilitate the work. 
For the c,, we obtain the formulas 


N 


CsA oy (2 ie pag | (A, A—1 EY, y=] eA eee if |A| = 2N, 
y=—N 
NV. 
Cy.q = (20) aA, (An ay) : Seay 
if A= 6, (2) orl Pa ae ee ee eee (5.3.9) 


Here the ines are constants which can be determined once and for all since they 
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os derived only from the coefficients of F,,. Therefore, to obtain a general solution 
@, of (5.1.5) in an NV, situated in a transonic region we write 


’ GY — @ + GO), (5.3.10) 
where OD) — 7X) gh Ff, exp (iv 0) is given by (5.2.9) with B, replaced 


vy=—co 
by arbitrary constants X“ which depend upon the specific solution and whose 
determination will be discussed below. On the cther hand, in the solution ®® which 
does not involve any arbitrary constants depending upon the specific solution, the 
B, are completely determined (see 5.2.12), since we start from a known solution. 
@®” is determined as follows: Using the given B,, we compute c, , (see 5.3.9), which 
are substituted into (5.3.8) to obtain M,, ,. Using these values, we finally, according 
to (5.3.5) calculate OD, 
The coefficients are given by the formulas: 


BO Say DO Dee vane |, (5.3. 11a) 


ot 
where Lege Des (5.3.11b) 


co 


co 
ED SBE IT BF AL Cie Seaeza) 


j=0 N=max(j,o—1 


if |y| + |u| = 20 +1, o integer, 


co co 2 27)! fore) 
Lo = a ALY Saas i) os (v, 1) (N—p). 
D; wes 2*—F (4 — 4)! F Di Bae ey? 


EAU AG ERY N = max (j,¢ —1) 
-(y+1)-*(F) je, é (5.3.12b) 


if |y| + |u| = 20, o integer, 
n+ 3 


Bye = 3 Oy, 5 ns 3 ( — [ol) (H* — Lal) + 9H (loa + Dol + lel +1) + 
"2 


+2(n +3 —7) pe + lel) (el — 1) +7 (er — 2 |e] — 1) (of — 1) + 
iG —i2n+ 9 -~ hls Dp 7 oe Se Ale aa 0, Bes) 
and 


Bo = oy Ko [( — |r|) ( — lel) Hee (rel] + lel + lel +1), (6.3.14) 
BY”) = (649 & yy + O%y1 Myo) * (0 — |v) (w? — ||) +e (lr m| + | + [el +1) 4 


$2 ot y.0 Oy (MY + [te]) (Ie] — 1) + 2 o%y1 yo (wy — 2 |u| — 1) (lr|— 1), (5.3.15) 

ow o9=1, (5.3.16) 

de nei=inioo tity {ot-(t)oet (2) +4]-[s)e@+(g,) +44]. 
olf) + (Gp) +4 >] -[(a)e+ Gy) 4) a) eta) -4 + 
Ee (5.3.17) 
A= (Ss) PP +IAOF LPH) + Ay (5.3.18) 


All the proofs are given in the Appendix to the version of Bergman’s paper?. 


2 §. Bergman: Determination of axially symmetric flow patterns of a compressible fluid. 
With an Appendix. Harvard University, Division of Engineering Sciences, Technical Report 
No. 16, Cambridge, Massachusetts (1951). 
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4. Determination of the solution in the subsonic and supersonic regions. 


In this and the next sections, we shall describe a method analogous to that in 
section V. 3, but where the subdomain J, lies completely in the subsonic or in the 
supersonic region. A difference which now arises is that we have to operate with 
the new variables, «, p: 

a=6+H(q), B=9—H(), (5.4.1) 
(instead of the variables q and 6). Here, 
gz} ars 2 L Dynan 
(1 ee : (+42)" dete we, [2 le e's fey 2 be 


A (q) = H (q) = — arctan |= 


TG et") (LE OG") 
in the supersonic case; and (5.4. 2) 
L —1 Vly a0 Viren 
A= H (q) =*n (q); na) = (4) (low / pmo ae 
7 — lla ys 5.4.3 
e=| 7 (@—@) | (5.4.3) 


in the subsonic case. The variables « and £ are in the subsonic case conjugate imaginary 
quantities, while in the supersonic region they become two independent real variables. 

Furthermore, in order to solve the non-homogeneous equation (5.1.5), instead 
of a direct attack such as discussed in section V. 3, we shall use the method of so-called 
fundamental solutions. 

The expression (5.1.2a) assumes, in the variables «, 6 the form 

2 
L(®)=|2(1—-4)— 21H @Pe] O59 +{H" @) — (HP a} (©, — ®). (5.4.4) 
Here q and a are known functions of « and f. 

Besides some minor differences, the considerations in both sub- and supersonic 
cases proceed in almost the same manner; therefore we carry out our discussions 
simultaneously. 

The problem of solving the system (5.1.5) reduces to solving, for every n, a non- 
homogeneous equation. As we mentioned before, this will be done by using so-called 
fundamental solutions. We shall first explain this method in the general case of a 
linear equation. 

Let a partial differential equation 


M (®) = ®,, | a®, | b@, ; CO 2). (5.4.5) 


be given, where a, b, and ¢ are functions of « and B which are analytic in « and f in 
a sufficiently small domain, say B. A solution of the non-homogeneous equation 


M (®) = N («, 8), (5.4.6) 


where N («, 8) is a given function which is analytic in B, can be obtained by using 
a fundamental solution 


W («, B; A, B,) = U (a, £3 A;.B) [log (« — A) log (8 —B)} 4+ V (ob ABs ee 


of M (®) = 0 (see). Here U and V can be written in the form of infinite series which 
converge in B: 


sU te Ue Ores (5.4.8) 


where 
B A 
U,= (5) lexp { a(x, Y) a) + exp (( b (X, B) ix) (5.4.9) 
B a 


* R. Courant and D. Hilbert: Methoden der Mathematischen Physik, Vol. I and II. Berlin: 
Julius Springer. 1931 and 1937. 
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and U,,, is given by 


AB 
ER ipeu,, aU, 
— Une. =\\ (am 40-2 400, \da, ap, 
a B 
U, =U, (a, B,, A, B). 5.4.10 
Furthermore, ee 
V chia Ve inl (5.4.11) 


where V, and V,,,, are given by 


(aU—U,)  (U—U,) 


eS. 3 : le 
: (4 — a) (B— Bp) ’ eae 
oe oV. Vi 
em n 0 n 
a ee | | (a 0x, t b op, c Vn) doy dB, 
ap 
Van = Vile By ASB): (5.4.12a) 


A fundamental solution of the equation (5.4.5) satisfies (5.4.5) and possesses the 
property that at the point (A, B), U («,B; A, B) =1. 


The method of a fundamental solution is usually applied to the equations of 
elliptic type (subsonic case), in which case « and f are conjugate. 


In order to be able to apply it to the subsonic and the supersonic case, we continue 
the variables to complex values, and consider our functions in a four-dimensional 
domain Bt. Let 

a +B 
Di > 


and 
eae 


_ ’ 
24 


where y = y, +77, 6 = 6, +76,. At a given point y, y, 6, 6, B(y, +i, 
6, + 7 6,) satisfies simultaneously the equations 
bile + Fr) +A ln +i”, 5, $i) S + Bly +i, b +48) Se + 
+E (yy +t yh, 6, +062) = N (yy +iyy”, 6, +762), (5.4.13) 
as a function of y, and 6,, and 
1\/26 © Sapien wa? a1 73 a cwevan De) tanoalhag 
(3) ($5 — FA lr tiv, of + 463) oy, | B ; 28, 
+O (yy, +92), 0 +452) = N (yy +i y2”, OY? + 4 55), (5.4.14) 
as a function of y, and 0g. 
In the fundamental solution (5.4.7), the quantities x, 6, A and B now have the 
following significance: 


ae = yy + typ” + 2 (0; 4-4 02”), 

B= ye typ = (dy 1 0); 

A=P+id— 6 +i yO, (5.4.15) 
B= +id+ 6 +i, 


where J’ and A are real integration variables. We shall now denote it by W (71) 043 
yO), 6; LP, A) = W (a, Boel tb): 
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According to classical results (see, e. g., [29]), 
® (yy + ivy”, 6, +1 6g”) so 
= JN W Gps coat N(+iy?, A +469) drdA. (5.4.16) 
ID ISE CRIN 
The integration has to be carried out over a circle [? + A? < 1, in the plane y, = y, 
6, = 6. The expression (5.4.16) can be used in the supersonic case. We wish, 
however, to stress that in the subsonic case, the integration domain K lies in the 


plane of the variables of the differential equation y,, 6,, while in the supersonic case, 
K lies in the plane which is perpendicular to y, and 0, of the equation. 


5. Determination of ®” in WN,’ belonging entirely to the subsonic or 
supersonic region. 


In this section, we shall apply the formulas derived in section V. 4 to our specific 
case. Since the explicit expressions for gq and 6 in terms of the variables « and # are 
rather involved, it is convenient to rewrite the formulas of V. 4, in terms of g and 6. 
This will form the main content of the present section. 

In (5.4.2), (5.4.3), the function A (q) has been introduced. In the following, 
we shall need the inverse function 2 (A), 

q= 2 (A), (5.5.1) 

or, even more, 
q=2(A,), A, +tA+iaA. (5.5.1a) 
The function 2 is a function of a complex variable A,. In the present application, 
A, is real in the supersonic region (i. e., 4, = A), and purely imaginary in the subsonic 
region (41, ==27 A)p ie. A= for Sy ane Se, 
Table 5.5.1 for M <1. We note, however, that in our treatment 
we need the values of 2 for complex values. 2 (A,) 


ilar in can be represented in the form of power series. In 
1 1 1 the subsonic case, the power series of 
3 .23333 —.23333 K “ A 2 
5 11619 + 047143 D = pea (la gt [2 Lyie? — o [2, (5.5.2) 
7 07477 | —.010829 and 7'-1 in terms of exp (2A) have been given in 
9 -05435 + .002721 4 49. 75. 76 
11 04234 | —.000669 [> PP. 49, 75, 76]. : 
13 03448 | +.000158 We note that tables for 2 (A,) are available for the 
15 .02896 | —.000038 subsonic region, i. e., for A, purely imaginary. 
7 02488 | + .000009 In the neighborhood of the point 4 = 0, for y = 0.2 
19 -02176 —.000002 (k= 1.4 f l tati jeldsth : 
21 01930 = 1.4), a formal computation yields the representation 
23 01731 A =iCy Tee Te 5.5. 
25 | .01568 ae Be Daas Orbe) 
T = dS ds Gee Pes (5.5.4) 


where § = (— 3.6 A)'* and the c, and d,, are given in the table 5.5.1. We assume 
in the following that tables for the 2 (A,) for complex values of the argument A, 
are available. 


Let us substitute into the expressions (5.4.9), (5.4.10) and (5.4.12), the variables q 
and 6. If 


“a 


1\ A,” (9) @— [g Ay’ (@P 
G (q) = ) 1 oe gy @P (5.5.5) 
Le ele 


* 8. Bergman: On supersonic and partially supersonic flows. N. A. C. A., T. N. No. 1096 
(1946). 
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then we obtain 


Uo(q, 9’; 0, T; A, 6) = 
ao =2 (i A—a(0)) 


= (5) le, [2 | G (o)A,' (co) do] + 


P Py 


(c=2 iA+ fee) 


2 
+ exp [2 \ G (a) Ay’ (a) dol" (5.5.6) 


6 =2 (6A —4(0)) 


pe ge aad. sil ee 


Here and in the following, A, = A, (q), 4=A(q’). Let 


. ° (0) e . j 
(A, 2, 65, 6, P, A) = i[A (r+ 4 0) decal Sache) Baht 
t Se! 


Then 


o. ay? (8) 6, ADs) a AD (a2) A®) (9) 


eens 1 Agee AO OP) ey ee ee tele Pee cee 


= 6 AMC) m1 APO) AMC) A CW) 
Ay” (%) wi? oU, 


dq, d9,, (5.5.8) 


2 
L a — [Ay (¢;)]7.¢;? on 
where 
ray + (0) (1) ae ee | 
0, = 9, +7263”, A; Wp an a dace past 0, 
L 
0 yf 5 ‘ — 
By =F +169 + i(454), APH) = 0, +46 +i (A 42) 8, 
6 IB : a 
I= 3 +40 —1 (AF), AP 0) = 8 — 6 — 10M +4 (4 — 61%, 
B=—T +162, A® (8) =P 4104-4 (A — 14) — 9, 
AP (4) =O — T 4 0D 45 (A 11). 
Further, 
Ay” (gq) U— (g Ay’ (g))2 U U 
Peeseapn sate , Pkseiorrxate Atay) + P= 8) Ug 
i pei OP { 
For (70) + (4— AP] sa REARS) 
and V,,, follows the law of (5.5.8). 
We then get 


W (9,93 91, Ptr Ace a) a 
= 6 (gq, 4/3 Os, Ps A, 0) log [( — 6,8 + (4 — A) + F (g, g's 8 Ps 1, 6. 
(5.5.10) 


In order to obtain the solution of the non-homogeneous equation (5.3.1), where 
N, = N, (©), in a neighborhood N,, which lies completely in a sub- or a supersonic 
region, we apply the method of fundamental solutions, and according to (5.4.16), 
we write 

eee ; 
OO | |W (args Os, £919, 0) 
BM babe CO) ees 1 
Gq’) Ay’ (7) [41 (0) = @ (Av ('))8 17 dg’ aD + @Y. (5.5.11) 
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where N, (©), is considered as a function of q’ and J’ rather than q and 6. Here, 
@®) is again an arbitrary solution. 

By the procedure described in sections V. 3 and V. 4, we can, starting from a 
solution © of an axially symmetric flow, determine in every neighborhood NV, a 
solution, say ®, of equation (5.1.5) for n = 1. As we stressed before, in doing this 
we can use tables of functions (5.5.5), (5.5.6), (5.5.7), (5.5.8), etc., and expressions 
(5.3.9), (5.3.15), (5.3.18), which can be computed once and for all. We remind 
the reader that in every NV, the function N (6) is replaced by Ny (O), see (5.1.7), 


and that naturally c, = (oe + 6.) will change from region to region. 


On the other hand, the tables mentioned above can be prepared for c = 1, and the 
values obtained from the tables have to be divided by ¢,,. 


co 
Each of the solutions ®® involves an arbitrary function BO = YS”? XW gle’? PF, 
y=—co 


[see (5.2.9)]. We note that in (5.2.9), the symbol B, appears instead of x, and 
we wish to indicate how to determine the constants X‘“ in order that the solution 
obtained in this manner will represent but one solution in the whole domain 


N= 2) N, in which has originally been defined. 


We assume here that we have only one domain N,,, say N,, which lies in the 
transonic region, and the remaining V,, say N., N3, Ny,..., belonging to the super- 
sonic region. Since we wish to determine a solution of the axially symmetric flow 
which differs as little as possible from the plane solution @, we determine the X‘” 
in such a way that 

| | [GO — OP dqdo = min. (5.5.12) 

Ny; 
This condition permits us to determine the unknown X‘”. Let us denote (5.3.10) 
into which X‘” have been substituted by ®{?. Let us denote the function obtained 
by &° + @ = op. The function g{? and its derivative assume certain values 
on the line A B of intersection N, and N,, and now we determine the constants X@) 
in the solution ®) so that on the segment AB the obtained value of yg = O©) + 
+ © and its derivative differ as little from (6) and its derivatives as possible. 
This procedure can be continued through all NV, so that we get a solution g™ in the 
whole domain. We should like to note that if the domain involves neighborhoods V 
belonging to the subsonic region, we can determine the constants X™ by using the 
method of analytic continuation, which we shall not discuss in detail at this time. 

The second approximation ~@ can be obtained, in principle, in the same manner 


as yg, We note, however, that presumably in most instances it will be sufficient 
to determine y™. 


6. General remarks. 


We described above a method of starting from a given solution of (5.1.3) to 
obtain functions which can be considered as approximations of (5.1.2). Every 6™ 
‘ 5 7 

can be inserted into (5.1.2); we obtain then a non-homogeneous equation 


L (®) — N (®) = M (u,v), (5.6.1) 


and from the magnitude of M (u, v), one can make some conclusions about the accuracy 
of the approximation in every J,. 


On the other hand, it is of interest to give at least some sufficient criteria when 
the &” converge, as m — oo, to the solution of (5.1.2). A discussion of this subject 
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is beyond the scope of the present work, and the reader interested in this item may 
refer to Bergman’s original paper’. 


Prag: ab Ve 
Singularities. 


1. General remarks. 


As mentioned above, singularities play an important role in hodograph trans- 
formations. In a brief résumé of Bergman’s method given in the present work, it 
is impossible to discuss all the details and items on the subject of singularities. Hence, 
below only general and more important features will be briefly outlined. A reader, 


advanced in mathematics, may find more details in the original literature, particularly 
in °, 


2. Properties of Bergman’s operator in the subsonic range. 


It is obvious that investigations in fluid dynamics require, in addition to the 
study of regular solutions, the investigation of singularities, for example, sources, 
sinks, vortices, doublets, and so on. Since, further, the image in the hodograph (or 
an allied) plane of a flow pattern is not necessarily uni-valent (uni-valent functions 
are sometimes called schlicht functions) various aerodynamical problems lead to the 
study of branch points as well as to the study of the corresponding Riemann surfaces. 

As mentioned in previous chapters we have several procedures for generating the 
functions ® and YW in compressible fluid flow. Chaplygin, who introduced the 
hodograph method into the theory of compressible fluid, applied the method of 
separation of variables. If, in the power series development for the stream function 
of an incompressible fluid, the powers of the speed, qg, are replaced by certain hyper- 
geometric functions of g, then, as Chaplygin has shown, the series obtained in this 
manner is a solution of the compressibility equation. We designate this procedure 
as the Chaplygin operator, P,, sometimes denoted as P,. In 7 and * Bergman 
introduced a new operator which generates solutions ® and VY of the compressibility 
equations. This is defined in general as the operator Pz or P,. Pg can obviously be 
applied to finite sums and in some cases to infinite ones, producing potentials and 
stream functions of a compressible fluid flow. As mentioned above, in a joint 
investigation, Bers and Gelbart®, independently of Bergman, found the same operator, 
which they applied in®. They term the functions obtained 2-monogenic. 

The operator Pg and Pz, can be applied only to a power series development of 
an analytic function of a complex variable. They transform these series into potential 


5 §. Bergman: Determination of axially symmetric flow patterns of a compressible fluid. 
J. Math. Physics 29, 133—145 (1950). 

6 G. S. S. Ludford: The behavior at infinity of the potential function of the flow of a perfect 
gas in two dimensions. Report 23a, Nav. Res. Contr. Ndori 76/XVI NR 043 046, Harvard 
Univ. (1950). — Two notes on the singularities of a compressible flow in the pseudo-logarithmic 
plane. Report 23b, Nav. Res. Contr. Nb5ori 76/XVI NR 043 046, Harvard Univ. (1950). — 
The behavior at infinity of the potential function of a two-dimensional subsonic compressible 
flow. Report 23, Nav. Res. Contr. N5ori 76/XVI NR 043 046, Harvard Univ., J. Math. Physics 
30, (3), 117—130 (1951). ; 

7 §. Bergman: The hodograph method in the theory of compressible fluids. Supplement 
to “Fluid Dynamics” by R. v. Mises and K. O. Friedrichs, Brown University (1941/42). 

8 §. Bergman: A formula for the stream function of certain flows. Proc. Nat. Acad. Sci., | 
U.S. A. 29, 276—281 (1943). 

® TL. Bers and A. Gelbart: On a class of differential equations in mechanics of continua. 
Quart. Appl. Math. 1, 168—188 (1943). — On a class of functions defined by partial differential 
equations. Trans. Amer. Math. Soc. 56, No. 1, 67—93 (1944). 
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and stream functions of both subsonic and supersonic flows. The potential and stream 
function of a supersonic flow each satisfy an equation of hyperbolic type. Thus P¢ 
and P, transform solutions of an equation of elliptic type (Laplace’s equation) into 
solutions of a hyperbolic one. It is well known that the solutions of both equations 
have a quite different character, and therefore these operators cannot preserve various 
properties of the functions upon which they act. 

Operators Pg and Pg can be applied only to functions which in the hodograph 
plane are defined in a circle with the center at the origin. Often we are interested in 
considering stream functions in domains different from circles. The methods for 
analytic continuation of solutions obtained in this manner are studied comparatively 
little. As is known, a power series is defined only for a certain value of the independent 
variable, and hence is convergent inside the circle with the center at the origin of 
radius equal to that value. 

Moreover, complex potentials of an incompressible fluid flow have branch points, 
poles and logarithmic singularities, around points other than the origin, in the — 
hodograph plane. 

Operators, described so thoroughly in the present work, to be denoted henceforth 
as P,, or P3.,, k = 1,2 (operators of the first and second kind), when specialized 
to the case of functions ® and ¥Y which satisfy the compressibility equations, have 
the advantage that they can be applied directly to functions (and not only to their 
power series developments), and they generate solutions which are defined in an 
arbitrary simply-connected domain (and not necessarily in a circle with the center 
at the origin, which is a domain of convergence for a power series). 

These operators generate solutions which are defined in any simply-connected 
domain which includes the origin. They can easily be extended so as to act on functions 
possessing the singularities of various types. around points other than the origin, 
in the hodograph plane. 

The operators P;,, (or P;,) transform analytic functions of a complex variable 
into stream functions or potential functions of subsonic flows, and differentiable 
functions of one real variable into stream functions of supersonic flows, in both cases 
preserving many properties of the functions to which the operator has been applied, 
thus serving as a useful tool in the investigation of flow patterns of a compressible 
fluid. In contrast to operators of type P;,,, which have been introduced in order 
to “translate” properties of solutions of simpler equations into properties of more 
complicated ones, the operators Po and Pz (denoted sometimes as P, and P,) can 
be used only in a few cases as a tool for investigation of the properties of the functions 
which they generate. 


It should be indicated that by using the Chaplygin operator Pg it is possible to 
obtain various important types of singularities of supersonic flows and various 
singularities, at stagnation points. On the other hand, certain singularities, which 
can be obtained easily by using Bergman’s operator cannot be directly obtained by 
Chaplygin’s method. If we wish to represent these functions in a larger domain using 
Chaplygin’s hypergeometric functions, we need several series developments, each of 
which represents the function in another part of the neighborhood of the singular point. 
It should be stressed that in considering the potentials, ©, and stream functions, W 
there exists a basic difference between the incompressible and compressible fluid case: 
In the first case an analytic function f= © +% YW of a complex variable can be 
interpreted as a flow of an incompressible fluid, the flow being in general defined in 


the whole physical (a, y)-plane or the logarithmic (A, 6)-plane, A= log q, that is, for 
all values of 4 and 6. In the second case, the functions defined on Riemann surfaces 
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arise in a quite natural manner, representing flow in the physical plane which has 
physical Significance. But if we limit ourselves to the subsonic case, and attempt to 
develop in the (A, 6)-plane a theory of compressible flows analogous to the theory of 
analytic functions of a complex variable, and in particular to develop an analogue to the 
theory of algebraic functions and their integrals, then ® and ¥ are defined only for 4 < 0. 
To the value 4 = 0, there corresponds M = 1. For M > 1 the flow becomes supersonic, 
the equations for ® and Y become hyperbolic, and the functions have basically different 
properties. As mentioned and used in previous chapters, if one wishes to apply 
operators to the flows which are partially subsonic and partially supersonic, it is 
useful to introduce as the class of functions to which the operator is to be applied, 
functions of Z = 1 +106, orZ = 4 — 1.0, defined in the (A, 0)-plane or (M, 6)-plane. 
For M <1, Z is complex, and for M > 1 it becomes real. 

As has been shown in”, the operators P3,,, k = 1, 2, preserve the location of 
the singular points as well as the order of infinity. They transform branch points 
of analytic functions of a complex variable into branch points of the same order of 
functions satisfying the linear differential equations, so that the operator P;,, can 
be used successfully for generating and investigating these singularities. On the other 
hand, the operator P, , when applied to analytic functions g (Z) of a complex variable 
which posses a pole, yields solutions of differential equations which are infinite of 
the same order as g (Z) and at the same point, but which are no longer single-valued. 
If we apply the operator P;,,, to functions with logarithmic singularity whose real 
part is single-valued, we obtain functions whose real as well as imaginary parts are 
multi-valued. On the other hand, in connection with the transition to the physical 
plane the equation of single-valuedness of at least one of the quantities ® or Y is 
of great importance, and therefore the question arises of defining procedures which 
generate real solutions of the compressibility equations which are logarithmic or 
infinite of the n-th order and single-valued. Briefly, it is convenient to use operators P,, ;, 
only in the case of branch points, defining the correspondence in the case of poles 
and logarithmic singularities in a different way. 

Each of the operators, mentioned above, can be interpreted as a correspondence 
rule which associates with the complex potential of an incompressible fluid flow, 
defined in a domain B, a stream function of a compressible fluid flow which is defined 
in a domain B’. All three operators act primarily on functions which are regular 
in the domain in which they are considered, and they produce solutions which are 
regular in the domain in which they are defined. 

It is natural to require that, insofar as possible, the stream functions and the 
potential functions of compressible fluid flows behave in the pseudo-logarithmic plane 
at singular points in the same manner as the corresponding functions for incompressible 
fluids in the logarithmic plane; that is, they become infinite of the same order, the 
geometrical structure of the stream lines and potential lines in a sufficiently small 
neighborhood of the singularity is essentially the same, and so on. 

Let us now explain why single-valuedness of at least one of the components (that 
is, of the potential or stream function) is of importance for our purposes. Poles and 
logarithmic singularities in the logarithmic plane are, as a rule, images of doublets, 
vortices, sinks or sources at infinity in the physical plane. It is natural to require 
for compressible fluids that at infinity the speed is constant and that the motion 
has the same character as in the case of an incompressible fluid. If the potential or 
the stream function is known in the pseudo-logarithmic plane, we can come back 


AY Sy Boyan Certain classes of analytic functions of two real variables and their properties. 
Trans. Amer. Math. Soc. 57, 299—331 (1945). 
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to the physical plane by using known expressions for x and y, cited in previous sections. 
It follows from those equations, that if the derivatives of at least one of the functions, 
® or W, are single-valued, then, if we move around the singular point, the functions 
a—a(u,v), y= y(u,v) can each increase at most by a constant, say X, Y. In 
most cases we have at infinity a linear combination with constant coefficients of 
these singularities, and the constants have to be so chosen that the quantities X 
and Y vanish, and therefore x (u,v) and y (u,v) are single-valued. This is exactly 
the case for an incompressible fluid. If this property were not preserved, the flow 
would completely change its character. Indeed, suppose that x (u,v) and y (u, v) 
are multi-valued: then the functions wu (x,-y), and v (x, y) are periodic, that is, in the 
physical plane every value w, v, is assumed infinitely often. The question whether 
other types of singularities—in particular those for which both components are 
multi-valued—are of interest in hydrodynamical applications, remains open. The 
first step in answering this question would consist in studying the behavior of flow 
patterns in the logarithmic and physical planes in the neighborhood of these 
singularities. 


3. The behavior of certain types of flows in the case of an 
incompressible fluid. 


Following the line of our approach, that is, applying the principle of correspondence 
between incompressible and compressible fluid planes described in previous sections, 
we first have to describe the behavior in the logarithmic plane of various types of 
flows of an incompressible fluid, and in particular those possessing singularities such 
as sinks, sources, vortices, and so on. 


Let us consider the flow of an incompressible fluid around a closed curve, a profile, 
say, and let us assume that neither the velocity nor the circulation vanishes at infinity. 
The complex potential in the physical (z-plane) plane may be expressed for |z| 
sufficiently large, in the form 


> es : 
w (2) =q2+mlogz+ Dd) a,2-", (6.3.1) 
T 
2 
GW + 9, m + 0, 


= 
where qq is the velocity vector at infinity and the real constant m/i (m is an imaginary 


quantity) is the strength of the vortex at infinity. For the velocity vector q we 
obtain by differentiation with respect to z 


q (2) =4% + m/z = Snag (6.3.2) 
The expression for the complex potential in the logarithmic plane is 4 
W, (2) Sm[2 BF tee (6.3.3) 
+ Power series in 7 (Z Z»); 


where Z=1+76,4= log g, Zy = Ao +4 Oo; ho = log do. The independent complex 


variable Z may also be assumed to be equal to Z = 4 — i 6; theoretically this change 
has no effect upon the considerations presented below. 


In the most general case we can derive the stream function in a form: 
P+ = wD) YOY 4 y2) WA?) 4 YB) EV 4 44) yr, 2) 4 we, 
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If there is only circulation in the physical plane and no source, then yw, and ji 
are determined solely by the circulation. Similarly jo and uo are determined 
solely by the source in the physical plane. This means that + splits into two pairs 
of Singularities (YW), W,2)) and (Wa.2, wa, D). The first pair is absent if there 
1s no source in the physical plane, and the second if there is no circulation. In case 


there are both source and circulation, all the terms in (6.3.4) must be preserved. 
Thus we have: 


(1) (Yo, 89) determines location of singularities in pseudo-logarithmic plane, 

(2) source strength alone in physical plane determines fo and uo, 

(3) circulation alone in physical plane determines 2 and [p). 

Assuming the most general case, we can derive the stream function in a sufficiently 
small neighborhood of Ze iZ, in the form: 


4 
Pe Oe (6.3.4) 
iL 
Cy (Ae MS) (6 oy, (6.3.4a) 


etc., where the dots indicate terms which go to zero as Z — Z, and therefore can be 
neglected in considering the behavoir of Y* in a sufficiently small neighborhood 


of iZ =iZ,. This point corresponds to the point at infinity in the z-plane. 


Thus we see that W, (Z) and ¥*+ have singular points for = Aas One may 


verify that z=2(0,2)=2-+iy, is single-valued at the singular point 05, Ao. 
Following our principle of introducing (in the pseudo-logarithmic plane) for the 
compressible fluid, functions which have a behavior similar to that of the corresponding 
functions in the logarithmic plane for the incompressible fluid case, we see that we 
should consider flow patterns whose stream functions Y possess, as singularity at 
Z=Z,, a linear combination of a logarithmic singularity and two independent 
singularities of the first order. Suppose a flow, which symbolically we shall denote 
as y (®, Y) (of a compressible or incompressible fluid) has a singularity at the point 
(A, 95). If at least one of the functions, ® or Y, is single-valued at this point, the 
singularity of y at (Ao, 09) will be said to be of type S. For example, the flow 
D = log [(@ — 6)? + (A — A,)?|/2, W = arctan [(@ — 6,) (A — A,)+], represents a 
flow possessing a singularity of type S. As we indicated above, this property plays 
an essential role for our purposes. In one of the following sections we shall describe 
a method of generating required singularities of type S in a compressible fluid flow. 


4. The behavior of a subsonic compressible flow at infinity. 


To explain the behavior of a subsonic compressible flow at infinity, let us discuss 
some examples. At the point a = « +78, «, B, real, of the hodograph plane which 
corresponds to the point z = (% +7 y) — co of the physical plane, the stream function 
has a singularity. This fact leads to the study of the singularities of functions 
satisfying equations in the pseudo-logarithmic plane. Let us put down the equation 
of a compressible fluid flow in forms 


Wo + Woo* +4F,, (24) ¥* = 0, (6.4.1) 


or 
Ue ie 0. (6.4.2) 
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A solution of the equation (6.4.2) may be written in the form 


Zn—1 


Z — — > — 
w* (A, 0) = NS ee 2m)! Gin) (9 gy « )-[... [gZndZn...dZy}. (6.4.3) 
0 10) + 


g2n ny ! 
Since g (Z) is an arbitrary function we may choose it so as to fulfill certain conditions. 
If we know that the point corresponding to z — co should be a branch point, then 


we choose g = (a — Z)~: and we obtain a singularity which possesses the desired 
features. Indeed, let us insert g into (6.4.3) and the function 


W (2, 0) = HIm{(a — Zs —Q® [(a — Z's — alle] + 
+ 3Q® ((a—Zyn— ah — Fah} +...) (6.4.4) 


with Z = 4 — 70, is a stream function which is two-valued at point a and becomes 
infinite as (a — Z)~/2. If, however, the point a is not a branch point, that is, if (6.4.3) 
is applied to the function g = : —Z)-, (a pole), say, then 


WV (A, 0) = HI mi{(a — Z) + (1/2) Q™ [log (a — Zo loga] +...}, (6.4.5) 
is obtained, which is not a single-valued function because log (a — Z) is a many 


valued function. Similarly the function g = log (a — Z) will originate a many-valued 
function. The function (6.4.5) can, however, be made one-valued by replacing the 
many-valued term (0 = arg Z) by its mean value (in the sheet under consideration). 
Clearly, this new function will no longer be an exact solution of the equation in the 
pseudo-logarithmic plane, but in many instances it will not differ very much from 
an exact solution. Naturally, this procedure may be refined. 


5. Fundamental solutions in the subsonic range. 


It is, however, of interest from a theoretical point of view to determine exact 
solutions of (6.4.2) which are single-valued and have a logarithmic singularity at 
point a. Clearly, it is sufficient to find functions for equation (6.4.2) which possess 
a logarithmic singularity. A function 


W* (A, 93 dg, 09) = W(Z,Z; Zo, Zo) = 
= A (Z,Z;Z,,Z,) log |Z — Z| + B(Z,Z: Zo, Zo); (6.5.1) 


which considered as a function of Z, Z, satisfies equation (6.4.2) in the whole plane, 
except at the point Z = Lo, is termed a “fundamental solution’ of equation (6.4.2) 
with the affix at Z=Z). Then, clearly, Y = W (Z,Z;a, a) represents a desired 
stream function with a singularity at Z, =a. The corresponding flow can be said 
to have a pseudo-cortex (we may also call it a flow with a “vortex-like” singularity 
at Ao, 09). The names “vortex-like’, “source-like’’, etc., are used although there is 
no direct relation between a source-like singularity and a source in the physical plane. 
In the case of an incompressible fluid (and in the physical plane) the stream function 
of a vortex flow is given by the expression m log (z — 29), and in the case of a source 
the potential function is given by m log (z — 29), m being a constant. We note that 
in general A and B are not analytic functions of one complex variable Z (or Z,). In 


order to emphasize this, they are often denoted as A (Z, Zz sich; Zo); B(Z, Vis Lis, Z.), 
respectively. 
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Tn order to find the functions A and B a new variable will be introduced, namely 


; ay ny (6.5.2) 
Equation (6.4.2) then becomes 
A fundamental solution of (6.5.3) with the affix at Z, = 0, will be a fundamental 


solution of (6.4.2) with the affix at Z —Z,. The functions A and B are given by 
the equations: 


ZZ, “ TATA hi Tie ii un 
A=1—| | F,,dZ,dZ, + \ { Pal \ Fy ZZ, dZ,dZ,+..., (6.5.4) 
00 00 00 
ZZ, # ZZ, ZZ, as a 
Bia | Gd. da. | Pal ("dt dZ| az, dly +... (6.5.5) 
00 00 00 } 
G = (Az/Z,) + Az,/Z,). (6.5.6) 
Not only the function W (6.5.1) satisfies equation (6.4.2); the expression 
owWw* 1 
pe Oot eat ee Co) Mere re) nt (0 Oo)es at 
+ A, log ([A — A,)? + (6 — 6,)?] + As, (6.5.7) 


where A;, 7 = 1, 2, 3, are entire functions, is also a solution of that equation. 
In the case where the pressure-density relation is of the (simplified) form 
p=A-+ B/C, the functions W* and 0W*/06 are 


(1/2) log [(A — do)? + (6 — 0o)?], (6.5.8) 


and 

(8 — Oo) [(A — dg)? + (9 — O57], 
respectively. It is noted that in this particular case new singularities are obtained 
by differentiating W* with respect to A. 

Let us differentiate ¥ = YW = W (Z,Z;Z,,Z,) with respect to the real and 
imaginary parts of the parameter Z. We obtain single-valued functions ¥") and 
W2) with an infinity of the first order. Flow patterns generated by ¥,°% YO), 
wL2) have a simple physical meaning: Y“) yields a vortex at infinity (in the physical 
plane), YW, WC,2) yield components of a doublet at infinity. 

There are instances where it is more convenient to operate with the potential 
function ®, rather than with the stream function Y. Functions ® and W satisfy the 
system of equations 


= Vy, 1(f)¥o= — Og. (6.5.10) 
Introducing the variable 4, we obtain 
Og =02(1— M)hY,, et(l— My) Yo = — @, (6.5.11) 


Differentiating the first equation (6.5.11) with respect to @ and the second with 
respect to A, and applying a few transformations yields the equation 


Gai a diate ae a (6.5.13) 
N =(— 1/8) (6 + 1) 4 — Mf?) — 3/2. (6.5.14) 


pL 


Equation (6.5.12) possesses a fundamental solution of the form (6.5.1); in this case 
a flow with a ‘“‘source-like” singularity is obtained. In the physical plane it yields 
a source (or sink) at infinity. 
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6. Modified integral operator in the subsonic range. 


As mentioned above, in connection with the transition to the physical plane, 
it is of importance to have singularities for which at least one of the functions ¥ 
(stream function) or ® (potential function) is single-valued at the singular point. 
Sometimes it is necessary to modify the operator described in previous chapters. 
If we have a function g (for example, a complex potential) which possesses poles and 
logarithmic singularities, we decompose it writing 


g= gt 2 >A, (2-7) to ae os 2 (6.6.11) 
k=1 k=m 

where A, are constants and the », are positive integers. The function g* possesses 

as its only singularities algebaic branch points. 


Consider the equation (6.4.2). The integral operator may be chosen in the form 
of equation (6.4.3). Applying the integral operator to the first part of (6.6.1), we 
obtain a solution of (6.4.2), say Y%, = I m [p (g*)]. By introducing the fundamental 
solutions Y“ — W in form of (6.5.1), we obtain solutions of (6.4.2) which have 


the same behavior as log (Z — Z,). Adding to WY, the functions with singularities 
described above we obtain a solution of the compressibility equation, possessing 
the same characteristic features as g. We shall refer to this procedure for determining 
stream functions of a compressible fluid as the “modified integral operator’. In 


the case of non-symmetric flow patterns with the velocity vector do eXp: (i Oy) at 
infinity (z — oo), we obtain as a rule the stream function (in the 4, 6-plane) in the 
form 


P= A, PO4 A VOD + A, WD 4 Ww, (6.6.2) 


where Y, = YW, (Z) is a function which has branch points of positive order as its only 
singularities, and which can be obtained by the operator from a suitably chosen 
function g* (Z) of one complex variable. In the case of symmetric flow patterns the 
representation (6.6.2) simplifies. In the latter case we can directly apply the operator 
to the function 


g(Z) = Dd) A(F—Z,)"2. (6.6.3) 
y=—l1 

The introduction of the modified integral operator enables us to develop a theory 
of compressible fluid flow which in the subsonic case strongly resembles that of an 
incompressible fluid (when we use the hodograph method). In the following chapters 
a generalization of two results of the latter theory will be discussed, namely: (a) the 
condition that a stream function defines a flow pattern around a closed curve; and 
(b) a generalization of the Blasius formulas. 


7. Subsonic compressible fluid flows with singularities of type S. 


Classical methods in the theory of partial differential equations yield almost 
immediately two independent flows with logarithmic singularities of type S. Assuming 
the equations of the stream function and the velocity potential of a compressible 
fluid flow in Z-plane or 1 Z = 6 + i A-plane, one can find the two pairs of fundamental 
solutions of those equation in the pseudo-logarithmic plane: 


WE) 2) G&L) G2, 


The superscript L symbolizes that either the stream-function or the potential 
function has a logarithmic singularity. By means of a simple procedure, one may 
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show that the flows y4) (Oi, wu, *)), k=1,2, are flows with singularities 
of type S and that both components O4"), WU of the pied). a. a=) oo are 
single-valued. Infinitely many new flows may be obtained from the first two linearly 
independent flows by a differentiation, namely 


vin, b) ( a” ps *) on pi, =| = 

oo” a 
The components of those flows are derivatives of n-th order of the first potential 
and stream-functions, and they are infinite of n-th order at the point (95, 49). The 
question may arise as to whether for some index n (and hence for all higher indices) 
the two flows y™ and y™2) are not essentially the same: that is, whether functions 
YY) and ¥™2 (and likewise O”) and @) are not linearly dependent. The 
proof that this is not the case, as well as all the other proofs of statements cited in 
the present section, are given in}, 

A problem of considerable interest is that of determining a flow of a compressible 
fluid around a given profile, or at least around a profile the shape of which approximates 
the given profile. Since, in many instances, by reasoning from the incompressible 
case, the approximate image in the hodograph (A, @)-plane is known (the image in 
the logarithmic plane of an incompressible fluid flow around a profile is supposed to 
be used as a first approximation of the image in the (A, 6)-plane of the flow of a 
compressible fluid around a similar profile, it is possible to consider, instead of the 
above problem, the question of determining a flow for a given hodograph. The 
behavior at the point of the hodograph corresponding to z = oo is prescribed. Clearly, 
instead of the image in the hodograph plane, the image in the (A, 6)-plane may be 
used. By means of formulas given in previous sections, it is possible to determine 
a function VY, (A, 6) satisfying the equation of a compressible fluid flow which possesses 
the required behavior at z = oo. Naturally, WY, (4,6) for the point (A,, 0.) 
corresponding to the point z = oo of the physical plane, must have a singularity 
which satisfies the conditions indicated in section VI. 8 (below), in order that the 
flow in the physical plane will be a flow around a closed curve. Function Y, (A, 0) 
is not as yet the required stream-function, since it does not assume constant values 
on the boundary of the domain. But it is possible to determine this function in an 
approximate way by means of a procedure given in ?’, p. 45. 


a Ae (6.7.1) 


8. The transition from the pseudo-logarithmic to the physical plane 
in the subsonic range. 


We shall discuss below the conditions that W (A, 6) defines in the physical plane 
a flow around a closed curve. As has been indicated above, the introduction of 
singularities of type S enables us to define in the pseudo-logarithmic plane, flow 
patterns which have in this plane a behavior similar (in the sense indicated above) 
to that of the flow pattern of an incompressible fluid in the logarithmic plane. The 
problem arises of determining the conditions that the flow defined in this manner 
in the pseudo-logarithmic plane will represent a flow around a closed curve in the 
physical plane. In the following discussion we shall present necessary and sufficient 
conditions to assure this. 

If a flow has a singularity of type S at (A, 09) then either ® or ¥ is a single-valued 
function, and, according to the statements and equations cited in previous sections, 


11 §, Bergman: Two-dimensional subsonic flows of a compressible fluid and their singularities. 
Trans. Amer. Math. Soc. 61, 452—498 (1947). 

122 §, Bergman: Methods for determination and computation of flow patterns of a compressible 
fluid. N. A.C. A., T.N. No. 1018 (1946). 
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a and y increase by constants if we go around the singular point (Ao, 09). The 
constants X, Y, by which the functions «x (A, 0) and y (A, @) increase (if we move 
once counterclockwise around a singularity of type S of a flow y), will be denoted 
by 7-periods. (The letter 7 is suggested by the expression “transition to the physical 
plane.”) These constants X, Y, can be determined for the flows y%, y%”, (for 
details see). Following the principle of correspondence, explained above, we 
introduce stream-functions which in the pseudo-logarithmic plane have a behavior 
similar to that of stream-functions of a incompressible fluid in the logarithmic plane. 
The images of.a flow of an incompressible fluid in the logarithmic plane and that of 
a compressible fluid in the pseudo-logarithmic plane extend in general to infinity, 
since at a stagnation point g = 0 and therefore A= log qg tends to — co. Below we 
shall consider only flows satisfying the following conditions: 

(a) The stagnation point (or points), if any, lies on the boundary (but not in the 
interior) of the flow; 

(b) the flow (that is, the stream and potential functions) possesses only a finite 
number of singularities (including branch points), so that for |Z| > #& the stream- 
function is regular at every point of the domain (call it P); 

(c) the stream-function Y as well as WY; and YZ exist and are continuous on the 
boundary except perhaps at a finite number of points. 


Assume that the stream-function Y of y is defined in a (not necessarily schlicht) 
domain P of the pseudo-logarithmic plane and can be represented there by 


YP = pg) MOY + yo 2) WA?) + 98) WD 4 yy (4) Ws 2) 4 (6.8.1) 


Assume that YW, satisfies hypotheses (a), (b), (c), and that the functions Y», 
WL), ete., are stream-functions possessing at the point A), 6) singularities of the 
type described in the section VI. 7. The point 4), 0) corresponds to z — co. Under 
the hypotheses (a), (b), (c) and some additional assumptions, the following theorem 
can be proved !! and °. 

Theorem. The necessary and sufficient condition in order that the stream- 
function (6.8.1) defined in a simply-connected domain P represents a flow in the 


physical plane around a closed curve is that the following equations are fulfilled 
(called ‘‘the consistency relations’): 


fg) XY + pg’?? XO 2) + Me’ XL YD + by? Xb; 2) — 0, (6.8.2) 
tig? YD 4 pyg®) YU) 4 y (8) YD 4 4) YG 2 — 0, (6.8.3) 
XL), YL" being T'-periods of y),.k = 1, 2, X%®», Y“4® those of y*), 


We may derive another form of the conditions in order that the obstacle in the 
physical plane (which is formed by the image of the boundary curve h in the simply- 
connected domain H in the hodograph plane) be a closed curve. Assume that the 
stream-function W defined in this domain H has a singularity at the point a (the 
image of z= co). It is assumed that the point g = 0 is not an interior point of 
domain H. Clearly, the necessary and sufficient condition in order that the image 
of H in the physical plane be single-valued is that 


| go exp ( 40) 07 {[— (1 — M?) q? Vo +1 Pq] dq + [(%, +2 oq] d0} = 0, (6.8.4) 


ois 1 is any simple closed curve lying entirely in H +h. The symbol (H + h) 
means that at most the curve J can coincide with h. Since the integrand of (6.8.4) 
is a complete differential, the value of the integral does not change if / is continuously 
deformed without leaving H +h and without passing through the singular point a. 
Thus, in particular, if the boundary curve hf is chosen for J, it follows that 
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d¥ = VY, dq + Pod = 0, that is dd = — (¥,/Po) dq along h, since h, being the 
boundary curve of the body, must be a streamline. If this expression is inserted 
into (6.8.4), the condition is obtained that the image of h in the physical plane be 
a closed curve. Notice that H must be simply-connected. 


On the other hand the relation (6.8.4) can also be written in a different form 
which is often more suitable for applications. Let us assume that the stream- 


function Y can be represented in the neighborhood of a (the image of z = oo) in 
the form 


k 
= > A, WO (6.8.5) 
n=0 
where A, are constants, Y, is a function which is regular at the point a and W™ 
are fundamental solutions [see (6.5.1)]: 


W (Gq, 93 qo, 90) = W* (A, 8; Ao, 9), (6.8.6) 


AW) (A, 05 Ags O%) 
W™ (q, 93 do, 99) = —— Fe ae (6.8.7) 
With every singularity W™, n = 0,1, 2,.. 


. there may be associated a complex 
number 


R(W™) =X, +7Y,, (6.8.8) 


which is equal to the left-hand side of (6.8.4) with Yo = Wo” and ¥, = W,”, 
and where / is an arbitrary simple closed curve around a, the direction of integration 


being such that a always lies to the left. The alternative form of conditions (6.8.2) 
and (6.8.3) iss that 


k 
Di An (WO yes 0: (6.8.9) 


n=0 


The evaluation of the quantities R (W™) in the special case, when k = 1, is given 
fis 


9. Generalization of the Blasius formulae in the subsonic range. 


As is known, the behavior of the flow at infinity determines the force and the 
moment acting upon a profile immersed in a fluid. As one may show, the results 
obtained in previous chapters enable one to obtain for the force F and the moment M, 
about a given point 2 = (Xp, Yo) the formulas analogous to those holding in the 
incompressible case. For an incompressible fluid these quantities are given by the 
formulas of Blasius: 


F= Fi 00\(G)Pde, (6.9.1) 
C 

M,o= — Re fae |@—=0)(a} a (6.9.2) 
) 


Here 0, is the density of the fluid, w is the complex potential of the flow and O a 
control contour. Formulae (6.9.1) and (6.9.2) can be expressed in terms of the 
coordinates A, 0 of the logarithmic plane. Bergman and Ludford generalized those 
formulas by taking into account the compressibility of the fluid. 


18 §, Bergman: On two-dimensional flows of compressible fluids. N. A. C. A., T. N. No. 972 
(1945). 
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Moreover, Ludford proved the validity of the generalization of the Kutta-Joukowski 
formulae to the subsonic compressible case: 


ff = Os (6.9.3) 
The moment M about the origin is: 
Mizsiogt? Ds (6.9.4) 


with D, given in ®. 


10. An integral formula representing a subsonic flow inside a domain 
in terms of its values on the boundary. 


The Cauchy integral formula, cited in one of the previous sections, is one of the 
most powerful tools in the study of analytic functions of a complex variable. In 
the subsonic case there exists an analogue of this formula, which enables us to express 
the values of © and ¥W inside a domain B of the (A, 6)-plane in terms of the values 
of ® and ¥ on the boundary of B. The following theorem will be given without the 
proof: 


Theorem: Let y (®, Y) be a flow which is regular in the bounded domain B 
(lying entirely within the subsonic region) whose boundary C is a simple closed curve. 
Then at every point (Ao, 05) of B 


(L, 2) 2Ple @ iL, 2) 
2 @ (Ay, 0) = LH (24))* |[— wre SO 4 yw PEPE) G6. 10.1) 
C 
( (L, alt) = Ye (L, 1) 
2 mY (Ayy 09) = H (2 2g) || — Hy A po ot ds, (6.10.2) 
C 


where 0/0, indicates differentiation in the direction of the inward normal, and 6:2) 
and Y“.) are the fundamental solutions. 


11. A proof of an auxiliary lemma in the subsonic range. 


In the following, the proof will be outlined that in the subsonic range the 
function F (2A) can be approximated by the polynomial 


Fy(24) = Dd) «exp (28), ap = 0, (6.11.1) 
s=0 
in every interval (— oo, 4), 44 <0, and the method of determining the F,, (2x) will be 
indicated. 


If 24 is replaced by log X, then F [2 4 (X)] is a continous function of X in the 
interval (0, X)), X 9 <1, and by classical results, it is obvious that it can be 
approximated by a polynomial in X. It will be seen that it will not suffice merely 
to approximate fF, but in addition to this it will be possible to require that any given 
number of higher derivatives of F be approximated in the interval (— oo, 2») by the 
corresponding derivatives of F,,,. It is, however, of interest to give a more explicit 
form of the approximating polynomials. This will make it possible to determine the 
corrections of Q™) given in previous sections on subsonic flow which have to be made 
in order to obtain functions Q™ corresponding to a given F,,,. 


If M increases steadily from 0 to 1, A increases from — co to 0. Since the relation 
M <-> is a one-to-one correspondence, to every do, A, <0 there corresponds an 
My = M(A)) <1. Therefore, if A ¢(— oo, dy) then M (A)e(0, M,). For M < M,, 
the function 
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f= [64 (1 — M)*)> (kh + 1) M*[— (8% — 1) M*— 4 (3 — 2h) M2 + 16}, (6.11.2) 
may be approximated by a polynomial F,, of the (2” + 8)-th degree in M?, 
F,, (M) = (64)-1 (k + 1) M*[— (3k — 1) Mt — 4(3 — 2k) M2 + 16]- 


: 3 (v3) (— 1y ar (6.11.3) 


Only a finite number of powers of M? appears in (6.11.3). It will now be shown that M? 
can be developed in the uniformly convergent series 


co 


fe ae (6.11.4) 


v=1 


Yo 
Cee ee 


(+ 1 + (e—1)7'2 
with 4 <0. Instead of considering M2, it is well to introduce 


gS 1 (l= yh, (6.11.6) 


Since M? = 1 — (1 — s), it will suffice to determine the series for s. As one may 
easily verify, by simple calculations one obtains 


Be ata ieee aca 
eo ee ls k>1. (6.11.8) 


It is possible to prove that the image of the area of the circle |s| < 1 is a schlicht 
domain in the X-plane. That is, one can show (38, p. 58 and 59) that the boundary 
curve of the image of |s| < 1 is a curve which does not intersect itself. By classical 
theorems of the theory of functions, the domain bounded by this curve is schlicht. 
Clearly, it includes in its interior the domain |X| <1. Since the image of |s| < 1 
is schlicht and includes |X| <1, the inverse functions s=s(X) is regular in 
|X| <1 and can be expanded in |X| < 1 in the form of an infinite series 


Dens 


exp (2 A), (6.11.5) 


MOOI ODEs (6.11.9) 
v=1 
With some additional proofs it is possible to show that (6.11.9) yields the required 
approximation. 


12. Singularities in transonic flow. 


In one of the previous sections we indicated a procedure for constructing a mixed 
flow (partially supersonic) around an obstacle the boundary of which is a closed curve. 
If the boundary is prescribed, then probably, in many instances, no solution of the 
problem exists. Unfortunately, no definite results in this direction are known. 
However, investigations of problems of similar type seem to indicate that the above 
is the case. Therefore the problem arises of finding necessary and sufficient conditions 
in order that a mixed flow pattern may exist around an obstacle. In considering 
this problem, two cases have to be distinguished*: 

1. The behavior of the flow at infinity is completely prescribed (for example, a 
uniform flow with no vorticity), or, 

2. the flow at infinity possesses certain properties, but any solution will be 
considered as admissible if its behavior at infinity is of a certain type. 

The second problem is to determine when the solution is unique and, further, 
under what conditions it is “stable”. This situation suggests considering a problem 
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which is much simpler mathematically (but still quite difficult), and, in certain instances, 
can give the type of answer desired by a research engineer, namely, the question 
whether to a given hodograph there exists a flow possessing the prescribed behavior 
at infinity. Often it will suffice to know not whether a flow pattern exists around 
some prescribed (in the physical plane) profile, but only whether such a flow exists 
with prescribed (or approximately prescribed) velocity distribution. The whole 
investigation may be shifted to the hodograph plane, and the problem considered 
as formulated above. As an example, a flow will be considered in order to determine 
whether to a given hodograph there corresponds some mixed flow. While, as is known 
in the case of a subsonic flow, it is more convenient to operate in the (A, 6)-plane 
(pseudo-logarithmic plane), in the case of a mixed flow it is convenient to return 
to the (H, 6)-plane, since 1 (M) becomes imaginary for M>1. As A and H are 
connected by the relation dA = [1 (H ye dH, | (H) = (1 — M?) o-, it does not 
present any theoretical difficulty to replace the variable A by the variable H, in 
functions S;, (A, 6; 4», 99), & = 1, 2,3, which are assumed to represent a pseudo- 
vortex and two pseudo-doublets, respectively (see section VI. 5). Thus, functions 
are obtained which shall be denoted by 


T', (H, 0; Ho, 99) = Sx (A, 9; Ao, 90). (6.12.1) 


Let us now assume that there exists a function which at the point (Hp, 09) (which is 
the image of the point at infinity in the physical plane) behaves like 


A, T, (A, 63 ho, 69) + Ae Fo (H, 6; Ho, 0) + Ag Ts (A; 0; Ho, Oy) + & UT, 8); (6.8222) 


where R(H,6) is a regular function of H,6. Naturally, A; 1 = 1, 2,3, must 
satisfy the conditions discussed above, in order that in the physical plane a flow 
around a closed curve is obtained (6.8.9). Now, the question arises whether or not 
a function R (H, 0) can be found such that the expression (6.12.2) vanishes on the 
boundary line of the hodograph. Without going into too many details of a more 
advanced mathematical nature, let us note that one can find a necessary condition 
for the solution of this problem (for more details see +, pp. 40—44, or 14). 

Obviously, if the function VY (x, y) is to yield a flow pattern around a closed profile 
in the physical plane, it is necessary that 


(a) the streamline ¥ (x, y) = 0 (or a part of this streamline) yield a closed curve 
in the physical plane; 

(b) that every point of the physical plane be covered at most once (that is that 
the mapping onto the (xz, y)-plane be schlicht). It is remarked, however, that if the 
whole image in the (x, y)-plane is not schlicht, the flow pattern obtained may 
nevertheless be of physical interest since in some instance it is required to determines 
flow patterns covering only part of the (a, y)-plane or to use only a part of the flow 
pattern. 

The necessary and sufficient conditions in the subsonic case that the boundary 
in the hodograph plane when mapped onto the physical plane yield a closed curve 
in the latter plane have been established in previous chapters for flow patterns that 
are not necessarily symmetric. In the subsonic case, the determinant of the mapping 
from the hodograph to the physical plane 6(x, y)/0(u, v) is always positive. It is therefore 
sufficient and necessary, in order that the mapping be schlicht, that the boundary 
curve in the physical plane be schlicht, and that the mapping satisfy, at infinity, 
certain regularity conditions which will be stated below. In some cases the schlichtness 


1 S. Bergman: On tables for the determination of transonic flow patterns. Hans Reissner 
Anniversary Volume: Contributions to Applied Mechanics. J. W. Edwards (1949) 
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of the boundary curve follows from the symmetry of the image in the hodograph 
plane and a few other conditions“. The mapping must satisfy the following conditions 
at infinity : As a doubly covered circle about the point A (which represents the image 
of infinity in the hodograph plane) shrinks to a point, its image in the inverted physical 
plane x’ +i y' =(x +7 y)— must be schlicht and also shrink to a point. 


In generalizing these results to transonic flow patterns, one must be on guard 
against the following difficulties: 


(a) the functional determinant 6(w, v)/é(x, y) ist not necessarily positive, 


; (b) the mixed character of the equation of flow may give rise to hyperbolic 
singularities which do not appear in the subsonic case. 


(Eingegangen am 9. Mérz 1954.) 


Zur Bestimmung der Trigheitskoeffizienten unsymmetrischer 
Rotationsk6rper". 


Von W. Miiller, Miinchen. 
Mit 5 Textabbildungen. 


Zusammenfassung. In Fortsetzung einer Arbeit iiber die Integraldarstellung der Tragheits- 
koeffizienten eines Rotationskérpers wird nach Ableitung einiger allgemeiner Satze iiber die 
erzeugenden Quell- und Dipolsysteme die vollstandige Berechnung der Tragheitskoeffizienten 
einer besonderen Formengruppe mit Hilfe von Ellipsoid-Koordinaten durchgefiihrt, die eine 
geschlossene Lésung der vorkommenden Integrale erméglichen. Der Vergleich mit dem Ellipsoid 
fuihrt zu mehreren, auch fiir andere Formen verwendbaren Gesetzen. 


Summary. In a continuation of a paper on the integral representation of the coefficient 
of inertia of a solid of revolution, some general theorems are derived on the generating expanding 
and dipole systems and the complete calculation of the coefficients of inertia is shown for a special 
group of forms with the aid of ellipsoid coordinates, which enable a closed solution to the integrals 
occurring. The comparison with the ellipsoid leads to several laws applicable also to other forms. 


Résumé. Comme suite & un travail sur la pose de lVintégrale des coefficients d’inertie d’un 
solide de révolution, et aprés la démonstration de quelques théorémes généraux sur les systémes 
générateurs de source et bipolaire, les coefficients d’inertie d’un groupe de formes particulier 
sont entiérement calculés a l’aide de coordonnées ellipsoidales qui rendent possible la solution 


totale des intégrales qui se présentent. La comparaison avec l’ellipsoide méne & plusieurs lois, 
également applicables 4 d’autres formes. 


. Inhalt: 
. Kinleitung. 
. Ellipsoidkoordinaten. 

. Darstellung von geschlossenen Rotationskérpern. 

Anwendung auf den Fall eines unsymmetrischen Rotationsk6rpers. 
Berechnung des Triagheitskoeffizienten k der Langsbewegung. 

Allgemeine Satze tiber die Potentialfunktion der Querstrémung. 

. Bestimmung einer weiteren Naherung fiir das Potential der Querstrémung. 
. Ermittlung des Triagheitskoeffizienten k, der Querstrémung. 

Ergebnisse der Rechnung. 

10. Beziehung der Koeffizienten zu den Geschwindigkeiten. 

11. SchluBbemerkung. 


I. Einleitung. 


In einer vorbereitenden Arbeit? iiber die Hydrodynamik des Rotationskorpers, 
im besonderen iiber die Bestimmung des von den Luftkraften herriihrenden instabilen 


1 Die vorliegende Untersuchung geht zuriick auf eine frithere Zusammenarbeit mit der Firma 


Messerschmitt. 
2 Osterr. Ingenieur-Arch. 8, 171 (1954). 


Ingenieur-Archiv VIII, 4. . 18 
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Momentes, das in erster Linie von den auf die Lings- und Querrichtung bezogenen 
Trigheitskoeffizienten abhingt, habe ich zur Verdeutlichung der theoretischen Er- 
gebnisse das Beispiel des Ellipsoids nach der Quellsenkenmethode entwickelt und 
dabei die Ellipsoidkoordinaten zugrunde gelegt. Im folgenden soll diese Methode 
einschlieBlich des Gebrauches der Ellipsoidkoordinaten auf andere Korper ausgedehnt 
werden, die eine unsymmetrische Gestalt haben, das hei&t durch eine verschiedene 
Ausbildung von Bug und Heck gekennzeichnet sind. 


Nach einigen einleitenden Bemerkungen iiber diese Koordinatendarstellung werden 
mehrere allgemeine Sitze iiber Quell- und Dipolsysteme entwickelt, die fiir den Aufbau 
von praktisch brauchbaren Rotationskérpern in Frage kommen. Vor allem soll hier 
zum erstenmal die vollstandige Berechnung der Tragheitskoeffizienten einer besonderen 
Formengruppe durchgefiihrt werden, die insofern bemerkenswert ist, als die auf- 
tretenden Integrale in vollstindig geschlossener Form lésbar sind. Der Vergleich 
der Ergebnisse mit den Werten des Ellipsoids gleicher Lange und gleichen Volumens 
liefert einige bisher nicht bekannte GesetzmaBigkeiten, die als Grundlage einer Unter- 
suchung zusammengesetzter Falle dienen kénnen. 


II. Ellipsoidkoordinaten. 


Unter Zugrundelegung rotationssymmetrischer Anordnung bezeichnen wir die 
Koordinaten in einer Meridianebene, das heiBt die axiale und radiale Koordinaten 
in einer Meridianebene, mit x, und wenn wir dann von einer in der Achse liegenden, 
von x= —c bis « =c reichenden Grundstrecke ausgehen, deren Endpunkte C, 
und C, die Brennpunkte der zugeordneten Meridianellipse sind, und die Radien- 
vektoren von diesen Punkten nach einem beliebigen Feldpunkt mit 7, und r, bezeichnen, 
so sind die Ellipsoidkoordinaten definiert durch die Beziehungen 


und die Radien 7, und r, sowie die Koordinaten x, 7 driicken sich dann in folgender 
Weise durch ¢ und uw aus: 


m=e(S+u); m=eb—p); e=celu; n=cV@—-HO—p. (2) 
Man sieht zunachst, daB ¢ = konst. = ¢) bzw. w = konst. = wy einer Ellipse bzw. 
Hyperbel in der Meridianebene mit den Gleichungen 


6 


22 2 ae 2 
By H Ge y 


ete? ete) ue Da) (3) 


und den Halbachsen a = cl), 6b =cVC,2— 1 baw. a’ =Cm, 6b’ =cV1 —wue ent- 
sprechen, die bei Drehung um die Achse ein Rotationsellipsoid bzw. Rotations- 
hyperboloid ergeben. Wenn ¢, bzw. wo alle méglichen Werte in den Bereichen 
1<¢€ < oo; —1<yp <1 durchlaufen, so erhalt man in jeder Meridianebene ein 
konfokales Netz von Ellipsen und Hyperbeln, das als Koordinationsnetz verwendet 
werden kann. Statt ¢ und w kann man auch manchmal bipolare Koordinaten mit 
Hilfe der Transformation 


i UAE 


C=Col a; = cose (4) 
einfiihren, womit die rechtwinkeligen Koordinaten die einfache Form 
x=cojxcosB; 7=cCinasinp (5) 
erhalten. Diese Darstellung hat den Vorteil, da&B man gelegentlich die komplexe 
Schreibweise einfiihren kann wegen . 


z=" +in = cOoj (a + +8), 
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wodurch manche Formeln eine einfache Gestalt erhalten. Wir wollen jetzt einige 
einfache Funktionen der Koordinaten ¢ und y» betrachten und die entsprechenden 
Bildkurven in der Meridianebene verfolgen. Setzen wir etwa Proportionalitiét von ¢ 
und w voraus, also eine Gleichung 


das heifit eine durch den Nullpunkt gehende Gerade in der ¢ u-Ebene, so folgt mit 
Hilfe von (2) 


hola i d Eas 
2 rey 4 pl t#=aele——) ") 
das heiBt ein Kreis mit der Gleichungsform 

(@ — 2) +1? =a" 


hat, dessen Mittelpunkt in # = a = : C (4 + ) liegt und dessen Halbmesser 


. 33 1 il ; : ‘ 
die GroBe a = > ¢ (9 —- =| hat. Diese Bestimmungsgré8en gentigen der Beziehung 


ty — a = ¢. (8) 


Man findet also einen Punkt des Kreises als Schnitt des Halbkreises iiber der Strecke 
OM (= xo) mit dem Kreis um O mit dem Radius c. Geht man umgekehrt von dem 
Kreis aus 


(x — 2%)? + Pp = a, 
wo jetzt x) und a beliebig sind, so erhalt man als Bildkurve 
e7 (f? + u?) —2caCu=—a+c — 2,2. (9) 
Diese Gleichung stellt, wie man leicht feststellt, eime Ellipse oder Hyperbel dar, je 


nachdem ¢= wx, wird, deren Hauptachsen mit den Grundachsen Winkel von 45° 
bzw. 135° bilden. Bemerkenswert sind folgende Sonderfalle: 


a) Wenn 


to — a= c? 
ist, erhalten wir wieder den oben angegebenen Fall mit der einfachen Gleichung: 
% +a 
C= 4 = gb. 


b) Ist a) = ¢, fallt also der Mittelpunkt M des Kreises in den rechten Endpunkt C, 
der Grundstrecke, so entsteht die einfache Beziehung 


C—w= 45 (10) 


als Gleichung zweier Geraden durch die Punkte € = + = , & = 0, die mit der ¢-Achse 


Winkel von 45° bilden. 
¢) a = —c gibt in der a 7-Ebene einen Kreis um den linken Endpunkt C, der 


Grundstrecke und in der ¢ u-Ebene zwei Gerade durch die Punkte ¢ = + = ti =a, 
die mit der ¢-Achse Winkel von 135° bilden. 


d) x, = 0, dem Kreise 
a? + 72 = a? 
um den Nullpunkt entsprechend, gibt in der Bildebene den Kreis 


4+ ye=144. (11) 
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Diese Bemerkungen sind fiir die Konstruktion gewisser Grenzkurven von Rotations- 
flichen von Nutzen. 


III. Darstellung von geschlossenen Rotationskérpern. 


Um geschlossene Rotationskérper zu gewinnen, zu denen also auch geschlossene 
Meridiankurven entsprechen, geht man am besten aus von einer Verteilung von 
Quellen und Senken lings der Achse, deren Gesamtergiebigkeit gleich Null ist. 
Hine einzelne Punktquelle im Punkt x = a der Ergiebigkeit q hat das Potential 

q q é 
= — = 12 
@ 4nr 4a (a — a)? + 4?” (12) 
aus dem sich dann, wie bekannt, die Geschwindigkeiten durch Differentiation ab- 
leiten. Bei einer beliebigen, von der « = —c bis « = + ¢ definierten Verteilungs- 
funktion f(«) erhalten wir daher das Potential 


+e 
Dee eas (13) 


wo A eine Konstante bedeutet, die von der Quell- bzw. Senkenstairke und daher 
indirekt auch von den Abmessungen des zugeordneten Rotationskorpers abhangt, 
dessen Axialgeschwindigkeit mit v, bezeichnet ist. Die SchlieBungsbedingung oder 
die Bedingung dafiir, daf der Rotationskorper geschlossen ist, lautet dann 

+e 


| f (a) da = 0. (14) 


Es 1aB8t sich nun zeigen, da8 die Quellverteilung ersetzt werden kann durch eine 
gleichwertige Verteilung von Dipolen, die nach der x-Richtung orientiert sind: die 
Beziehung beider Funktionen ist derart, daB die negative Ableitung der Verteilungs- 
funktion g der Dipole nach der axialen Koordinate die Verteilungsfunktion f der 
gleichwertigen Quellen ergibt. Setzen wir naémlich 


f(a) = — (15) 


@—A Vy | _ S 


so wird das Potential 


Durch partielle Integration ergibt sich aber 


®= Av, 2"! + Ao, | al dr 


r J) 7? 


(16) 
g (x) |® ( B do 
a Agel Og cee 
Da das Potential eines einzelnen Dipols im Punkte x = « die Form 
m cos? 
@ Prey Tee (17) 


hat, so ergibt sich, daf der letzte Ausdruck in (16) identisch ist mit dem Potential 
einer Dipolverteilung mit der Verteilungsfunktion g. Wenn also im besonderen 
9: = 92 = 9 wird, also die Werte von g in den Endpunkten der Grundstrecke ver- 
schwinden, so ist ® nur von den Dipolen abhingig. Diese Bedingung ist aber ver- 
traglich mit der Bedingung (14), die in unserem Falle lautet: 

- a 
| #(a)da = — | Fda = — [9 (a)],. (18) 


1 a 
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Diese Bedingung ist z. B. erfiillt in dem Beispiel f= 4, ¢ = 
sich zeigen lit, zum Ellipsoid fiihrt. 

Der Satz la8t sich verallgemeinern auf den Fall, da& wir etwa in den Endpunkten 
der Grundstrecke eine Einzelquelle und Einzelsenke haben neben einer stetigen Senken- 


verteilung langs der Grundstrecke. Wir erhalten dann, da wir jetzt die positive Ab- 
leitung zu nehmen haben, 


(c? — «?), das, wie 


: qo q dg (x) 
he 4ary iar | A », | r (19) 
on 1 (Ge - % g g " g (x) cos bdo 
4x Gee 2) + Av, (2 —2)_ do,(! rae ee 
Nun ist nach der SchlieSungsbedingung 
r 1 
A\ f(a) da = A (92 — 91) = Gea (Ge — Mi). (20) 
1 
Wenn also die Bedingungen 
ee ate Re per oa S 
A 92= Gano, ADT can can, 
erfiillt sind, so 1iBt sich in der Tat das Potential in der Form 
2 
" g (x) cos Fda 
1 
darstellen, das heiBt von einer Verteilung von Dipolen allein ableiten. Wenn im 
besonderen die Einzelsenke in x = — c verschwindet, so wird mit q, = q 
Ag=z; 91 = 9. (21a) 


Dieser Fall ist im besonderen geeignet zum ersten Aufbau von einfachen unsymmetri- 
schen Rotationskérpern und soll daher den Rechnungen in den folgenden Abschnitten 
zugrunde gelegt werden. 

Zuvor noch einige Bemerkungen zu der Potentialfunktion und der zugeordneten 
Stromfunktion Y. Die Potentialfunktion gibt durch ihre Ableitungen die absoluten 
Geschwindigkeiten in der Umgebung des axial bewegten Rotationskérpers. Das 
Potential fiir die Relativstromung lautet dann, wenn v, die Bewegungsgeschwindigkeit 
bezeichnet, 

OD’ = —,x44+ @. (23) 


Spaltet man den Faktor v, ab, so erhalt man bei Zugrundelegung der Form (13) die 
fir die Folge maBgebende Quellfunktion 


2 
@ > f (a) da 
P=, = i aS (24) 
der die Quellfunktion der Relativbewegung 
By alge © (25) 


entspricht. Die Meridiangrenze des so bestimmten Rotationskorpers ergibt sich aus 
der Stromfunktion Y, die aus der Potentialfunktion durch eine einfache Integration 
gewonnen werden kann. Allgemein hat man 


ob rr oF ( )n dad 
C— 1 ez 
fi eee ae oh (26) 

03 


w= |\F nig =—Av e 


0 
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Wenn man Ellipsoidkoordinaten einfiihrt, so hingen ® und WY durch die beiden 


Gleichungen 
ov 
Ou 


BS op aw a® 
= —¢c(@—1) a FE i) aie (27) 


zusammen. 


IV. Anwendung auf den Fall eines unsymmetrischen Rotationsk6rpers. 


Wir setzen die Verbindung einer in x = c gelegenen Punktquelle mit einer Strecken- 
senkung gleicher Ergiebigkeit voraus. Das Potential der Punktquelle in Ellipsoid- 
koordinaten wird 


ep ae en Ss he 28 
a= 4n%, 42¢ (€—u)° (28) 
Fiir die von x = —c bis x =c reichende Streckensenke haben wir mit f(«) = 1 


ota 
= Av,\—, 
i 


wo die Konstante A aus dem Wert der Gesamtergiebigkeit der Strecke, das heiBt 
aus der Gleichung 


Zu 
ig Se (29) 


sich ergibt. Zur Ausfiihrung der Integration verwenden wir voriibergehend als Variable 
den Winkel & zwischen dem Radiusvektor nach dem Achsenpunkt « = « und der 
x-Achse. Dann wird 


ctg } = ak xr—a=neted, 


ad e 
hie ee 71 — sin %. 


sin? } ? r 


Das unbestimmte Integral geht dann iiber in 


r da r «on ao r a u 1 — cos @ 
| nes | sink Or | snd 2 Ig 1+ cos 8° (30) 


Gehen wir bei der HKinsetzung der Grenzen zu Ellipsoidkoordinaten iiber, so wird 
fiir die Endpunkte der Grundstrecke 


x—e Cu—l te+1 
cos 0, = = ; 9, = Se 
2 = Sg ee vy Sear ey (31) 
also 
1 1—cos 0]? €é+1 
2 1 + cos r= ar ee (32) 


Wir erhalten daher als Gesamtpotential fiir das System Punktquelle + Strecken- 
senke 


q q ¢+1 1 1 C 
= =-7o | -Fske Hl 6s) 


anc(¢—m) | @ne SE—1 4xc|F—u 1|° 


Nach dem angegebenen allgemeinen Satz erhilt man denselben Ausdruck, wenn 
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von einer von « — — ¢ bis x = ¢ reichenden und linear zunehmenden Dipolverteilung 
(kurz von einem Dipolkeil) ausgeht. Das kann direkt bestatigt werden. Setzen wir also 
g(a)=c+a, 


so ergibt sich als Potential des Dipolkeils 
P=— Av, [eee ate*) SRE 
Durch Auflésung der dabei auftretenden Integrale erhalt man 


2 


; 1 — cos & 2 + 3 2 
P= Avy, lg a face ee A ¥, [cos Clo Av, = [sin 7] . 


Mit Benutzung der Formeln 


2 2¢(1— z 2n V(e—1) (1 — #2) 
[eos #} — — ihe ws. [sn ? 2m fein oy — =A = #2) 


a re LS a 
(34) 


sowie der Gl. (32) erhalt man in der Tat den obenstehenden Ausdruck fiir das 
Potential. 


Wenn man die Konstante g durch die Geschwindigkeit v, und die Abmessungen 
des zugeordneten Rotationskérpers ausgedriickt denkt, so haben wir damit gleich- 
zeitig das Potential fir die absolute Stromung in der Umgebung des axial bewegten 
Rotationskorpers. Spaltet man v, ab und setzt 


62 — —2_ 


TU,” 
so erhalt man die Funktion F, in der Form 


6? 1 +1 


Mies ee ae lg eer (35) 


F', hat die Dimension einer Lange. Die GroBe c gibt ein Maf fiir die Langenabmessung 
und die GréBe b ein MaB fiir die Quer- oder Breitenabmessung des Ko6rpers. Die 
Gleichung der Meridianbegrenzung erhalt man aus der Stromfunktion WY, die sich 
mit Hilfe der Beziehungen (27) und (33) berechnen 148t. Man erhalt durch eine 
geringe Rechnung 


[ c2— ] 1 eae 
Y= - [é bf (eae SEO ie ae (36) 
daher fiir die Relativbewegung 
yee - vin + Ww. (37) 


Daraus leitet sich dann die Gieichung des Meridiananschnittes C in Ellipsoidkoordinaten 
ab in der Form 


u2 


1 1 | — 
ET CR ay iee a (38) 


Mit Einsetzung des Wertes fiir 7? ergibt sich daraus 


b? b2 
biiak = Fenteedy) ehins Hy PME Be 
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Diese Gleichung gestattet in sehr einfacher Weise, die Werte von « aus den ¢-Werten 
und damit [vgl. (4) und (5)] die Schnittkurve punktweise zu berechnen, waihrend 
man bei rechtwinkeligen oder Polarkoordinaten irrationale Ausdriicke und quadratische 
Gleichungen in Kauf nehmen mu8. In der Tabelle 1 haben wir eine Reihe von Koordi- 
naten ¢, « mit den zugehérigen rechtwinkeligen Koordinaten, und zwar fir sechs 


b ; 
Werte von -~ wiedergegeben. 


Tabelle 1. Berechnung des Meridianschnittes des Rotationskérpers. 


vw =01 aa =02 
€ c 
x n x n 
(3 B ea ak c & | cs "el 
1°00125 —1°0 | — 1:00125 0 1°00497 —l — 1°00497 0 
1°0015 — 0°665 | — 0°6655 0°041 1°005 — 0°989 — 0°994 0°015 
1002 —0°248 | — 0:2485 0°0612 1°006 — 0°655 — 0°659 0°083 
1°004 H 0°380 0°3815 0°0829 1°008 — 0°2373 — 0°239 0°123 
1°01 0°761 0°7689 0°0920 1°01 0°015 0°0152 0°142 
1°02 0°896 0°9141 0°0892 1°02 0°525 0°536 O171 
1°04 0°9787 1°0178 0°0587 1°04 0°794 0°826 0°174 
1°0494 1:0 1°0494 0 1°06 0°898 0°952 0°154 
— — — — 1°08 0°960 1°037 0°115 
a = — — 1°09 0°984 1°076 | 0°078 
pee = ws = 1:0976 1:0 1:0976 0 
2 Lo 2 =06 
x n x n 
: Rite eee a : r; 2 
1°0196 —]l — 1°0196 0 1°04315 —] — 1°04315 0 
1°02 — 0°96 — 0°98 0°056 1°045 — 0°9115 — 0°9525 0°1247 
1°025 — 0°555 — 0°57 0°187 1°05 | ——0°706 — 0°7413 0°1247 
1°03 — 0°284 | —0°293 0°236 1°08 — 0°0017 —0°0018 | 0°2263 
1°05 0°27 0°284 0,308 1°10 0°243 0°2673 | 0°408 
1°08 0°599 0°648 0°326 115 0°592 0°681 0°444 
1°10 0°719 0°792 0,320 1°20 0°791 0°949 0°458 
1°15 0°903 1:04 0°244 27 0°976 1°24 0°406 
118 0°976 1°152 0°136 1°2809 1 1°2809 0°169 
1719115 ] 1°1912 ) — — — 0 
| 
b f b 
c Mu =~ a c u = a 
1:0744 —] — 1°0744 0 Lo —]1 — 11121 0 
1°08 — 0°843 — 0°910 0°2195 1:12 | — 0°842 — 0°943 0°273 
1:10 | = O42 aa = O-4Ge 0°416 1°14 O58 — 0°605 0°464 
19 5) 0°158 0°1817 0°561 1°16 — 0°287 — 0°333 0°564 
1°16 0°234 0°2714 0°571 118 — 0°09 ase 0°1062 0°622 
1°18 0°3645 0°430 0°583 1°20 0°063 0°756 0°664 
1:20 0°472 0°566 0°585 1°25 0°362 0°453 0°699 
sO) 0°836 1087 0°457 1°30 0°576 0°749 0°678 
1°3676 1:0 13676 0 1°40 0°88 1°23 0°466 
-- -- — = 1°4515 1 1°4515 0 


Die beiden Extremwerte ¢, und ¢,, den beiden Scheitelpunkten mit w = + 1 
entsprechend, ergeben sich aus je einer Gleichung dritter Ordnung, die am besten 
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mit Hilfe von Kreis- oder hyperbolischen Funktionen aufgelést werden kann‘. Aus (2) 
und (39) bestimmt sich dann die Lange des Potentialkérpers zu 


1= 2¢(0, +¢,). 


Die groBte Breite, die von 6 abhangt und angenihert 6 proportional ist, liegt 
um so weiter nach vorne, je kleiner b ist, und riickt mit wachsendem 6 mehr und 


mehr in die Mitte. Die Breitenabmessung 
selbst ist stets kleiner als die GréBe 6, die 
den Grenzradius des unendlich langen 
Halbkorpers darstellt, der entsteht, wenn 
der eine Endpunkt C, der Strecke ins 
Unendliche riickt, wie man leicht durch 
Grenziibergang zeigen kann (Abb. 1). 
Wir haben die Extremwerte ¢, und ¢ , 
sowie die Lange / und den gré8ten Durch- 
messer d fiir spateren Bedarf fiir einige 
Werte von b bei festgehaltenem c zusam- 


mengestellt (Tabelle 2), ferner d und / in 
Abhangigkeit von 2 graphisch aufge- 
zeichnet (Abb. 2). 


Abb. 1. Meridianschnitte von Rotationskérpern 
erzeugt aus Punktquelle + Streckensenke oder 
Dipolkeil. 


Fiir die Berechnung des Triagheitskoeffizienten k, der Lingsbewegung interessiert 
vor allem der Verlauf von F, lings des Meridianabschnittes des Rotationskorpers. 
Wenn wir € — w aus (39) in den Ausdruck fiir F, einsetzen, erhalten wir 


Be Nae 
Le 


1 ¢+1 


(40) 


Der Verlauf dieser Funktion lings C ist fiir den Fall = = 0°2 in der Abb. 4 wieder- 


gegeben. 


4 Kir die obere Grenze hat man z. B. 


62 


Cos ee a CS = 0 
1 ; 
Mit ¢.=t+ > ergibt sich die Normalform 
— 3 pit—2q=0, 
b? 8 
ee gee ce OF 
Da hier g? < p® ist, so kann man setzen 
soa cae - b? 
PORE PSE ape pow 4@ 8 1 3°375 Le il 


Dann ist die fiir uns in Betracht kommende Losung 


= u 4 U 
t=2)p C08 5, re SO ae 


Fiir die untere Grenze erhalt man in ahnlicher Weise 


— 


4 8 b2 , 3 
Pars A omelet) TPs 
q gee ia ae or 
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Tabelle 2. Berechnung der Groen ¢,, 6, und J, K, L. 
I 


b : Qa Bars pie 
2 =01 ~ = 02 == "4 
ear ah pene by e di te 
Eee ee eee 
| 
Le 4 at oth Sepa A 100125 | 1:0494 100497 10976 101961 1°19115 
FARE IRE © 0:0025 0°10124 0:009965 | 0°20472 0°0396 0°41884 
5 fs oe 400°562 10°3654 | 100°748 5:3614 25°747 2°8439 
Ar Ty = ee 369 1°8626 3°00 1°5338 2°3173 1°2195 
OnE en 26946 0'9546 2:0149 0°6835 1°3627 0°4526 
OEE Serre 396'872 85028 97°748 3°8284 23°417 16244 
Eee NS eeoeery 0:00968 0:03763 0°1430 
Cc 
et oer 0:0003196 0:002378 0°01705 
Cc 
Uh Ae aa 4 12904 12207 
(ox 
206 sO eit 2 iG 
c c ce 
cy Co at Ge oy | fo 
hannlela ahead 1:04315 | 12809 10744 13676 11121 1°4515 
al Rig, maa 0:0881 0°6408 0-15433 | 0°87033 0°2367 1°1069 
a. oMarene 11°840 1:9989 69034. | 15714 4°698 13113 
Ur Iq = werent! 0°9586 0°7808 1664 0°9315 1468 0°8460 
Oe > eek 1012 0°3416 0°7878  0°2739 0°6326 0°22797 
ted | A Rane 9°9110 0°9515 52977 | 0°6399 3-230 0°4653 
! ~ = SSS 
= ae oe 0°3084 052964 0°80792 
1 
eS eee 0°051384 0°11034 0°19456 
1 
EO ny eV 1170 11358 11030 


V. Berechnung des Tragheitskoeffizienten k, der Lingsbewegung. 


Das instabile Moment, das entsteht, wenn die Bewegungsrichtung um einen 
Winkel 6 von der Langsrichtung des Rotationskérpers abweicht, ist bekanntlich 
von den Komponenten des scheinbaren Volumens abhingig, die man erhilt, wenn 
man das wirkliche Volumen V des Korpers mit den Tragheitskoeffizienten k, und k, 
fiir die Langs- und Querrichtung multipliziert. Beide Koeffizienten lassen sich, wie 
ich im ersten Bericht ausgeftihrt habe, durch einfache Integrale lings der Meridian- 
grenze des Korpers berechnen. Fiir k, hat man z. B. den Ausdruck 


rn 


n \ F, dn? \ Fy dy? K 
x | 2 dx i | 72 da hae 


is (41) 
wo z J = V das Volumen des Rotationskorpers bezeichnet. Es soll nun gezeigt werden, 
da k, durch einen geschlossenen Ausdruck dargestellt werden kann, in dem nur 
die beiden Scheitelwerte ¢, und ¢, auftreten. Wir fiihren in dem Integral die un- 
abhangige Variable ¢ ein, Aus 


PC ai) ee Be) 
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erhalten wir 
1 
Cf OAS. Ic (1 — pw?) — met (¢2 — 1)| dt. 
Nun folgt aus der Gl. (39) der Meridiankurve 
du ck 
ey a ey 

Damit wird 

te 2 bo 

Fe eG fe otek (ees ah) Se ce an 
Mit Einsetzung der Werte 


Ce ee | bec é aoe Eo peer eee ee 
f DER AC” 1en ere ee cad) 


erhalt man daraus nach Vereinfachung 


b? bre 
dnt = — 208 |2¢(@ — 1) — 9 e2 rear e- 
Daher wird 
1 402 1 
Fy orf = 08/5 o( — 18 —(@ — 1) —o(@ — les tt 
be C ¢+.1 C Cet 
Peers eee Pipcies (aac een iy) 
Benutzt man nun folgende Integralformeln: 
1 1 1 
f@-pydg= ees; \e-ynat=4ce@—-3); |\a ta =Zle(@—v, 
1 lesa 
[e@—le ey df= 5 (@— lest, + FZ 6G — 3), 
(43) 
1 
Jig =F at = cig ++ + 1g ( — 1), 
¢ Cart 1 Tea! 1 ¢ Li ped at 
lesa Bea“ = — ey BS tele 3 8h 
so ergibt sich fiir K folgender Ausdruck: 
1 1 gue 
i es oe sli ioe G I)lg--— 4 
c+ - 1 bt Peet lpg eal pe ey 
rie ra bles even yer) ee eas — gig ays (44) 


Der Ausdruck kann weiter vereinfacht werden, wenn man beriicksichtigt, daB fiir 
die obere und untere Grenze, das hei&t die Scheitelwerte ¢, und ¢,, folgende Beziehungen 


bestehen : 
, b? ae b? 
ieee Praag hy, cans BEG aes 


(45) 


Dann erhalt man 


2 ib gfe & Owaiae 
K=e|t—1 +g 568-O+g eae eta), 46) 


Wenn man schlieBlich noch mit Hilfe von (39) c? durch b? ausdriickt, so kann man 
die Werte fiir die obere und untere Grenze auf die Form bringen: 


Smians Cc. [6G C*) ae (Ce ph laren #44. =, (47) 


K,-e(e SE -M-20-M +a ar Oh, 
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wo zur Abkiirzung 


4 too ee we 
gesetzt ist. Dann wird 


K —— Ky a ne 
Ebenso 1a48t sich das Volumen in geschlossener Form darstellen. Man erhalt zunachst 
$e 
V Wiawoe DEG Poe (G aA) F 
a fee =|) oo + aie se (— 1) ares 
oy 
b2 
Sa (a igs. 


Mit Benutzung der oben angefiihrten Integrale und der weiteren Formel 


PoE? ec de ga cle ee ee 
| (cee ee 1 2 (@2— 1)?’ 
ergibt sich 
1 62 1 bt 1 bo ¢ p6 2 
= Co 3 Cc" ae 1); 2 cf (C21) 3 (C4 1) | 16 cé (2. 1)2 39 6 Qi], (48) 


oder mit (39) 


Rael es¢-yTe-»2-d-gw Gb. 


(49) 
{1 Be i 8 
J=elag C+ DT(C—2) +01 —arn Gh. 
Damit wird der Tragheitskoeffizient k, 
po Fee re 
i = ere aie (50) 


Auf die Abhangigkeit dieser Koeffizienten von e und das Verhaltnis . komme 


ich an spaterer Stelle zuriick, nachdem die Berechnung des Tragheitskoeffizienten k, 
durchgefiihrt ist. 


VI. Allgemeine Satze tiber die Potentialfunktion der Querstrémung. 


Die Auffindung der Querstr6mung bietet zunachst einige Schwierigkeiten. Da 
aber die hier zu betrachtenden Rotationsformen mit dem Ellipsoid vergleichbar sind, 
kann man von folgenden allgemeinen Uberlegungen ausgehen. Wir haben oben ge- 
sehen, daB die Langsstromung aus einer Dipolverteilung gewonnen werden kann, 
bei der die Achsen der Dipole in die x#-Richtung gestellt sind. Bei einem einzelnen 
Dipol erhalt man bekanntlich als Stromgrenze eine Kugel und die Drehung des 
Dipols ohne Anderung seiner Starke hat eine Anderung der Bewegungsrichtung der 
gleichbleibenden Kugel zur Folge. Die Langsstr6mung um ein Ellipsoid gewinnen 
wir aus einem Quell- und Senkenkeil, die durch die Verteilungsfunktion f («),= « 
(—c <a Sc) zusammengefaBt sind oder einer damit gleichwertigen parabolischen. 


Dipolverteilung mit der Verteilungsfunktion g («) = = (?? — «?). Die Querstr6mung 


des Ellipsoids infolge Bewegung in der y-Richtung wird auch in diesem Falle da- 
durch erhalten, da die Achsen der Dipole bei gleichbleibender Verteilungsfunktion 
in die y-Richtung gedreht und gleichzeitig die Konstante, das heiBt die Gesamt- 
starke, der Randbedingung entsprechend, abgeindert wird. Dadurch wird die Quell- 
verteilung / (x) fiir die Langsstrémung der Ableitung der Dipolverteilung h(«) = Ag («) 
fiir die Querstrémung nach der axialen Koordinate proportional oder beide Funktionen 
f(x) und h (x) verhalten sich im Falle des Ellipsoids, von einem Proportionalitats- 
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faktor abgesehen, wie die Querkraft- und Biegungsmomentenverteilung bei 
emem senkrecht zu einer Achse belasteten linearen Balken. Obwohl dieser Satz 
streng genommen nicht fiir jede Verteilungsfunktion giiltig ist, so kann man den 
Vergleich doch benutzen, um sich eine erste Anschauung von der gesuchten Dipol- 
verteilung g(x) zu bilden, die denselben Rotationskérper als Grenze ergibt wie die 
Ausgangsverteilung f(«) der Quellen und Senken. 


Insbesondere li8t sich die Querstrémung fiir den aus Quellpunkt und Senken- 
strecke abgeleiteten Fall des ovalartigen Rotationskérpers in erster, allerdings nur 
roher Annaherung aus einem Dipolkeil, das hei8t aus einer linearen Verteilungs- 
funktion von Dipolen ableiten, die in die y-Richtung gestellt sind, wenn man nach- 
traglich die noch unbestimmte Konstante etwa durch die Forderung festlegt, daB 
der Staupunkt der Querstrémung auf die Schnittkurve des Kérpers mit der Meridian- 


ebene » = 0 fallt. Fiir die Querstrémung eines Dipolkeils haben wir aber die 
Potentialfunktion 


2 2 

; in Oda = ‘ 
| (c + a) sin bdo x + c—yn ctg Linda (51) 
1 


®, = — Bv,cos@ / = — Br, cosy | jj 
ui 


wenn ¢ den Polarwinkel in der Parallelkreisebene gegen die y-Achse bedeutet. Wenn 
man die Integration ausfiihrt, so erhalt man 


®, = — Bv,{*~** [cos 6], + [sin 8] | (51a) 
7] uf iJ 
oder nach Hinsetzung der Ellipsoidkoordinaten und nach Abspaltung von »v, 
RH! yi a Vi@— 1) i—-) ‘ 
ah we eo ee CSN roe (om 
Nun ist aber auf der Meridiankurve des Rotationskorpers 
b2 
Cowie a= = 3. 
Man erhalt also langs C 
268 ae Be , 
(F,) = —2B=--V(@—-)0-)=—4 >Re. (53) 


F, nimmt also auf der Meridiankurve C Werte an, die mit den negativen Ordinaten 7 
verhialtnisgleich sind. 


Der entsprechend verallgemeinerte Satz kann dahin formuliert werden, da die 
potentielle Funktion F, eines quergestellten Dipolsystems mit der Verteilungs- 
funktion g (x) in dem Meridianschnitt y = 0 eines Rotationskérpers Werte annimmt, 
die den negativen Abstinden 7 von der Achse verhiltnisgleich werden, wenn die 
den Kérper (bzw. die Lingsstrémung um den K6rper) erzeugende Quell- bzw. Senken- 


: d 3 
verteilung / («) als negative oder positive Ableitung f(«) =F oo der Dipolver- 


teilung darstellbar ist. 
Wir nehmen etwa eine Quelle in C, und eine Senkenverteilung tiber der Grund- 
strecke (, C, an, dann lautet nach (13), (26), (27), (28) die Stromfunktion der Absolut- 
strémung nach Division mit v, 
1 q 
Sin ail _ TG ED +A 


Da der Parallelstrémung in Richtung der negativen x-Achse die Funktion 


(54) 


Vi 


fg) or eakde 


1 
YY — Sent 
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entspricht, so erhalten wir als Stromfunktion der Relativstromung 


1 A 1 q f (x) (~@— a) dx 54 
= Sop iia cos By +A\? - ‘ (54a) 
Damit ergibt sich aber als Gleichung des Rotationskorpers 
f (a) (w— a) do 
ie 


j= — Tyo 008 8, +24 | 

oder nach Division mit 7 

( f(a) (w©— a) da q 

n=2A\- cae Fav, 7 008 De. (55) 
dg (x) 


Setzt man jetzt f(«) = —q7,_> 80 kommt durch partielle Integration 


9A 2 - 9 (x) iT (a — a)? q 
n=" fg (a) cosy, —2.4 | 2 | 4 Bo lda — a 008 
oder wegen 7? = 7? + (x — «)? 

2A Due: a Slate 
n =~ [g (x) 008 8], — ga5 a 008 Os +24|# 


Wenn g, = 0 wird, so folgt aus der scuisiesea ete 


g (x) sin & eB 
7 


g 
ye 


42, 
die einfache Gleichung 


- g (x) sin 0 da 


n= 2A \ = . (56) 


Das Integral auf der rechten 
Seite stimmt aber formal bis 
auf eine Konstante tiberein mit 
der potentiellen Funktion einer 
Verteilung g («) von Dipolen, 
die in die y-Richtung gestellt 
sind, also mit der Funktion 
g (x) sin da 


72 


? 


Abb. 2. Verlauf der Funktionen #, und F, langs der #, — — B cos@ \ 
Meridiankurve C. . 


(57) 
wenn man cosg = 1, also y= 0 setzt. Der damit eingefiihrte Hilfskérper, ent- 
sprechend der Quellverteilung f («) = — weicht um so weniger von dem urspriing- 


lichen Rotationskorper ab, je weniger dieser vom Ellipsoid sich unterscheidet. Man 
kann also sagen, daB bei geringer Abweichung des Rotationskérpers vom Ellipsoid 
die Werte der potentiellen Dipolfunktion /’, auf der Meridiankurve in erster Annaherung 
den Abstanden von der Achse proportional sind, daf also die Kurve (/’,)c, in Abhangig- 
keit von x aufgetragen, der Kurve 7 (x) ahnlich wird (Abb. 2). Ob der Satz weiterer 
Verallgemeinerung fahig ist, um insbesondere zur Auffindung der Funktion F, bei- 
zutragen, soll hier nicht untersucht werden. 


Vif. Bestimmung einer weiteren Naherung fiir das Potential der Querstrémung. 


Um eine gute Naherung fiir die Funktion F, der Querstrémung um den aus Quell- 
punkt und Senkenstrecke erzeugten ovalen Rotationskérper aufzustellen, wollen wir 
eine quadratische Dipolverteilung zugrunde legen, die sich zusammensetzt aus einem 
linearen Dipolkeil von der oben bereits bestimmten Art und einem rein quadratischen 
Keil, dessen Starke vom Werte Null (bei « = —c) bis zum Maximum (bei x = c) 
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zunimmt. Mit zwei zunichst unbestimmten Konstanten B, und B, haben wir fiir 
diese beiden Verteilungen die Potentialfunktionen 


ry ’ (c+ x) sin Odd H 1 
(F.)r = — B, | ; PL aT B) * - 
(Fa) = — B,\ "+ eRe 


Oh 
Setzen wir wieder wie oben vortibergehend 
c+ta—x+c— netg #, 


so erhalten wir fiir die zweite Funktion 


2 2 
(a = ¢)? cos? 8 dd 


| sin 99 +n | “5 
1 ab 


4 


(Ff); a B, 


2 (x +c) | cos 9.as| 
, : 


oder nach Ausfiihrung der Integration und Einsetzung der Ellipsoidkoordinaten 


Cpe ee ae C 1 I 2, “ 
(F'.)77 = 2 Bec V (2 1d 4) a4 7 Ig ae (Gena) mar (59) 
Im ganzen haben wir also 
B 2 
F, aad (F); ate (Py)r = 2 - (a= TES TE ae Bin Q, (ee Tic A) | (60) 


Da die Dipolverteilung bereits der Form des Rotationskérpers angepaBt ist, so geniigt 
es, die peiden Konstanten B, und B, so zu bestimmen, da in zwei Punkten, das 
heiBt fiir zwei Wertepaare ¢ uw der Meridiankurve C die genaue Tangentialbedingung 
fiir die Geschwindigkeit erfillt ist. Ftir die Geschwindigkeiten erhilt man mit Hilfe 
der Umrechnungsformeln 


a B(C—1) . ac CVC@—)A—-) . 
po 2 yt) ? me pS, /2 > 
Ox Ce) on ¢ (2 ==?) (61) 
om f(1~y*) om wV (2 —1) (1 — ) 
OG CNC? = yr)’ Cine e (2 — p?) 
die folgenden Ausdriicke: 
Oye Pep IBS 7) il ree , $A *}4 
eS Ce (Sn Rr ie fan ae aa 
7] 1 cu—l| 
+43, Gaye Ht ee me 
On 2B, I cos ES) + 2 BLQ 
pee. he (Cc a ee 8 Gr 201 
4B, [wuid— Cp) La 
seer hea | 
Die Tangentenneigung in einem Punkte der Kurve C' wird mit der oben angegebenen 
Formel fir a 
= 1 Mere ores = 
tee =>) at ettaS C Snyest Saye Pe 
Die allgemeine Tangentialbedingung lautet dann 
Un = Ve 
ax a A V» tg é. 


Da es in der Folge vor allem auf die Auswertung des Integrals 


\ Fy dz 
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ankommt, wahlen wir zwei Punkte mit gréBeren Abstiinden von der Achse, also 
etwa den Punkt mit der gré8ten Ordinate, also tg 6 = 0, und den auf der y-Achse 
gelegenen Punkt «= 0. Fir tgé—=0 erhalt man eine quadratische Gleichung 
fir «= mit der Lésung 


ee tae Ose W\F . 
‘ os ( 2e ee 


und damit die Méglichkeit, zusammengehérige Werte von cy za berechnen. Die 
entsprechenden Kurven X (¢,) bringen wir in der ¢ u-Ebene zum Schnitt mit den 


——— = 


o4 
O6 


sees 
——$ || Ce 
AJ 
S' 
S 
NN 
) 


NG) 


10 47 le ¢ - 43 14 


Abb. 3. Zur Bestimmung der Punkte an der Meridiankurve. 


Tabelle 3. Berechnung der ausgezeichneten Punkte der Meridiankurve C. 
a) Kurve X (€ uw) als Ort der Punkte tu (0 = 0). 


1 1°01 1°02 1°04 1°05 1°06 110 | 1°20 1°25 11633 14 | 15 
1 0°820 | 0°757 | 0°673 | 0°643 | 0°618 | 0°544 | 0°429 | 0°391 0°366 | 0°361 | 0°358 
b) Koordinaten € uw der Punkte. 
b | | 
= Ol 02 o4 06 os 10 12 
(ea iia) IW) teres io a 10494 1°:0976 119115 | 1°:2809 1°3676 14515 1°533 
at (i A Nicer ole c 1°00125 | 1°00497 | 1°01961  1°04315 1°0744 11121 1°1551 
Cette eee eeeee ees 1012 1°03 1078 {= Tei les | 1193 1°255 1°32 
Wiech erator eteke 0°806 0°706 0°582 | 0°498 | 0°435 0°388 0°350 
ar (i 0) ieetirnereenen 1°00705 | 1:00985 | 1°0378 | 1:081 ; 171325 | 1°1915 | 1°2546 
| | 


Bildkurven C, (¢) der Meridianbegrenzung, die sich nur durch den Wert der 


Konstanten b unterscheiden (Abb. 3). Damit ergeben sich die Wertepaare ¢j fiir 
die Koordinaten der Aquatorpunkte, Einsetzung in die Gleichung . 


pre ba or (65) 
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liefert dann eine Bestimmungsgleichung fiir die beiden Konstanten. Fiir « = 0 er- 
halt man ferner 


ao 2 B v BO (CP 1 
Vo lee: 2 B,|; ee = 2 B, (Q,):, * a ad ofl 
wo ¢, aus der Gleichung 
62 
bo (C2 — 1) = as 
zu bestimmen ist [vgl. (39)]. Die zweite Bestimmungsgleichung wird daher wegen (63) 
bigeOe case iy 
tg 0) = (tg 6), = a (66) 
nach Vereinfachung 
Ese) seine) 1 il 
a Fa] + B, x 160 (3 fo" Da? er 4 Co (3 fq" — 1). (67) 


Mit Hilfe der Gl. (65) und (67) erhalt man dann die in der Tabelle 4 zusammen- 


gestellten Werte der Konstanten fiir eine Reihe von Werten des Verhiltnisses 8, 


Tabelle 4. Berechnung der Tragheitskoeffizienten k, und ky. 


4. 4x B, “. k= — Ky | hig — hy 

02 0°01168 — 0°001975 34°289 0°063 0°908 0°845 

04 0°04505 — 0°00864 8°5362 07119 | 0°842 0°723 

0°6 0°0957 — 0°01953 3°7942 0°167 | 0°778 0-611 

0°8 0°1493 — 0°02582 2°1444 0°208 0°721 | 0°513 

10 0°2144 — 0°03408 1°3652 0°241 0°678 0°437 
| | 


Wenn man andere Punkte wahlt, ergeben sich nur geringe Abweichungen von den 
berechneten Konstanten, so daf die Naherung als praktisch ausreichend angesehen 
werden kann. Eine weitere Naherung kann man dadurch erhalten, da man durch 
Hinzufiigung etwa einer Dipolstrecke mit konstanter Verteilungsdichte eine weitere 
Konstante aufnimmt und zur Bestimmung der drei Konstanten B, B, B, die 
Tangentialbedingung fiir drei Punkte der Kurve C in Ansatz bringt. 


VIII. Ermittlung des Traigheitskoeffizienten &, der Querbewegung. 


Auf Grund der Entwicklungen im ersten Bericht ergibt sich der Koeffizient k, 


aus der Formel 
Fy dx 
ky = — ' - ep (68) 


wobei die auftretenden Integrale iiber die Meridiankurve zu erstrecken sind wie im 
oben behandelten Falle der Langsbewegung. Nun erhalt man aber fiir den Wert 
der Funktion F, in den Punkten der Kurve C bei der Benutzung der Gl. (39) 


(F.) = 2 BinQi = (=: 2 B,) 7. (69) 


Der Ausdruck setzt sich aus einem mit 7 proportionalen Glied und einem vom Produkt 
7 Q, abhangigen Glied zusammen. Das fiir die Berechnung von k, maigebende Integral 
wird daher 


—|Fyndx == (2 1 2.B,) | ade — 2B, | 4Q,dz, 


Cc 


d 
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womit der Koeffizient k, den Wert 


| 77 Q, dx 
| 72 da 


402 / B, 
hy = b? is 


+ 2 B,| — 2 By (70) 


erhalt. Da das mit dem Volumen des Rotationskérpers verhiltnisgleiche Integral 


| n? dx 
bereits bestimmt ist, kommt die weitere Rechnung auf die Ermittlung des Integrals 
( G 1 Gus aa!! 
One welea ele da (71) 


hinaus, das, wie man zeigen kann, in geschlossener Form ausgewertet werden kann. 
Setzt man fiir 7? da die frither gefundenen Ausdriicke ein, driickt also alle Gro8en 
durch die Koordinate ¢ aus, so kommt 


2 
bae\ (PP e+ tae — spr — 2 IP ee FD) 
1 

diperer siete ile (72) 

Benutzt man die friiher angegebenen und die weiteren Integralformeln 

| = i) ; 1; | ae << : 0-47 a 

\ Se -pe  2 | SS =e 42-0: 

ego = ee (73) 


so erhilt man nach Zusammenfassung und bei Einfiithrung von 


Legge cate t ; dQ f hinge ie 
Qo gy BS es Oe Gee eee 
den folgenden Ausdruck fir das unbestimmte Integral 
il ( Diab? te | ax. e 1 bf (3241 
of a on I Ge et oa ay 4 
' bs 8 1 ] 1 2 1 “4 
DAR Sh(G2ESSn)2 16 c& | (@— 1? Qo 4 4 07 ace (2? — 18 Q, + goll Saer Ty ae 


der schlieBlich noch vereinfacht werden kann, wenn man die fiir beide Grenzen 
giltige Beziehung 


5c(@—-N)-@—1s- 77 =0 (75) 
benutzt. So entsteht 

: jl, GG ea Ca Sete ph ete 

3 L eee C2 Ig KG 1) 4 ct ee | Q, | 12 ¢! (C2 — 1)2 


be 1 ] 1 : z 
= 16 6 (¢# — 1)2 Qo 4 Q. miners ee a. 18 Q, +36(I eee zo. (76) 


Damit gewinnt man fiir den Tragheitskoeffizienten k, den Wert 


ee 


4c? / B, 
b? Cc 


L 
+ 2 B,) 2 (77) 
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IX. Ergebnisse der Rechnung. 


Wir haben die fir die Bestimmung der Tragheitskoeffizienten notwendigen Hilfs- 
groBen in den Tabellen 2, 3 und 4 zusammengestellt, insbesondere die Koordinaten 
der benutzten ausgezeichneten Punkte der Meridiankurve C, die Integrale J, K, L 


und die aus der Tangentialbedingung gewonnenen Konstanten ey B, in dem 


Potential der Querstr6mung. Die Tabelle 4 gibt dann die Werte der Koeffizienten k,, 
k, und k, — k,. Bei der Bewertung der Ergebnisse interessiert vor allem der EinfluB 
auf die Tragheitskrifte bei beschleunigter Bewegung und auf das instabile Moment 
um die Hochachse bei einer Kursabweichung von der Lingsachse. Bei Auftreten 
einer Beschleunigung in der Langs- bzw. Querrichtung werden die Tragheits- 
widerstande 


We hey Oa: W250. V (78) 


wo k, @ V bzw. k, @ V auch als die den beiden Richtungen zugeordneten scheinbaren 
Massenkomponenten bezeichnet werden. Bei einer Horizontalbewegung unter dem 


Winkel 6 gegen die als horizontal vorausgesetzte Langsachse entsteht ferner das 
Drehmoment um die Hochachse 


] : 
M ==; 9 U?sin 2B (ky — hy) V, (79) 
wenn 
0; = U_COs 6, 1-05. Using 


die Bewegungsgeschwindigkeiten in der Langs- und Querrichtung bezeichnen. Fiir 
den dimensionslosen Beiwert hat man 


Cm = (ky — ky) sin 28 = =——. (80) 
op Q U2 V 


Aus den Ausdriicken von W, und W, und M ergibt sich, daB es am zweckmAafigsten 
ist, als Vergleichskorper ein Ellipsoid zu wahlen, das mit dem ovalen Rotationskérper 
das Volumen und etwa die Lange / gemeinsam hat. Dann ergibt die Differenz k, — hk, 
direkt eine Handhabe fiir den 
Vergleich der bei gleichbleiben- 74 
dem Winkel 6 entstehenden Mo- 
mente beider Korper. Die Lange 
der benutzten Rotationskorper, 
bezogen auf die Grundstrecke 
c, ist 


14 


egy 


10 


l 
go oi + Se 


Ge 


Benutzen wir das bereits be- 
stimmte Volumen V, so ergibt ge 
sich aus 


4 l a4 
= Hig ; 
V = 3 0b 3 


G2 


als kleine Achse 0’ des Vergleichs- 


ellipsoids é 
UWP Er Clams OF) Ob UL > C6 HOS, 


10 


ey ES oa a a ch 
C i Cc Qntl ¢ Dah (81) Abb. 4. Abmessungen u. Koeffizienten der Rotationskérper. 
19* 
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Tabelle 5. Werte fiir das Ellipsoid gleicher Lange und gleichen Volumens. 
a SS oe 


. ion a ie (ide (dg (is — kg 
0-2 0°02685 0°16385 0°1560 0-038 0°920 0°882 
0-4 0°09703 0-315 0°2818 0-097 0°849 0°752 
0°6 0°19901 0°4461 0°3839 0°148 0°783 0°635 
0°8 0°3253 0°5738 0°4701 0°193 0°726 0°533 
1-0 0°4727 0°6876 0°5370 0-228 0°683 0°454 


Aus dem damit bestimmten Achsen- 
verhaltnis des Ellipsoids erhalt man dann 
aus den friiher angegebenen Formeln die 
Werte fiir die Tragheitskoeffizienten (k,), 
und (k,),. Das Ergebnis ist in der Tabelle 5 
wiedergegeben. Ferner haben wir in 
Abb. 4 die Abmessungen und Koeffizien- 
ten in Abhangigkeit von b/c und der Abb. 5 
die Werte der Koeffizienten fiir das Ovaloid 
und das zugeordnete Ellipsoid in Abhan- 


: eee 
gigkeit vom Verhaltnis 2 zur Anschau- 


ung gebracht. Man sieht, daB die Koeffi- 
zienten k, bzw. die Tragheitswiderstande 
bei der Liingsbewegung groBer und die 
Koeffizienten k, bzw. die zugehorigen Trag- 
heitskrafte bei der Querbewegung etwas 
kleiner ausfallen als die entsprechenden 
GroéBen des Vergleichsellipsoids. Der Unter- 
schied ist am starksten bei den schlanken 
Formen, bei denen gleichzeitig das Breiten- 
Abb. 5. Vergleich der Koeffizienten k, und k, maximum nach vorne verlagert ish, Wal 
mit den Werten beim Ellipsoid gleicher Lange Tend mit wachsendem 6 die Differenz in 
und gleichen Volumens. dem Mae verschwindet, wie die Rota- 
tionsform sich dem Ellipsoid annahert, was 
gleichzeitig als Bestaétigung der Rechnung anzusehen ist. Wichtig ist fiir uns vor 
allem, da die dem Moment verhiltnisgleiche GréBe k, — k, gegeniiber dem ent- 
sprechenden Wert beim Ellipsoid verkleinert wird, besonders bei ovalartigen 
Rotationsformen, die im grofen und ganzen, das hei®t im Schlankheitsverhaltnis 
und in der ungefaéhren Lage des gréBten Querschnittes den gangbaren Rumpfformen 
des Flugzeuges entsprechen. Theoretisch und ohne Riicksicht auf die Reibung wiirde 
sich also das merkwiirdige Resultat ergeben, daf diese Rumpfformen eine geringere 
Instabilitat besitzen oder leichter zu steuern sind als ellipsoidische Formen gleicher - 
Lange und gleichen Volumens. Das ist auch insofern ein bemerkenswertes Ergebnis, 
als gerade Rotationskérper dieser Art, also mit ,,stromlinienartigem‘‘ Umri8, auch 
durch ein relatives Minimum des Formwiderstandes ausgezeichnet sind. Wie der 
Verfasser in einem weiteren Bericht zeigen wird, bleibt das Ergebnis im wesentlichen 
auch fiir Formen giiltig, die sich von den ovalartigen Rotationskérpern durch eine 
starkere Ausbildung des Kopfes unterscheiden. Wenn sich also in der Wirklichkeit 
oder beim Versuch andere Wirkungen zeigen, so kann nur die Reibung dafiir ver- 
antwortlich gemacht werden. 
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X. Beziehung der Koeffizienten zu den Geschwindigkeiten. 


Wie ich bereits im ersten Bericht bemerkt habe, lassen die Koeffizienten k, und k, 
im Falle des Ellipsoids eine einfache Deutung zu, die zu einer interessanten Ab- 
schatzung im allgemeinen Falle fiihrt. Es li8t sich namlich zeigen, da beide 
Koeffizienten mit den Faktoren identisch sind, mit denen man die Geschwindigkeiten v, 
und v, der Lings- bzw. Querbewegung multiplizieren mu8, um die absoluten Ge- 
schwindigkeiten am Aquator bzw. an den Polen oder Scheiteln des Ellipsoids zu 
erhalten. Abgesehen vom Vorzeichen, hat man daher bei der Langsbewegung an 
der breitesten Stelle die Relativgeschwindigkeit 


0. sao) (1 Se ky), 


Bei der Querstrémung entstehen ferner in den Scheitelpunkten des Ellipsoids in der 
Symmetrieebene p = 0 die relativen Maximalgeschwindigkeiten 


[v,’| = v2 (1 + ky). 
Um nun z. B. das Verhaltnis 


ly 
ev es OF 
4 | 

Vv} 


Tmax 


fir den Fall der Langsbewegung des ovalen Kérpers zu berechnen, gehen wir aus 
von dem allgemeinen Ausdruck 
ee dees ald ia ee Mu 
| (¢ — wu)? C+ yu 
und benutzen die unter (64) gegebene Beziehung fiir die Koordinaten ¢ « des Punktes 
der Meridiankurve mit groé8tem Abstand von der Achse. Man erhalt dann nach 
Vereinfachung 


P (82) 


CNM al Oe ors uareeseay 7, 


ease Uae ike (1 

etic eS aN RON ine 
Ebenso erhalten wir aus der grof%ten Geschwindigkeit langs des Aquatorkreises bei 
der Querbewegung 


is == 


(83) 


do ie 
Ui 


Ve 
Ve 


Pied 53 
— 2B, (Qi +a (=P + 2 Ba. (84) 


Tmax 
Weiter entstehen aus (62) fiir die auf die Scheitelpunkte bezogenen Geschwindigkeits- 
verhaltnisse bei Benutzung von (39) die einfachen Ausdriicke 


ce / B 
fest! => — iB, (Q1)c, [. 4 b2 ( - | Y B,), 
renee (35) 
aaa ¢ ¢ 1 ¢ 
hy!” = —2 Ba (Qr)e, + 4G (+ 2 Ba). 
Tabelle 6. Vergleich mit den Geschwindigkeiten am Aquator und an den Polen 
(Scheitelpunkten). 
a ee 
b Sie y In nem icahs A k.’ pire |e. 
© : 1% amxa _ | x i: Yaly = 1 # F i : V2 Iu =—1 
0°72 0°090 0°063 0°038 0°789 0°908 0°832 1°160 
0°4 0°152 07119 | 0°097 0°722 0°842 -0°776 1:099 
6°6 0°197 0°167 | «+0°148 0°666 | 0O°778 0°719 1°015 
0°8 0°234 0208 07193 | 0°644 | 2O°721 | 0693 | 0°883 
1:0 0°266 0°241 | 0°228 0°617 | 0°678 | 0°659 | 0°805 
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Aus der Tabelle 6, in der die sieben Werte k,, ky’, (hy)z, ba, hy’, hea", he” zusammen- 
gestellt sind, erkennt man, da® z. B. die Triagheitskoeffizienten k, zwischen den 
Werten k, und (k,), liegen oder daB die Ungleichung 

(ki)p < ky < hy (86) 
besteht, so daB man berechtigt ist, bei Kenntnis dieser Werte ihr arithmetisches 
Mittel als Anniherung fiir k, zu benutzen. Was die Querstrémung angeht, 80 stellen 
die beiden Scheitelpunkte k,’’ und k,'”” zwei Extremwerte dar, die zugleich die obere 


und untere Grenze jsowohl fiir die Tragheitskoeffizienten /, (und k,') wie auch fiir 
den entsprechenden Ellipsenwert (4), bilden. Da8 z. B. 


heh ks OR (87) 
wird, sieht man iibrigens ohne weiteres ein, wenn man die Ausdriicke (84) mit (70) 
vergleicht und beriicksichtigt, daB (Q)z, der kleinste Wert von Q, und (Q,);, der 
groBte Wert von Q, auf der Meridiankurve darstellt. Aus der Ungleichung 


\ 1? Q,d 


(Qi), — [dx —< (Q1)c, 


: pe Oo arare 
ergibt sich dann sofort die obige Ungleichung fiir die Koeffizienten. Wenn -- wachst, 


so nahert sich der Rotationskérper mehr und mehr der ellipsoidischen Form und 
schlieBlich der Kugel mit unendlich groBem Halbmesser. Dabei riicken die oberen 
und unteren Grenzen immer weiter zusammen, um schlieBlich in den Grenzwert 


be he SS ae = OS 


tiberzugehen. 
XI. Schlu8Bbemerkung. 


Mit Hilfe der angegebenen Methode diirfte es nicht schwer sein, weitere Rotations- 
formen durchzurechnen. Um z. B. den Einflu8 einer besonderen Buggestaltung auf 
das Moment zu bestimmen, kann man den bisher benutzten Singularitaten weitere 
einfache Quell- und Dipolelemente hinzufiigen. Bei Anbringung eines einfachen 
punktformigen Dipols in dem Endpunkt C, der Grundstrecke erhalt z. B. der Kopf 
des Rotationskorpers eine kugelformige Ausweitung, wahrend bei entgegengesetzter 
Richtung des Dipols der vordere Teil der Begrenzung verengt bzw. eingezogen wird. 
Wenn es auch im allgemeinen nicht moéglich sein wird, bei solchen zusammengesetzten 
Fallen die auftretenden Integrale wie im vorliegenden Fall in geschlossener Form 
darzustellen, so bietet doch die allgemein begriindete Zuriickfithrung der Momenten- 
berechnung auf einfache Flachenauswertung eine wesentliche Erleichterung der 
gesamten Zahlenrechnung. Ich hoffe, in einem weiteren Bericht einige praktisch 
bemerkenswerte Ergebnisse mitzuteilen, insbesondere auch die Frage nach dem 
Einflu8 einer Ausweitung des Rumpfkopfes auf die Stabilitat zu behandela. 


(Eingegangen am 3. Juni 1954.) 


Kegelige Uberschallstrémung in Schallnihe. 
Von K. Oswatitsch und L. Sjédin, Stockholm. 
Mit 4 Textabbildungen. 


Zusammenfassung. Die tblichen Tabellen oder Kurvenblatter fiir die Uberschallstromung 
an axial angestr6mten Kreiskegeln reichen in Schallnihe nicht aus, weil sie nicht geniigend dicht 
gerechnet sind. Die Herleitung entsprechender Tabellen und Kurven unter den in Schallnahe 
zulassigen und iiblichen Vernachlassigungen fiir schlanke Kegel in dieser Arbeit gestattet die 
Riuckfihrung der Resultate auf-einen einzigen aus Kegeloffnungswinkel und einer Machzahlfunktion 
gebildeten Parameter fiir alle Verhaltnisse der spezifischen Wirmen. Die Darstellung enthalt 
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die Ahnlichkeitsgesetze fiir Schallnaéhe und gestattet gewisse’ formelmaBige Aussagen, wie etwa 
fiir die Machzahl der Kopfwellenablésung. Mit dem Aquivalenzsatz fiir Schallnihe kénnen die 


Sileetis auf Kegel beliebiger Querschnittsform und auf die Spitzen von Deltafliigeln tibertragen 
werden. 


Summary. The usual numerical tables and charts for axial supersonic flow around circular 
cones are inadequate in the sonic range because there is a lack of calculated data in this region. 
By deriving tables and curves under such assumptions as are permissible and usual for slender 
cones in the sonic range, the results can be reduced to a single parameter which is derived from 
the cone aperture angle and a function of the Mach number for all specific heat ratios. The 
similitude laws for the sonic range are included, and some formulae are given, as, for instance, 
for the Mach number at the bow wave separation. By means of the equivalence theorem for 


the sonic range, the results may be transferred to cones of arbitrary cross-section, and also to 
delta wing tips. 


Résumé. Les tableaux et graphiques usuels pour un courant plus rapide que le son le long des 
cones circulaires exposés & un courant axial ne suffisent pas au voisinage de la vitesse du son, 
car ils sont calculés avec des intervalles trop larges. On trouve dans ce travail la déduction des 
tableaux et graphiques appropriés, les termes négligeables usuels pour les cénes aigus au voisinage 
de la vitesse du son étant omis. Cette déduction permet de ramener les résultats & un paramétre 
unique calculé & l’angle au sommet du céne et une fonction du nombre de Mach, pour tous les 
rapports des chaleurs spécifiques. L’étude comprend les théorémes de similitude au voisinage de la 
vitesse du son, et permet d’énoncer certaines formules comme par exemple celle relative au nombre 
de Mach de séparation des ondes de choque. Avec le théoréme d’équivalence au voisinage de 


la vitesse du son, les résultats peuvent étre appliqués A des cénes de toute section et a la pointe 
@une aile delta. 


I. Einleitung. 


Die Uberschallstrémung an Kreiskegeln ist vielfach berechnet worden. In diesem 
Zusammenhang braucht nur auf die in allen einschligigen Lehrbiichern zitierten 
Arbeiten von G. J. Taylor-J. W. Macoll, A. Busemann und auf die Tabellen 
von Z. Kopal hingewiesen werden. Dabei handelt es sich stets um die numerische 
Integration einer gewohnlichen, nicht linearen Differentialgleichung, die keinerlei 
grundsatzliche Schwierigkeiten bietet, wenn auch zwischen Kopfwelle und Kegel- 
mantel in Schallnahe Uberginge von Uberschall- zu Unterschallstrémungen auf- 
treten. Wenn das Problem der Uberschallstromung von Kreiskegeln fiir beliebig 
diinne Korper und fiir Machzahlen der Anstromung MM. beliebig nahe an Eins 
(M., —~ 1) demnach grundsitzlich gelést erscheint, solange nur die Kopfwelle anliegt, 
so besteht dennoch ein besonderer Bedarf fiir die Lésungen bei schlanken Kegeln 
in Schallnahe. In diesem Gebiete reichen die vorhandenen Tabellen und Diagramme 
in keiner Weise aus. Aus diesen lassen sich keine allgemein giiltigen Formeln und 
Ahnlichkeitsgesetze fiir Schallnahe ableiten und es fehlt der AnschluB an die tiblichen 
Methoden fiir schlanke Korper in Schallnahe. 


Was diesen ersten und langeren Teil der Arbeit angeht, handelt es sich also um 
die methodische Aufgabe der Darstellung der kegeligen Strémung an Kreiskegeln 
in Schallnahe. Dabei unterscheiden sich die Resultate nur in ihrer Darstellungsweise 
von einer schon friiher ausgefiihrten, aber nur in Berichtsform vorhandenen Arbeit 
von v. Baranoff}. 


Ausgehend, von der Kreiskegelstro6mung ist es dann nicht mehr schwer, mit Hilfe 
des Aquivalenzsatzes von Oswatitsch® Aussagen tiber beliebige kegelige Stromungen 
in Schallnihe zu\machen. Uber den Inhalt vorliegender Arbeit wurde auf der GaMM- 
Tagung in Miinchen (1954) referiert. 


1 vy, Baranoff: Similtude transsonique pour les corps de révolution minces sans incidence. 
ONERA NT 1/1329 A 1952. g ; 

2 K. Oswatitsch und F. Keune: Hin Aquivalenzsatz fiir nicht angestellte Fliigel kleiner 
Spannweite in schallnaher Strémung (in Druck). 
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Il. Ausgangsgleichungen fiir die Kreiskegelstromung. 


Die Differentialgleichungen und die StoBgleichungen fiir Schallnahe haben unter- 
schiedliche Darstellungen gefunden, die sich jedoch nur unwesentlich voneinander 
unterscheiden. Der Einfachheit halber sei hier von der Darstellung im Buche ,,Gas- 
dynamik“ des ersten Verfassers* ausgegangen. Ist x die Kegelachse und y die Quer- 
koordinate und sind wu, v die entsprechenden Geschwindigkeitskomponenten, so seien 
mit u., und «,, als Geschwindigkeit und als Machwinkel in der Anstr6mung und 
mit c* als kritischer Schallgeschwindigkeit folgende schallnahe Groen eingefihrt: 


C8 00 = VMoF—1; Y= Y Cte X03 | 
U — Uco Vv 

Woo —c*? ack prnmcer Bl Soo: 

In diesen Groen schreibt sich die fiir Schallnahe vereinfachte gasdynamische Gleichung 

fiir Achsensymmetrie und die Gleichung der Wirbelfreiheit : 


(1) 


Mie 


@ 1 a(v by) 
(l + u) = = = 0; (2) 
ou ov 
Machzahl M = 1 entspricht dabei dem Werte wu = c*, das hei®Bt nach Gl. (1) u = — 1 


und das bedeutet das Verschwinden des ersten Koeffizienten in Gl. (2) oder einem 
Wechsel vom Hyperbolischen zum Elliptischen im System (2), (3). 

Bei einer Kreiskegelstr6mung hingen die Geschwindigkeitskomponenten nur von 
einer einzigen Variablen ab, als welche gewahlt sei: 


Wy) y 
1 8 as OLE ase (4) 
Nun gilt: 
ou du 4 ou its Lien 
6a -> Tidy ae an ~— ody’ 
ov dy ow dy 1 (5) 
Cr an a’ Cy nem 
Mit Gl. (5) geht das System (2), (3) ttber in 
du d(v 7) 
TL) gee ge ee (6) 
dit d(v 7) v” 
dn " eee bea 0. (7) 


An diesem Gleichungssystem sind die bekannten Eigenschaften von » und u an der 
Achse 7 = 0 leicht zu erkennen. Wenn u auf der Achse nicht stark singular ist, so 
ist das Produkt »-7 in der Nahe von 7 = 0 nach Gl. (6) offenbar konstant: 


(Mama es a) (8) 
Angenommen, auch dy/dy verschwindet fiir 7 —> 0, so folgt aus (7) fiir 
yO: we an(0) Inga a (9) 


Kin Kinsetzen dieser logarithmischen Singularitaét in (6) zeigt, daB das Verhalten der 
Komponenten an der Achse durch Gl. (8) und (9) wirklich richtig getroffen ist. 
Nun liegt die Oberflache schlanker Kegel in Achsennahe. Die Werte auf der 
Achse sind daher sehr erwiinscht, jedoch nicht von solchen Gré8en, welche wie u 
und » sich an der Achse stark andern, ja iiber alle Grenzen wachsen. Daher seien 
die Komponenten durch zwei neue Variable ersetzt, welche in Achsennihe konstant 


5’ K. Oswatitsch: Gasdynamik, 8. 342ff. Wien: Springer-Verlag. 1952. 
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sind, in welchen also das singulire Verhalten dort eliminiert ist. Als solche bieten 
sich mit Gl. (8) und (9) an: 

y(n) =v (n)> 9, 

%(n) = u(y) — y (q) ny. 

In x (y) steht vor dem Logarithmus nicht y (0) wie in Gl. (9), sondern das lokale 

y (y). Wenn man davon absieht, ist jedoch x in K6rpernahe im wesentlichen gleich 
jener GréBe, welche in anderen Arbeiten vielfach als ,,RaumeinfluB bezeichnet 
wird‘. In diesen Veranderlichen erhalten die Gl. (6) und (7) wohl eine etwas kompli- 
ziertere Form, welche jedoch nun frei ist von jeder Singularitat im Integrationsgebiet. 
Nach elementarer Umformung erhilt man: 


(10) 


dy miter coy ben) o 
dn ¥ ePil+yt+tyiny]—1’ (lt) 
d d 
a= ie (1 +Ing]. (12) 
Als Randbedingung auf der Achse sei dort der Wert von » vorgegeben: 
7=0: p=y (0). (13) 


Er kann mit Gl. (10) und (1) auf den Kegeloffnungswinkel zuriickgefiihrt werden, 
was spater geschehen soll. Die tibrigen Randbedingungen sind durch die Sto8- 
gleichungen gegeben. Die Werte unmittelbar hinter dem Sto8 werden durch ein 
Dach gekennzeichnet; dann gilt? mit y als StoBfrontwinkel: 


b= 1 +5 fl, (14) 
A it ere =e 
i = ctg oo tey = = = [tt yi] (15) 


Diese Bedingungen, in die neuen Veranderlichen tibersetzt, lauten nach einiger 
Umrechnung: 


A bake ana A A 
n=H=|1- 33] 4 = (be lH): (16) 


Wie immer ist die Bedingung am Sto dadurch gekennzeichnet, daB die StoBlage 7 
von der Sto8starke und damit von » abhangt. 


III. Die Integration der Differentialgleichungen. 


Wie schon im allgemeineren Fall, ist die Integration des Differentialgleichungs- 
systems (11), (12) allgemein nicht méglich. Sie muB numerisch oder graphisch durch- 
gefiihrt werden, wobei nicht von einer gegebenen Randbedingung auf der Achse (13) 
ausgegangen werden kann, da die GroBe dort ja gesucht wird. Es wird daher, wie 
schon friiher, der Wert hinter dem StoB 7 = 7 von y vorgegeben. Mit Gl. (16) folgt 
aus einem vorgegebenen Wert 7, die StoBlage 7 und der Wert y. Mit diesen Anfangswerten 
kann die Integration des Systems (11), (12) schrittweise durchgefiihrt werden, woraus 
sich die Werte von y und x auf der Achse fiir 7 = 0 ergeben. In Bereichen gro8erer 
Verinderlichkeit von 7 und y wurde Runge-Kuttas-Methode® angewendet, wahrend 
im ibrigen mit einer einfacheren Mittelwertmethode gearbeitet wurde. Die Rechnungen 
wurden fiir die Anfangswerte » = 0°10, 0°20, 0°40, 0°60, 0°80, 1°00, 1°20, 1°40 durch- 
gefiihrt. Héhere Werte von y haben keine praktische Bedeutung, da sich die maximalen 


4 fF. Keune und K. Owatitsch: Nicht augestellte Kérper kleiner Spannweite in Unter- 


d Uberschallstromung. Z. f. Flugwissenschaften 1/6, 137—145 (1953). 
a 2 Gs ee Scarborough: Numerical Mathematical Analysis. The John-Hopkins Press. 1930. 
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Kegelwinkel im Bereiche 1:00 < p < 1:20 ergeben. Die Ergebnisse sind in Tabellen- 
form am Schlu8 der Arbeit und in Diagramm 1 wiedergegeben. In letzteres sind 
keine 7-Werte aufgenommen. 

Ausgehend von der Sto8polaren 7 = 7, y = y steigen die Werte von y mit 
fallendem 7, wenn Uberschall 


M>1: 0<1+u=l1+y+ylnn (17) 


herrscht und fallen mit fallendem 7 im Unterschallgebiet. Im y, y-Diagramm 1 ist 
die Schallstelle durch eine deutliche Spitze in den Kurven fiir einen bestimmten 
StoBwert » gekennzeichnet. Ausgehend vom Wert y auf der StoBpolare steigen die 
Werte von y (bei fallendem y) an bis zum Erreichen der kritischen Geschwindigkeit 
und fallen dann zum Werte y (0) auf der Achse 7 = 0 ab. Die Kurve dieser End- 
punkte auf 7 = 0 wird als ,,Apfelkurve‘‘ bezeichnet, obwohl die Apfelform in Dia- 
gramm 1 nicht mehr zu erkennen ist. Jedoch werden aus dieser Kurve die Werte 
am Kegel entnommen. Die Apfelkurve schneidet in dieser Darstellung die StoB- 
polare bei etwa y = 0°56. 

Die lineare Losung entspricht in der vorliegenden Darstellung kleinen Werten 
von u und » GI. (1). Das bedeutet kleine Werte von y und y oder ein Streichen aller 
nicht linearen Glieder, das hei8t der y- und y-Terme im Bruch von Gl. (11). Die 
Gleichungen lassen sich dann integrieren. Ebenso gut kann man die bekannten 
Resultate der linearen Theorie fiir die Kreiskegelstroémung auf die neuen Variablen 
umschreiben. Man erhalt: 


y= (0)V1— 9; 


x= (0) (1 —V1 — 7?) Ing — (QL + V1 — 9), 

oder die Apfelkurve ftir die linearen Gleichungen: 
4 (0) = — p(0)In2. (19) 
Auch diese Gerade ist in Diagramm 1 eingetragen. Wie erforderlich, ist sie die 

Tangente an die schallnahe Apfelkurve im Punkte y= y = 0. 

Damit 1a8t sich auch gleich der Unterschied von schallnaher und linearer Lésung 
ablesen. Der Einflu8 der nicht linearen Glieder auf y ist wesentlich gréBer als aut 
die w-Komponente, weil das logarithmische Glied in Gl. (10) den Hauptanteil liefert. 


Man sieht in Diagramm 1 den Fehler der nicht linearen Theorie ,,unter der Lupe“, 
tatsaichlich ist letztere nicht so schlecht. 


(18) 


IV. Auswertung der Resultate fiir Kreiskegel. 


Angenommen der halbe Offnungswinkel des Kreiskegels 9 ist gegeben, so gewinnt 
man den entsprechenden 7-Wert mit Gl. (4): 
# =tg), yn=te ectga,,—ted-/M.2—1. (20) 
Mit einer in Schallnihe hohen Genauigkeit, die jene der Vereinfachung der gas- 
dynamischen Gleichung tbertreffen diirfte, kann nun y am Kegel dem Werte auf 
der Achse gleichgesetzt werden. Mit Gl. (10) und (1) gilt dann: 
Uv U Mo* U — Uco 


y (0) = » (tg ¥ Cte 6.5) =) UE co Ge — Ue Oe og f ET) 


Uco 


y (0) ergibt sich also aus dem halben Offnungswinkel # und der Machzahl der An- 
str6mung, die hier nicht selbst, sondern iiber 


M..* == 


co c* 
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eingeht. Ferner steht in der Klammer neben 1 als kleine Korrektur die Storung 
der u-Komponente. Dieses Glied kann iterativ beriicksichtigt werden, was im all- 
gemeinen sinnvoll sein diirfte. Jedenfalls steckt in der Vernachlassigung dieses Gliedes 
der gréBte Fehler, so da sich aus diesem Term ein oberes MaB fiir den Fehler ergibt. 
Fir schlanke Kegel kann die eckige Klammer durch 1 ersetzt werden, womit dann 
die Loésung im Sinne der Ahnlichkeitsgesetze fiir Schallnahe nur von einem einzigen 
aus tg? und M..* gebildeten Parameter abhangt. 

Mit Gl. (1) und (10) kann nun die w-Komponente angegeben werden, indem 7 
am Kegel wieder dem Werte x (0) auf der Achse gleichgesetzt und dieser der Apfel- 
kurve in Diagramm 1 entnommen wird: 


Uu — Uoco Uoo — c* 
igs wet ane | eh Deu ae 
mae age raperia el 0)-hy (Olin (te Oi cote cee). 


Ist die Lésung fiir das lineare Problem bekannt, so kann auch der Unterschied 
zu dieser aus Diagramm 1 abhingig von yp (0) entnommen werden (Index L: linear): 
U=3 6 ‘eo 
= ee [x (0) — x, (01: (23) 
Zur Erleichterung der Arbeit sei eine Tabelle fiir die 4 
in dieser Arbeit auftretenden Machzahlfunktionen fiir das + 
Verhaltnis der spezifischen Warmen x = 1-400 aufgestellt. 
Die Resultate dieser Arbeit gelten jedoch fiir alle Ver- 7% 
haltnisse der spezifischen Warmen x, wenn man von dieser 
Tabelle absieht. Die tibrigen Diagramme und Tabellen 74 
sind ohne irgendeine Voraussetzung tiber die GroBe von x 
berechnet. Dieses erscheint in den Formeln nur implizit, 
indem — wie in Gl. (22) — neben dem Werte von M,.* 
auch ctg a, auftritt. 
Von groBerem Interesse ist auch der StoBfrontwinkel y. 77 
Er 148t sich mit Gl. (15) aus 7 berechnen: 


WA 


’ Abb. 1. Reduzierter StoBfront- 
te y =n - te-og. (24) winkel. 


Machzahlfunktionen ftir x = 1°400. 


Beene = 1:00 1°02 1°04 1°06 1:08 110 1°12 1°14 116 118 1°20 


Nksseceiod Ww 1°017] 1°033| 1°049) 1°065) 1°O81) 1°097]) T1113) 1°128} 1143) 1°158 
COUKooe | 0 0°201] 0°286] 0°352} 0°408] 0°458]) 0°504] 0°547] 0°588] 0°626] 0°663 
1—1/Moo*| 0 0°016| 0°032] 0°047} 0°061) 0°075] 0°088]} O°101] O°114] 0°125] 0°137 
Oie5 7 oe 90° | 78°6 74°1 70°6 67°8 65°3 63°2 61°3 59°6 58°0 56°5 


Zu diesem Zweck ist in Abb. 1 7 als Funktion von y (0) aufgetragen. Man sieht, 
daB 7 im Bereiche der gré8&ten vorkommenden y (0)-Werte ganz erheblich vom Werte 
fiir den Machwinkel der Anstr6mung, 7 = 1, abweicht. 

In der Naherung des die schallnahe Ahnlichkeit erfiillenden Hauptgliedes von 
Gl. (21) ist der maximal zulassige Kegelwinkel ?,,,, durch den Maximalwert von 


wy (0) gegeben: Vestal sae é 
p ( ) geg tg acre Moo* y (@O)s.= == 0°86 1 adhe pea (25) 
Mit der Niherung im Bereiche von M,,—1 gilt [GI. (IX, 2) in*]: 
1 il 1 ‘ , ¥ 
ee es yO oF = 
1 


Pe 90" ea) aie 


— we--til 


a D2 D4 O6 0b WO) 
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Dies in Gl. (25) eingefiihrt, ergibt, indem fiir kleine Winkel die Funktion tg durch 
den Winkel selbst ersetzt wird, 
Benge? =e (902 Oi Vict wrnreeemet 02 (OO oa) lala (26) 
Sheed alae 
Der hochgestellte Index ° besagt, daB die Winkel in Grad einzusetzen sind. Der 
mittlere Teil der Formel gilt fiir beliebige Gase, der rechte nur fiir Luft. Fiir diese 
bewahrt sich die Naherungsformel (26) auch noch sehr gut bis zu M,, = 1°20, also 
bis zu einer Machzahl, bei welcher die tiblichen Uberschalltabellen mi -diagramme 
bereits ausreichende Dienste leisten. Auf die Tatsache, daB der maximale Kegel- 
winkel in Schallnihe proportional zu ctg «,, ist, hat bereits A. Busemann® hin- 
gewiesen, jedoch fehlte bisher eine Angabe iiber die GroBe des Faktors. 


Eine iiber die schallnahen Ahnlichkeitsgesetze hinausgehende Niaherung wiirde 
auf folgendem Wege zum maximalen Kegelwinkel fiihren. Mit Gl. (21) und (22) ware 
fiir einen beliebigen Winkel #?: 

oo* 0 
se! TW. 1 ee = Bae — Wo vid) : 
bagels a (Oe (tg @ etgaco)] 


Moo U—Uco 


ey 


Nun ist die Funktion In esi sehr empfindlich auf kleine Fehler ihres Argumentes, 
weshalb an dieser Stelle die erste Naherung verwendet werden kann: 


Moot _ tg = y (0) 2) 
1 


“Moo*® Moo*¥—1 1 Moo*¥ —1 
Me 1) + vO) Fn (HES -ctg? aco (0))) 


Aus dieser Gleichung ist nun der Maximalwert von tg # zu bestimmen, der sich 
jedoch vom Resultat der Gl. (25) nur wenig unterscheiden wird, weshalb der Auf- 
wand kaum lohnt. 


VY. Strémung an beliebigen Kegeln. 


Fir schlanke Kegel allgemeiner Querschnittsform, deren Spannweite iiberall 
erheblich kleiner ist als der 6rtliche Durchmesser des von der Spitze ausgehenden 
Machkegels der Anstrémung (Abb. 2), kann der Aquivalenzsatz von Oswatitsch 
angewendet werden’. Nur von solchen Kegeln soll in dieser Arbeit die Rede sein. 
Da der Offnungswinkel des Machkegels in Schallniihe ziemlich groB ist, sind ziemlich 
geringe Pfeilungen flacher Kegel zulassig. ~ 


Nach dem oben genannten Satz ist der Unterschied von gesuchter und linearer 
Lésung derselbe, wie von einem Kreiskegel gleichen értlichen Querschnittes. Das 
bedeutet, da auch in diesem allgemeinen Fall Formel (23) verwendet werden darf, 
wenn unter vw und w, nicht die Werte auf der Oberflache des Kreiskegels, sondern 
irgendwelche Werte in K6érpernahe, beispielsweise auf der Oberfliche des allgemeinen 
Kegels von Abb. 2 verstanden werden. Dabei gehéren die Werte von y (0) und 7, (0) 
zum Kreiskegel gleichen 6rtlichen Querschnittes Q (2), also: 


O(a) = a itg? 0: (28) 
Nach dem Aquivalenzsatz ist nur das Hauptglied in Gl. (21) zu nehmen. Mit 


Gl. (28) ergibt sich demnach der Wert von y (0), fiir welchen das zugehérige x (0) 
in der Apfelkurve aufzusuchen ist zu: 


OF M oo* 
v0) = $2 ey. (29) 


* A. Busemann: Die achsensymmetrische kegelige Uberschallstromung. Lufo XIX/4, 
137—144 (1942). 
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| Die Verschiebung der Geschwindigkeitsverteilung durch die schallnahen Effekte 
ist in Abb. 3 fiir eine angenommene Geschwindigkeitsverteilung wu, an der Ober- 
flache eines allgemeinen Kegels skizziert. Die Ablosung der Kopfwelle bei Anniherung 


Koptwelle 


alg. Kegel Krerskege/ 


Machkege/ 


QOH ] IY 
— u 
——— uy 
Abb. 2. , Abb. 3. 
Allgemeine Kegelstrémung. Verschiebung der w-Verteilung in Schallna&he. 


von M.,-—>1 wird infolge von Grenzschichteffekten, Endlichkeit der Kegel und 
ahnlichem nicht genau an der Stelle maximal zuldssigen Kegelquerschnittes erfolgen. 
Doch ist der maximal zulaissige Kegelquerschnitt des unendlich lang gedachten Kegels 
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Abb. 4. Apfelkurve fiir Schallnahe. 
sicher ein praktisch ausreichendes Kriterium fiir das Ablosen der Kopfwelle. ‘Damit 
ist mit Gl. (29) [in Analogie zu Gl. (25)] die ,,Ablosungs-Machzahl“ gegeben durch: 


1 Qn 4/2 % '. Qaex 
TS Toa OY (a: 118837 on” (30) 


ies 


Aus der Machzahltabelle kann das dazugehérige M,, entnommen werden. 


Fiir die Frage der Ablésungsmachzahl eines Deltafliigels beliebiger Dickenverteilung 
oder fiir die Frage der u-Verteilung an der Spitze eines solchen Deltafliigels erweist 
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sich die Darstellung der Formeln mit Q,,/2 a tiberlegen einer Darstellung mit Q/x?. 
Beide Werte sind ja fiir einen Kegel, Gl. (28), dieselben. Fur die Spitze eines Delta- 
fliigels braucht in (29) und (30) nur (See ), 
wert lim bilden zu miissen, wie es bei Q/2? erforderlich wire. Damit ist auch die 


, eingesetzt werden, ohne einen Grenz- 


Teaee eeh der Ablésung von Kopfwellen bei Deltafliigeln praktisch beantwortet. 


Tabelle 1. 


Schallnahe Kreiskegelstrémung. 


a 


y = 120 


y = 1:00 


? = 0°80 


1°8258) 2°2428 
1°800 | 2°2247 
1°700 | 2°1558 
1°600 | 2°0889 
1°581 


Ah) 


1°4000 | 1°6020 
1°3887 | 1°6020 
1°3446 | 1°6033 
1°3000 | 1°6068 


1°500 | 2°0244 


| 
| 


1°2552 | 1°6128 


29 


Wy 


y —x/¥ 


pod, 


—x/v 


1°4582 
1°4585 


1414 
1°400 
1°300 


1°9627 
1°9042 


12102 | 
11652 


16218 
(16343 


1°4607 
1°4650 


1°3466 
1°3466 
1°3469 


1-0000 | 
0°9999 | 
0°9945 


1°2911 
1°200 
1°100 
1000 
0°900 
0°800 
0°700 
0°600 
0°500 
0°400 
0°300 
0°200 
0°100 
0°000 


1°8491 
1°7980 
1°7514 
1°7097 
1°6736 
1°6437 
1°6206 
16051 
1°5978 
15994 
1°6099 
1°6281 
16470 


1°1202 
1:0754 
1:0309 
0°9870 | 
0°9439 | 
0°9019 
0°8613 
0°8227 
0°7867 
| 0°7544 
0°7271 | 
0°7068 
0°6973 


1°6507 
1°6719 
16989 
1°7322 
Ugigaalh 
1°8225 
1°8816 
1°9510 
2°0310 
2°1210 
2°2141 
2°3035 
2°3620 


14722 
1°4829 
1°4977 
1°5175 
1°5430 
15754 
1°6156 
1°6648 
1°7233 
1°7908 


| 1°8649 


1°9383 
1°9898 


13484 
1°3520 
1°3583 
1°3682 
1°3827 
1°4022 
1°4282 
14615 
175030 
1°5527 
1°6107 
1°6686 
1°7128 


0°9828 
0°9665 
0°9468 
0°9245 
0°9003 
0°8747 
0°8483 
0°8216 
0°7954 | 
0°7707 | 
0°7487 
0°7315 
0°7228 | 


1°0043 
1°0286 
1°0386 
1°0388 
1°0334 
1°0249 
10154 
1°0064 
0°9994 
0°9957 
0°9967 
1°0032 
1°0156 
10300 


1°2554 
1°2545 
1°2534 
1°2535 
1°2557 
12611 
1°2699 
1°2840 
1°3040 
1°3310 
1°3653 
1°4069 
1°4523 
1°4882 


Cae 4 


| 
1°1952| 0°7070 


y — uly 


0°6000 | 1°1783 


11181 
1100 | 0°7589 
1°0541 


0°6454 | 1°1759 


oy, 


— x/p 


—— 2 


—xly 


0°4446 
0°4639 


1°0260 


/0°2000  1°0530 


1000 
0°900 
0°800 
0°700 
0°600 
0°500 
0°400 
0°300 
0°200 
0°100 
0°000 


0°7855 
0°7982 
0°8029 
0°8029 
0°8006 
0°7977 
0°7958 
0°7964 
0°8006 
0°8092 
0°8199 


0°6706 
0°6840 
0°6897 
0°6895 
0°6854 
0°6781 
0°6686 
0°6578 
0°6468 
0°6372 | 
0°6318 | 


1°1712 
1°1670 
11645 
1°1645 
1°1681 
1°1764 
1°1902 
1°2107 
1°2378 
1°2699 
UePAS ZA 


0°5275 
0°5594 
0°5768 
0°5854 
0°5894 
0°5904 
0°5902 
0°5904 
0°5922 
0°5967 
0°6032 


1°1040 
10491 
1°0842 
1°0745 
1°0685 
1°0657 
1°0669 
1°0735 
1°0860 
11044 
1°1233 


0°2534 
0°3043 
0°3353 
0°3562 
0°3706 
0°3803 
0°3863 
0°3894 
0°3903 
0°3897 
| 0°3887 


1°0481 
1°0325 
10149 
0°9975 
0°9808 
0°9661 
0°9542 
0°9461 
0°9431 
0°9461 
0°9539 


0°1026 
0°1333 
0°1847 
0°2098 
0°2243 
0°2327 
0°2369 
0°2384 
0°2381 
0°2371 
0°2362 
0°2364 


071000 | 
0°1302 
0°1840 
0°2139 
0°2343 
0°2490 
0°2595 
0°2669 
0°2722 
0°2749 
0°2759 


0°2758 


1°0260 
1°0238 
1°0038 
0°9808 
0°9573 
0°9345 
0°9129 
0°8932 
0°8747 
0°8625 
0°8561 
0°8571 


(Hingegangen am 21. Juli 1954.) 
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Strémung in einer Expansionsdiise mit Reibung. 
Von C. Kammerer, Wien. 
Mit 7 Textabbildungen. 


Zusammenfassung. In den meisten Lehrbiichern der technischen Warmelehre! wird als 
Beispiel fur ee reibungsbehaftete Diisenstrémung die ,polytropische‘‘ Stroémung herangezogen. 
Dieser Fall wird hier eingehend behandelt und die zugehorige Form der Diise ermittelt, die sich 
aber als praktisch nicht brauchbar erweist. AnschlieBend wird die ginstigste Dusenform, namlich 
jene mit kegeligem Erweiterungsteil erértert sowie eine Beziehung zwischen dem Rohrreibungs- 
beiwert und dem Diisenwirkungsgrad abgeleitet. 


Summary. Most textbooks on heating engineering refer to the polytropic flow as an example 
for a nozzle flow with frictional resistance. This case is dealt with in detail in the following paper, 
and the corresponding nozzle form is determined, the latter being of no practical use, however. 
Furthermore the nozzle type having a conical widening, which is the most favourable form, is 


discussed, and a relation is established between the pipe friction coefficient and the nozzle 
efficiency. 


Résumé. Dans la plupart des traités de thermologie, on prend comme exemple de courant 
passant dans une tuyére et soumis & un frottement le courant «polytropique». Ce cas est ici 
traité a fond; la forme correspondante de la tuyére est déterminée, mais se révéle pratiquement 
inutilisable. On étudie ensuite la forme de tuyére la plus appropriée, c’est-a-dire celle dont la 
partie évasée prend la forme d’un céne, puis on établit un rapport entre le coefficient de frottement 
du tube et le rendement de la tuyére. 


J. Kinleitung. 


Es ist bekannt, da bei einer Laval-Diise und Nichtberiicksichtigung der 
Reibung eine bestimmte Form und Linge nicht ermittelt werden kann, wenn man 
von der EKinhaltung eines gewissen Erweiterungswinkels « absieht. Beriicksichtigt 
man jedoch die Flissigkeitsreibung z. B. nach dem tiblichen Ansatz: 

2 

Ap, =9 35 Sp (1) 
[wobei Ap, den Druckverlust durch Reibung (kp/m?), y das Raumgewicht (kg/m%), 
x die Diisenlange in Richtung der Rohrachse, 1 den Rohrreibungsbeiwert und r den 
Radius des kreisf6rmigen Querschnittes (m) bedeutet], so wird die Form und Lange 
der Diise vom Expansionsverlauf abhaingig oder umgekehrt. Der Expansionsverlauf 
[das heiBt die Abhingigkeit des Druckes, der Temperatur 7’ (° K), der Geschwindig- 
keit w(m/sec.) usw. von der Rohrlange x] laBt sich daher fiir eine gegebene Diisen- 
form ermitteln, wenn der Rohrreibungsbeiwert 2 bzw. die Wandrauhigkeit gegeben 
ist. Letztere ist bekanntlich fiir das glatte Rohr nur von der Reynolds-Zahl fiir mittlere 
Rauhigkeit von Re und dem Rauhigkeitsma8, fiir groBe Rauhigkeit jedoch nur vom 
RauhigkeitsmaB abhingig. Der Kinfachheit halber wird im folgenden mit gleich- 
bleibender Rohrreibungszahl 4 gerechnet. 


Fiir den divergenten Teil der Laval-Diise wurde eine solche Rechnung bereits 1905 
von Proell? durchgefiihrt, und zwar fiir eine lineare Abhangigkeit des Halbmessers r 
des Rohrquerschnittes von der Rohrlinge x (kegelige Diise). M. Ledinegg® legt 
eine logarithmische Abhingigkeit der Querschnittsflache zugrunde. In beiden Fallen 
wird von jenem Querschnitt ausgegangen, in dem die Schallgeschwindigkeit herrscht 
(vom sog. ténenden Querschnitt). 


1 Siehe z. B. R. Plank: Handbuch der Kaltetechnik, II. Bd.: Thermodynamische Grundlagen, 
S. 346 und 353. Berlin: Springer-Verlag. 1953. 

2 R. Proell: Z. ges. Turbinenwes. 1904, 161; 1905, 151. 

3 M. Ledinegg: Strémung kompressibler Medien in Strémungsmaschinen. Maschinenbau 


u. Warmewirtschaft H. 11 (1954). 
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In vorliegender Arbeit wird zunachst der besondere Fall der Diisenstromung mit 
polytropem Zustandsverlauf erdrtert, und die zugehdrige Diisenform ermittelt, die 
aber wegen des langen Verjiingungsteiles und schlechten Wirkungsgrades praktisch 
nicht brauchbar ist. Daher wird anschlieBend noch die fiir technische Zwecke allein 
in Frage kommende Diisenform mit kurzem Verjiingungs- und kegeligem Erweiterungs- 
teil behandelt. 


II. Allgemeine Beziehungen fiir die Diisenstr6mung mit Reibung. 


Aus dem Impulssatz* folgt die Gleichung: 


wdw+gudp + oe ay (2) 
Aus der Kontinuitatsbedingung folgt: 
pfw=pwra=GRT (3) 


[@ sekundliches DurchfluBgewicht (kg/sec.), f Querschnittsflache (m*)]. Ferner ergibt 
sich durch Differenzieren von (3): 


dp dT dw df 
piere aor Ta le (4) 


Im folgenden wird die Giiltigkeit des Gasgesetzes pv = R7', ferner konstante 
spezifische Warme c, bzw. konstantes x = c,/c, vorausgesetzt. Unter diesen Voraus- 
setzungen ist das is-Diagramm der Gase bekanntlich mit dem 7 s-Diagramm 
bis auf eine Mastabverschiedenheit identisch. 


Setzt man auBerdem Warmedichtheit der Wandungen, also adiabatische Str6mung 
voraus, so ist nach der Stromungsgleichung: 


A wdw A w : 
; =d( 29 )=-di=—c,dT (5) 
oder integriert 
Aw : 
5 aes eh tae (6) 
wobei 
Aw 
C; TARO) ra +r aii (6a) 
Aus (3) und (6) ergibt sich eine Beziehung fiir den Druck p: 
onlo.—42 
I &p 4 
p= SOYA . (7) 


Nachdem Gl. (6) auch fiir reibungsbehaftete Stroémungen gilt, so ist die Temperatur 
die gleiche Funktion der Geschwindigkeit, ob nun Reibung vorhanden ist oder nicht. 
Tragt man 7' tiber w auf, so gilt die 7'-Kurve nach Gl. (6) fiir alle méglichen Reibungs- 
verhaltnisse (Abb. 1). 


Da die Schallgeschwindigkeit bei Gasen bekanntlich nur von der Temperatur 
abhangt, letztere aber nach Gl. (6) eine eindeutige Funktion der Geschwindigkeit 
ist, so muB die Schallgeschwindigkeit stets bei derselben Geschwindigkeit auftreten, 
das heiBt ihrer GroBe nach unabhingig von den Reibungsverhiltnissen sein. Aus 
Gl. (2) wird mit v= RT/p und dp/p nach Gl. (4): 

dw RT (dT dw d A dx 
Ww 3 as if wens ; a 4r a: (8) 


* C. Kammerer: Stationére Gasstrémung durch ein gerades Rohr mit und ohne Warme- 
durchgang und Reibung. Osterr. Ingenieur-Arch. 5, H. 4, 340 (1951). 
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Ersetzt man hierin 7’ und d7'/T nach Gl. (5) und (6), ferner df/f durch 2 dr/r, so wird 
nach einigen Umformungen: 


ie R Aw 20 F A w w2 7. 
tof Cite ore (a, g) —ae 
x Ww 29 Cy, r 29¢, a 
oder 
w 2 dr 24dzx 
<7 (whe — 6) + (hot — b) = — SE = — gdy,/A (9a) 
mit 
“+1 *x—1 
dee 2% =h, gkT,=b. 
32 
37 
30 
Tit) 
| Tog “OT | | 
z: \ eygsler | fonenter ne P 
24 A Wuerschnit T)\ p 
= "ha 000 470 
coe ee | 
22+ 7 : 
27 
20 
79 | 
i 
A 
16 
75: 
a 00 200 300 400 500 600 700 600 900 mu ~" 
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Abb. 1. Druck (p), Temperatur 7 (° K) und Halbmesser des Diisenquerschnittes r (cm) in Ab- 

hangigkeit von der Geschwindigkeit w(m/sec.). Voll ausgezogen = polytropischer Zustands- 

verlauf, strichliert = verlustlose Strémung, strichpunktiert = Reibung nur im Erweiterungs- 
teil, punktiert = Diise mit zwei kegeligen Erweiterungen. 


Die Differentialgleichung (9) gibt einen Zusammenhang zwischen der Geschwindig- 
keit w, dem Halbmesser 7 des kreisformigen Diisenquerschnittes und der Diisenlange wx. 
Ist die Form der Diise also r = [ (x) gegeben, so ergibt sich aus Gl. (9) durch Inte- 
gration der Verlauf von x bzw. r tiber w, damit aber auch der Druck und Temperatur- 
verlauf iiber x. Bezeichnet man den jeweiligen Erweiterungswinkel der Diise mit «, 


so wird nach Gl. (9) 


r{dw\{ 6b 
Sessa dpe i (ae) fw ’ (10) 
aC 
oder mit Gl. (25) ; 
a t eae (/u*—o)| r 
MS aes (Bla® hha? (10a, 


wenn mit c, die spezifische Wéarme lings der Zustandsdnderung bezeichnet wird. 
Fir den ténenden Querschnitt bzw. fiir die Schallgeschwindigkeit ist: w? = gx RT,: 


2 
GRE OL I Lh Be ia Pa et 
Ra yy bzw. 7a =¢ T, 


20 
Ingenieur-Archiv VIII, 4. 
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ist. Somit wird der Ausdruck in der eckigen Klammer im Zihler gleich Null und jener 
im Nenner gleich 1/x, da c—h = 1/x ist, daher wird 
A 1 


tga = 4% oder B= Sa (10b) 


Die Schallgeschwindigkeit kann also nur im Erweiterungsteil (tg « > 0) durch- 
schritten werden und die charakteristische GroBe B betragt hierbei unabhangig 
VON 7p Oder n fiir jede beliebige Diisenform 1/x, also fiir Luft 0°715, vorausgesetzt, 
daB dw/dx endlich ,bleibt, das heift unstetiger Querschnittsverlauf ausgeschlossen 
wird. So z. B. kann w, am Ende einer konvergenten Diise oder eines geraden Rohres 
auch bei tg « = 0 oder tg « <0 auftreten. 


III. Die Diise mit polytropischem Zustandsverlauf. 


Bei polytropischer Entspannung nach der Beziehung pv" = kst ist der poly- 
tropische Wirkungsgrad nach der Beziehung?® 
(n—1) 

Npo = FRE iy (11) 
mit dem Exponenten n, der zwischen x und 1 liegt, verbunden. Anderseits ergibt sich 
der adiabat-verlustlose Wirkungsgrad 7,4 aus 7, wenn der Anfangs- und Enddruck p, 
und p, gegeben ist, aus der Beziehung 


u—1 
i (Bey * eo 
Pi 6 
Naa car x—1 3 (12) 


Bhi | 5 
; a 

Bei verlustloser Entspannung tritt die Schallgeschwindigkeit bekanntlich im 
engsten Querschnitt der Diise beim sog. kritischen Druckverhaltnis 


Ps 2. fe at 
Pr e + r) ue 
auf. Bei der reibungsbehafteten Diise ist das bekanntlich nicht mehr der Fall. Dort 


herrscht im engsten Querschnitt (f’) ein hoherer Druck p’ (p’ > p,) und das zu- 
gehorige Druckverhaltnis ist: 


p 2 n—1 
A Diya 14 
waihrend He (n ge wae 
pre ye 
DS ot 
ist. Die Schallgeschwindigkeit 
w= Vgu RT," (15) 


tritt hingegen bei einem kleineren Druckverhiltnis 


( p’ 2 Ps) 
Pi Pi 


auf, welches sich aus der Gleichsetzung: 


BIE EE / 29057 
Ws Fig = ] 


; i C. Kammerer: Wirkungsgrade bei polytropischem Zustandsverlauf in Diisen und 
Stromungsmaschinen. Maschinenbau u. Warmewirtschaft H. 5, 70 (1949). 
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ergibt, wihrend 7/7, = 2/(x +1) ist. Da der Exponent n/(n — 1) gr6Ber als 
%| (% — 1) ist und 2/(x + 1) <1 ist, so wird p” auf alle Fille kleiner als p,, das heiBt 
bei der reibungsbehafteten Diise, stetiger Querschnittsverlauf vorausgesetzt, tritt 
die Schallgeschwindigkeit nicht im engsten Querschnitt, sondern erst nach Durch- 
schreitung desselben, also bereits im erweiterten Teil der Diise auf. 


Durch den EinfluB der Reibung wird also gegeniiber dem reibungslosen Fall der 
engste Querschnitt zu einem héheren Druckverhiltnis und der ténende Quer- 
schnitt (f’’) zu einem niedrigeren Druckverhiltnis hin verschoben. Die mit (’) und 
mit (’’) bezeichneten Querschnitte wandern also in entgegengesetzte Richtungen. 
Die Geschwindigkeit im engsten Querschnitt w’ betriigt fiir die reibungsbehaftete 
Diise nach Gl. (6): 


OP La 2 gre, Dey 2 \ 
w =|) A aaa) (18) 
und die Schallgeschwindigkeit : 
Bt 2o¢,2, fee) 
an V eat L =1} ba 


Ersetzt man hierin c,/A durch R~x/(x — 1), ferner 7, durch Ta so wird: 
w" =VgxBT,, 

das hei®Bt die Schallgeschwindigkeit ist fiir die reibungslose und reibungsbehaftete 

Diise von gleicher GréBe, wie bereits aus Abschnitt II hervorging. Die Geschwindig- 


keit im engsten Querschnitt (w’) ist natiirlich kleiner als die Schallgeschwindigkeit w,, 
und zwar ist das Verhaltnis nach Gl. (18) und (19): 


On) (aa) (a1) 
a a ea iY 
Fiir die polytropische Diisenstromung ist der sog. polytropische Wirkungsgrad: 
Aw 
: Adit a 29 ) (21) 
Poll nA Rete Ph gp ris v4 Ane 
(ay) + 


[L.4 wirklich geleistete Arbeit (kpm), g, Reibungswarme (kcal)] eine fiir den ganzen 
Zustandsverlauf konstante GroBe>. Nach Gl. (1) ist: 


Ap, w* 4 dx . 
dq, = - = ae (22) 
Dieser Ausdruck, in Gl. (21) eingesetzt, ergibt die Beziehung: 
1 d(w?) Ada 1 2 dw 
= - = ——— 1 . y 
Sa 1) we 27 ES Ww (23) 
Nach Gl. (11) ist: 
il 1 \ = “x—n 
aes x (n—1) 
Ferner ist bekanntlich: 
Se (os) 
Cn Sr Cy (n— 1) ) 


Mit der Substitution 
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erhalt man schlieBlich die einfache, allgemein giiltige Gleichung: 
d 
2k, (S) == Sderot (=) = : (25) 


AC Nee r 


in welcher k, baw. c, fiir die Polytrope eine Konstante ist. Gl. (25) verkniipft die 
Geschwindigkeit w mit der Diisenlange x und dem Halbmesser 7. Anderseits lat 
sich y nach Gl. (3) auch schreiben: 


’ ao (3a) 


mW p 


Da fir die Polytrope 


1 
(po =p, 0")... = 0,(24)", 


= arta ral 
oe 2 (2) 
ek ee neta : (26) 


Tt WwW 7 W 


so wird mit T/p =v/Rk 


‘ Nun ist fiir polytropische Expansion die Geschwindigkeit : 


a 
Se oot re 
= 1 27 
v= | 224% [1-(2 7) 
n 

rae (1 a he 

Pr 29 Cy Ty 
Dies in Gl. (26) eingesetzt, liefert: 

(liz ae rae (28) 


oder 


(27a) 


Diese Gleichung gestattet, den Halbmesser des Diisenquerschnittes in Abhangigkeit 
von der Geschwindigkeit zu berechnen. 


Setzt man 7 aus (28) in (25) ein, so erhalt man schlieBlich folgende nur w und x 
enthaltende Differentialgleichung : 


G if 
2 ky / — 3 is fe =k (29) 
2 A w? oa 
Vie |e 
29 Cy 4BF 


die sich aber in geschlossener Form nicht integrieren liBt. Den kleinsten Halbmesser r 
des Diisenquerschnittes erhalt man durch Einsetzen von w’ aus GI. (18) in Gl. (28) zu: 


Gv, 1 

min ~~ se | af . 

4 eS n—l 2 (eet 

Va Vaaa 

Der veranderliche Erweiterungswinkel der Diise mit polytropischem Zustandsverlauf 

betragt nach Gl. (10a), (25) und (24) 

ho oe (b/w? —c) 
8 (b/w? — h) 

Beispiel: Es ist eine polytrope Diise fiir Luft von 10 ata Anfangsdruck, 1 ata Gegen- 


druck ¢t, = 600°C und y,4= 0°7 zu ermitteln. 4 wird als konstant mit 0°04 an- 
genommen. Dann wird: igh; 


Nya = 0°623, hy = 3025, und | = 1-216 


r (28a) 


tga= (30) 
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Mit R = 29°27 wird v, = 0°2555 m3/kg. r laBt sich nach GI. (28), 7 nach Gl. (6) 
und p nach (3) als Funktion von w ermitteln (Abb. 1) und Tabelle 1. (Zum Vergleich 
sind in Abb. 1 auch die betreffenden Werte fiir die reibungslose Diise eingezeichnet.) 
Um die Rohrlinge «2 zu erhalten, kann man Gl. (25) beniitzen und schrittweise fiir 


Aw z. B. = 20 oder 10 (m/jsec.) das 
zugehorige Ax ermitteln, da der zu 
jeder Geschwindigkeit w gehdérige 
Radius r bereits bekannt ist*. So 
ergibt sich r als Funktion von z, das 
heiBt die gesuchte Form der poly- 
tropen Diise (Abb. 2). Fiir kleine 
Werte von w (etwa w < 40 m/sec.) 
ist nun /z. B. (= 0°04) gegen den 
Klammerausdruck im Zahler der 
Gl. (30) zu vernachlassigen. Man 
erhélt dann, weil die Ausdriicke c 
und / ebenfalls zu vernachlassigen 
sind, néiherungsweise: 


een a 


R 
N 
Pa 
S 
Py 
CN 
N 
i 
SS 


Wshistebe aur 
30+ 70 verkleiner? 


ES engs 


20-933 
Zz(m) 


- 100 -$0 pe a Md 


Abb. 2a. Diise mit polytropischem Zustandsverlauf 
und 70% Wirkungsgrad (¢, = 600°C, p, = 10 ata, 
P, = 1 ata). 


als konstant. Der Verjiingungsteil der polytropischen Diise nimmt daher fiir kleine 
Geschwindigkeiten die Form eines geraden Kegels an (also gerade umgekehrt wie bei 
der normalen Diise). Zahlt man die Diisenlange « vom engsten Querschnitt, also von 


engsier — ronenger 750 


o 
Lisenlénge z{(Merer) 


200 +2 


=) 
70 00 -7 


Abb. 2b. Druck- und Geschwindigkeitsverlauf fiir die polytrope Diise Abb. 2a itber der Diisen- 


lange x. 


w — 414 m/jsec. aus, so ergeben sich die «-Werte nach Tabelle 1. Bei weiterer Ver- 
kleinerung von w unter 0°1 m/sec. bewegt sich sowohl r als auch x gegen oo. Der 


* Cenauer ist der Verlauf von w iiber 2 durch Reihenentwicklung aus Gl. (29) VAVN erhalten. Fur 
das vorliegende Beispiel wurde diese Entwicklung von Dipl.-Ing. Mayr-Harting durchgefuhrt, 


wofiir ihm an dieser Stelle gedankt sei. 
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konvergente Teil wird linger als der Erweiterungsteil, letzterer konvex nach innen. Bei 
einer solchen Diise tritt bereits im Verjiingungsteil ein erheblicher Reibungsverlust auf. 


Tabelle 1. 
a ere ee 
| p (ata) | 1A Ga 6) | w (m/sec) | 7 (em) | a (m) 

IBehialteriahys atcalsciceee recites 10 873 0 co 6S 

10 873 01 90°25 — 107°47 

10 873 1:0 28°50 — 32°89 

10 873 10 9°10 — 9°33 
Engster Querschnitt ...... 5°75 788 414 1°79 + 0 
Ténender Querschnitt ..... sri 730 540 1°875 + 0°293 
Endquerschnitt........... 1:00 450 765 2°145 + 0°747 


Das adiabat-verlustlose Warmegefille betrigt (873 — 450). 0°24 kcal/kg, die 
Temperatur beim engsten Querschnitt 790° K, daher ist der Verlust bis zum engsten 
Querschnitt (873 — 790)/(873 — 450) = 83/423 = 19°6% oder 19°6/0°3 = 65°5% des 
Gesamtverlustes, also unnétig groB. Auch die konvexe Form des Erweiterungsteiles 
ist ungiinstig, weil der giinstigste noch zulassige Erweiterungswinkel von zirka 4° 
nur an einer Stelle, nimlich am Diisenende angewendet werden koénnte. 

Fir praktische Zwecke kommt bekanntlich gerade umgekehrt nur ein kurzer 
Verjiingungsteil und ein Erweiterungsteil in Frage, der an allen Stellen den giinstigsten, 
das heiBt gro8tméglichen Erweiterungswinkel « aufweist. Da dieser Winkel geschwindig- 
keitsunabhingig ist, ergibt der gerade Kegel den besten Wirkungsgrad bzw. die beste 

Diisenform. Wiirde bei kleinen Uberschallgeschwindigkeiten ein 
eroBerer Erweiterungswinkel als etwa 4° zulassig sein, dann ware 
etwa eine nach innen konkave Form zu bevorzugen. 


IV. Die kegelig erweiterte Diise. 
1. Der Ubergangsteil. 


Abb. 3. Die Schallgeschwindigkeit tritt nach Gl. (10a) immer dort auf, 


Diise mit unstetigem 
Querschnittsverlauf. 


wo die GréBe B = A/8 tg « von einem Wert unter 1/x (= 0°715 
fiir Luft) in einen solchen iiber 1/x% entweder stetig oder unstetig 
tibergeht. Fiir konstantes 2 mu8 dabei der Tangens des Erweite- 
rungswinkels den Wert 4/8 B iiberschreiten. Dieser sog. ténende 
Querschnitt kann bei unstetigem Ubergang, wenn also die Kontur 
der Diise einen Bruchpunkt aufweist, auch mit dem engsten Quer- 
schnitt zusammenfallen (z. B. Abb. 3). Bei stetigem Verlauf der 
Diisenbegrenzung mu8 er jedoch in einem gewissen Abstand 42, 
nach dem engsten Querschnitt liegen, wobei es auf den Aus- 
rundungsradius des Ubergangsteiles ankommt. Eine Ausrundung 
ist schon deshalb erforderlich, damit die Strémung sich von der 
Wand nicht ablost. Ist z. B. beim kleinsten Durchmesser der 
Ausrundungsradius @ (Abb. 4), so betragt der Abstand des toénen- 
Abb. 4. den vom engsten Querschnitt: 


Diise mit kegeligem 
Erweiterungsteil und 


stetigem und mit tg «, = 1/8 B= 1x/8 
Querschnittsverlauf. 


Ary = 0 sin x, 


Ax, = osinarc tg (Ax/8) = 0 Ax/8. 


So ist z. B. Ax, fiir eine Diise mit 3m engstem Durchmesser, 9 = 1 em, x = 1°4, 
A = 0°04, Ax, = 0007 cm = 0°07 mm = 70 4 (70 Mikron). 
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Der Abstand des ténenden vom engsten Querschnitt ist demnach so gering, dah 
er fiir praktische Rechnungen auBer Betracht bleiben kann. 


Der geradkegelige Teil der Diise beginnt erst im Abstand 
Aut O'sil os 
vom engsten Querschnitt, also z. B. fiir a, = 4° (tg « = 0°07) und 9 = 1em im 


Abstand von 0°697 mm. 


2. Der kegelige Erweiterungsteil. 


Die allgemein giiltige Gl. (9) laBt sich mit dx = dr/tg ~ auch in folgender Form 
schreiben: 


1 AR 
ae: || ee == 
2dr dw Alea ste gre, ss 
a es +42) RT)’ ) 
8 tga Oe. nae 


wobei tg « nun als konstant betrachtet wird. 


Fiihrt man der Einfachheit halber folgende Substitutionen ein, 


1 AR x+1 


x oF 2 Cy way Porta 2 
yh AR x—l1 
8tg x j 2¢,» =e Bi Deere 
und 
Gly Ca Oe), 
so wird: 
2dr dw{w?c—b|_ dw |w?—bd/c| ¢ 
ii elas 6 | = w = Gh (31a) 
Integriert, ergibt Gl. (31) nach einigen Umformungen: 
C a—c 
(ie (a 1b) nt (33 «) 
w Vx 
oder 
Oo} .3 
a ieee oe (33 f) 


wobei, wie leicht einzusehen ist, Fall « fiir Unterschall (w? > b/a und a >), hin- 
gegen Fall B (w? <b/a und a <c) fiir Uberschall gilt. Fiir die Schallgeschwindig- 
keit ist nach Friiherem w,2 = b/c. Damit nun r mit w zunimmt (tg « > 0), miissen 
Zahler und Nenner der rechten Seite von Gl. (31a) ungleiche Vorzeichen haben. 
Daher muB gelten: 


entweder bla <u* =< bic Fall «, Unterschall, 
oder bla > w® > b/c Fall 6, Uberschall. 


Im ersten Falle tritt die Schailgeschwindigkeit w, am Ende der kegeligen Erweiterung, 
im zweiten Falle (8) am Beginn derselben auf. Aus Gl. (31a) geht ferner hervor, 
daB die Funktion r? = f (w) fiir w? = b/c, das heiBt fiir die Schallgeschwindigkeit 
ein Maximum oder Minimum aufweist (waagrechte Tangente der Kurve r tiber w in 
den Fallen « und f). 


302 C. Kammerer: 


Man kinnte z. B. an den engsten Querschnitt der polytropischen Diise, also bei 
Unterschall, einen Kegel bis zur Erreichung von w, (Fall «) und beim ténenden Quer- 
schnitt, also bei w,, einen zweiten Kegel mit gréBerem Erweiterungswinkel (Fall 6) 
ansetzen, um auf Uberschallgeschwindigkeiten zu kommen. Siehe Abb. 5 und die 
punktierte Kurve in Abb. 1. 


Fiir den ersten Kegel mu8 nach Friiherem b/a < w, bzw. a > b/w? im Beispiel 
also gréBer als 250500/171500 = 1°46 sein, damit a w* — 6 eine positive GroBe ist. 
a konnte auch beliebig groB gewahlt werden. Z. B. co (tg « = 0, gerades Rohr, s. 
FuBnote 1). Wahlt* man z. B. a= 1643, so wird B=a— — == 15 une 
tg « = 1/8 B = 004/12 = 0°008, das hei®t 333mm pro m. Ferner wird b= 
— 9°81 - 29°27 - 873 = 250500, c = (x + .1)/2 x = 0°857, (a — c)/2 a = 0°786/3°286 = 
= 0°239 und die Konstante C, = 143°5- 107%. 


Damit wird fir 


w, = 414 m/sec. fT, SSO cm p = 5°55 ata 
w == 500%, Yized Gea, T= oD hhg 
16, == DAO 7p Mi 99s qo NOD = ebro. 


Die Linge des ersten Kegels betragt 60°1 cm. 
Eine héhere Geschwindigkeit als w, ist mit diesem ersten Konus nicht erzielbar. 


Fiir den zweiten Kegel (Fall 8) mu8 nach Friiherem b/a > w,Z bzw. a < b/w, 
gewahlt werden. 


Wegen w = b/c muB also a <c oder B < 1/x (z. B. fiir Luft < 0°714) sein. Die 
mit einem bestimmten Wert a erreichbare Hochstgeschwindigkeit ist w?,,, = b/a, 
da sonst der Klammerausdruck in Gl. (35) negativ wird. Setzt man z. B. w”,,,, fiir 
das vorliegende Beispiel mit 750 m/sec. fest, so wird a = b/w*,,,x = 250500/562000 = 
= 0°446, B= a — 0143 = 0°303, c — a = 0°857 — 0°446 = 0°411 und (c — a)/2a = 
= 0°461. 

Wirde man w”,,, = w, wahlen (keine Geschwindigkeitssteigerung), so ware 
a=c bzw. B = 1/x und der Exponent aber auch der Klammerausdruck in Gl. (33) 
wird Null. Beispiel: Mit w, = w, = 540 m/sec. und r, = 1:99 cm wird die Konstante 
C, = 577.540 - 10-4. 


Damit ergeben sich folgende Halbmesser und Driicke: 


Fir w = w, = 540 m/sec. r = 1°99 cm p=319 ata 
WU ==000 oh; Fe 202i5, pp 2 BCS 
wie= TOO. he 2385, p=Trb2T. y 
RU Ty. Sy fee Deen. p == 0987 - 5s 


Die Linge des zweiten Kegels betraigt 50°3 cm. 


3. Ein Sonderfall ist jener, wo a=c bzw. B= 4/8 tg « = 1/x gewahlt wird. 
Dann geht (33 «) oder (33 6) (den Fall w = w, ausgenommen) iiber in: 


C4 
wVa 


= 


oder 


ips 


(2) bzw. 7?w—=G@-v- /x = konstant. (34) 
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Da tg « positiv sein muB, gilt Gl. (34) nur fiir den Erweiterungsteil. Fiir den Er- 
weiterungswinkel « = arc tg x 4/8 nimmt die Geschwindigkeit nach Gl. (34) mit zu- 
nehmendem r wieder ab. Man gelangt sonach von der Schallgeschwindigkeit auf 
Unterschall. Wegen r?2w, = 7,2 w, =Gv/x ist fiir diesen Sonderfall nach Gl. (34) 
das spezifische Volumen v bzw. das Produkt r?w konstant. Nachdem aber w in 
Strémungsrichtung abnimmt, mu8 nach Gl. (7) der Druck (p) zunehmen. Es handelt 


sich daher um einen Unterschalldiffusor bei konstantem Volumen, beginnend bei 1,. 
Nach Gl. (7) ist: 


G (b—hw 
pao ba hw) 


g ur * 


(34a) 
wobei 7?w als konstant zu betrachten ist. 


Um auf Uberschall zu kommen, muB offenbar « groBer als vorhin bzw. B <1 |x 
gewahlt werden, wie bereits fiir den zweiten Kegel des vorhergehenden Beispieles 


ei ys : ; : 
angefiihrt wurde. « = arc tg bildet die untere Grenze des Erweiterungswinkels, 


bei welchem eine GeschwindigkeitserhOhung méglich ist. Bei dieser unteren Grenze, 
also bei B = 1/x = 0°857 ist, ausgehend von w,, keine Geschwindigkeitserhohung 
mehr moglich. Die Schallgeschwindigkeit bleibt bei sinkendem Druck langs des ganzen 


Kegels erhalten (Gi. (34) gilt fiir diesen Fall nicht, da 2,“ — (F)° das heiBt 
2 2 

unbestimmt wird). Es ergibt sich also eine isotherme Entspannung bei gleichbleibender 

Geschwindigkeit und Temperatur (vollstiindige Drosselung), worauf bereits Ledinegg* 

verwiesen hat. Hierbei ist der Druck p wegen Gl. (3) umgekehrt verhiltnisgleich, 

r? bzw. das Produkt p-r? = konstant. 


Ausgehend von der Schallgeschwindigkeit sind also fiir eine kegelige Erweiterung 
mit B= 1/x zwei Moglichkeiten vorhanden, und zwar, wenn der Gegendruck hoher 
ist, Drucksteigerung und Geschwindigkeitsabnahme, wenn der Gegendruck kleiner 
ist, isotherme Expansion bei konstanter Geschwindigkeit und Druckabnahme. 


4. Druckverlauf bei verinderlichem Erweiterungswinkel. 
Nach Gl. (7) und (33) wird fiir Unterschall (Fall «) 
GVa (b—hw*) 


Bis ae ew) 
xg CO» (aw*—b) 2a 


und fiir Uberschall (Fall f): 


Ga (b—hw?) (b— av?) 2a 
p= SAGs 3 (34c¢) 
Fiihrt man Verhaltniswerte ein, so ist fiir Fall « bzw. nach Gl. (34b) 
pr, (b— hw?) (awP—b\ 2 
Poe eee aa a 
und fir Fall B bzw. nach Gl. (34c): 
1 b—hw,?\/b—aw,?\ 24 . 
pe aatieeat)” pe 


Nach diesen Gleichungen la&t sich das Verhalten des Druckes berechnen, wenn man 
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z. B. von der Schallgeschwindigkeit w, bei p, und 7’, ausgeht. So erhalt man beispiels- 
weise fiir den unter IV, 2 am SchluB angefiihrten Ausgangspunkt: 


w, = 540 m/sec., p, = 319 ata, r, = 1:99 cm, 

verschiedene Kurven, je nachdem tg « bzw. B= 4/8 tg « gewahlt wird (Abb. 7). 
Gezeichnet sind die Faille B = 1/x (schon unter IV, 3 behandelt), ferner B = 0°303, 
B= 0 (A= 0 oder tg a = oo), B= — 0°786, B= — 1 und B = — 1000 (das heit 
tg « nahezu Null, also dem geraden Rohr entsprechend) (bereits im Aufsatz* ausfiihr- 
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Abb. 5. Diise mit polytropem Verjiingungs- und doppelkegeligem Erweiterungsteil. 


lich behandelt). SchlieBlich ist noch der Fall B = + 1°5 eingezeichnet (unter IV, 2 be- 
handelt bzw. auch in Abb. 5 gezeigt). Der médgliche Durchlaufungssinn ist jeweils 
durch Pfeile angegeben. Rechts oben und links unten sind keine Verlaiufe méglich 
bzw. befinden sich nicht durchschreitbare Gebiete. Die Bereiche der konvergenten 


Kegel sind durch Schraffur gekennzeichnet. Die Kurve B = 0 entspricht der reibungs- 
losen Strdmung, fiir welche bekanntlich 


Pi/P, = (T/T 2)? = = ” Ww? 


ist. Die Kurve B= 
Durchmessers. 


+ co entspricht wegen tg « = 0 dem geraden Rohr gleichen 


5. Fir praktische Zwecke kommt nur eine Diise mit kurzem Verjiingungsteil, 
abgerundetem Ubergangsteil und kegeliger Erweiterung nach Fall 6 in Frage. Im 
Gegensatz zur polytropen Diise ist der Verlust bis w, sehr gering. Der Hauptverlust 
tritt hier im kegeligen Erweiterungsteil auf, weil dessen Erweiterungswinkel wegen 
der Ablosungsgefahr nicht groBer als etwa 4° (tg « = 0°07) gewahlt werden darf. 
Praktisch kann man daher den konvergenten Teil der. Diise als verlustlos annehmen 
und die Rechnung auf den kegeligen Erweiterungsteil beschranken. Dann fallt auch 
der tonende mit dem engsten Querschnitt zusammen. 
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Legt man wieder eine Anfangstemperatur von 600° C, 10 ata im Behalter, 1 ata 
Gegendruck und, wie iiblich, einen adiabaten Wirkungsgrad von 94% zu- 
grunde, so wird w, = 892 m/sec., 7’, = 475°5° K, v, = 1°39 m3/kg und r, = 2:226 cm. 


Es entsteht nun die Frage, wie groB der Reibungsbeiwert / sein darf, damit sich 
der vorhin angenommene Wirkungsgrad von 94% beim ausfiihrbaren Héchstwert 
von tg « (= 0°07) ergibt. 


Nach Gl. (33 £) ist fiir kegelige Erweiterung im Uberschallgebiet: 


oes: 
te W, (b—aw,2\ 24 © 
Toh 0) Wy be —aw,? (35) 
oder 
Tre MOD 6==4 
log (Fr aot) = “aq Log (b — a we) — log (b — aw]. (35a) 
on 
\020 
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Abb. 6. B= A/8tga bzw. A in Abhangigkeit vom adiabaten Wirkungsgrad 7,4 fiir Luft, 
t, = 600° C, p, = 10 ata und tg « = 0°07 (a = 4°). 


Mit c = 0°857, b = 250500 laBt sich a berechnen, wenn beachtet wird, daB a < b/w? 
sein muB. Durch wiederholtes Probieren findet man a = 0°2027 und kann sonach 
den Halbmesser 7 auch fiir Zwischenwerte von w bestimmen. Auf die gleiche Weise 


ergibt sich 


fin 4yg2 (0°94 voce docs et a = 02027, B= 0:0597, 
Ors eee a = 02400, B = 0-097, 
ii ey a = 0°2850, B = 0°142, 
2 ae a = 0°3250, B= 0°182, 


da B=a — 0°143 ist. 


Fiir ein anderes Druckverhiltnis p,/p, ergeben sich andere Werte von B, die 
fiir p,—05ata in Abb. 6 ebenfalls eingezeichnet sind. Fiir andere Anfangs- 
temperaturen als 600°C miissen die a-Werte neu bestimmt werden. Man entnimmt 


aus Abb. 6, da z. B. fiir 7,g= 0°94 und tg « = 0°07 
1 = 8tg a: B, also 8- 0°07 + 0°0597 = 0:0334 


betraigt. Ist 2 kleiner, z. B. 0:0116 fiir die glatte Diise, so sind noch hohere Wirkungs- 
grade erzielbar. Bei Gegendriicken gréBer als 1 ata bzw. bei Druckverhaltnissen 
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Po/p, > 01 muB der Reibungseinflu8 auch des konvergenten Teiles beriicksichtigt 
werden. 
V. Die kegelig verengte Diise (B < 0). 
Liegt der Vorgang im Unterschallgebiet, so gilt Gl. (33 «). Wird der Anfangs- 
radius mit 7, und der Endradius mit r, bezeichnet, so besteht nach Gl. (33 «) die 


Beziehung : 
a—e 


Caladium yale (36) 


Ts Ww.) | (b—a w,2) J d 


wobei nun tg « als negative GréBe einzusetzen ist. Die GréBe B= 1/8 tg « wird 
dann auf alle Falle negativ, die GréBe a aber wegen a = B +h nur dann, wenn 
B>h ist. 


oe 2 . ae y, Dw-Liegrann 
= = + + 
ce me Bg ‘ y E53 fonvergener Teil 
iN L- p=£ -9775 
Ble +04 SSaA4 Y a 
SS U 
SS y nicht durchsthreveares 
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4 201 = + & SN y. 
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Abb. 7. Druckverlauf tiber der Geschwindigkeit (p w-Diagramm) fiir die kegelig erweiterte oder 
verengte Diise. 


Ist z. B. r, = 0°9025 m, w, = 0°1 m/sec., pj = 10 ata (wie im Beispiel fiir die 
polytrope Diise), ebenso tg « = 0°00825, A = 0°04, so wird B = — 0°606, h = 0°143, 
somit a = — 0°462 und (a —c)/2a = + 1°425. 


Nach Gl. (36) wird dann 


fur Wy 300 Ms Se0 5 rw eae f= 187) om 
We SOU 35) Donatos fas L785", 
bbs <4) Sy NE, Be ro = 1715 ,, 
MeeOLO + saps tnie See ae r, = 1675 ,, | p, = 4°485 ata 


t 


Bei einer solechen Diise wird w, schon bei einem wesentlich héheren Druck erreicht 
als bei der polytropen Diise, was sich daraus erklirt, daB die kegelig verengte Diise 
um zirka 1:2 m kiirzer ist als die polytrope (bis w, gerechnet), somit einen kleineren 
Druckabfall aufweist. 

(Eingegangen am 30. September 1954.) 


Persénliches. — Buchbesprechungen. 307 


H. Mache gestorben. 


Hofrat Prof. Dr. Heinrich Mache, langjahriges Mitglied des Professoren-Kolle- 
giums der Technischen Hochschule Wien und Mitglied der Akademie der Wissen- 
schaften, starb am 2. September 1954 im 
Alter von 78 Jahren. reeee 

Mache studierte bei Mannern wie Boltz- 
mann, F. Exner und Mach, promovierte 
1898, habilitierte sich 1901 fiir Physik an der 
Universitat Wien, 1906 wurde er auBer- 
ordentlicher Professor an der Universitit 
Innsbruck, 1908 an der Universitat Wien und 
schlieBlich 1911 ordentlicher Professor an der 
Technischen Hochschule Wien, wo er bis 
1947 wirkte. Er war mehrmals Dekan und 
1925/26 Rektor. 

Mache war als Physiker weit iiber die 
Grenzen Osterreichs hinaus bekannt und 
geachtet, was in zahlreichen Ehrungen zum 
Ausdruck kam. Er hat unter anderem an 
der Erforschung der Radioaktivitaét erfolg- 
reich mitgearbeitet. Neue Erkenntnisse auf 
dem Gebiet der inneren Ballistik sind ihm 
ebenso zu verdanken wie zahlreiche Ver- 
dffentlichungen auf verschiedenen Gebieten 
der Physik. Besonders beschaftigten ihn die 
Verbrennungserscheinungen, die Struktur 
der Flammen, wie tiberhaupt die Theorie der Warme, wie sein ausgezeichnetes 
Lehrbuch ,,Theorie der Warme“ vom Jahre 1921 sowie seine Schrift ,, Vorlesungen 
tiber Warmelehre‘‘ vom Jahre 1944 bezeugt. 

Als Experimentalphysiker war Mache der Forschung im Laboratorium, als Professor 
aber auch der Lehre im Horsaal und der Sorge um die Studierenden hingegeben, 
denen er in meisterhaftem Vortrag sein Wissen vermittelte. Seinen Schiilern, Mit- 
arbeitern und Kollegen war er ein hilfsbereiter Freund und Berater, an den sie sich 
stets mit Hochschatzung und Dankbarkeit erinnern werden. 


C. Kimmerer, Wien. 


Buchbesprechungen. 


Mathematics for Students of Engineering and Applied Science. Von L. B. Benny. Mit 177 Textabb., 
VII, 7838S. London-New York-Toronto: G. Cumberledge, Oxford University Press. 1954. 
Geb. 35s. net. 


Das vorliegende Werk mu im Hinblick auf die Eigenart des englischen Hochschulwesens 
gewiirdigt werden. Es ist als Behelf fiir die Studierenden der technischen Wissenschaften und 
der Naturwissenschaften gedacht, die ein ,,Higher National Certificate‘ anstreben und stellt 
daher in gewissem Sinne eine Fortsetzung der vom gleichen Verfasser stammenden Biicher 
,Mathematics for students of technology”, ,,Junior Course“‘ und ,,Senior Course‘ dar. Dem- 
gemaé® kénnen daher sehr rasch bedeutende Fortschritte in der Erarbeitung des Stoffes gemacht 
werden, wobei physikalische Gesichtspunkte oftmals im Vordergrund stehen. Die Darstellung 
ist klar und iibersichtlich und verliert sich niemals in unwesentliche Details, fiir die sich aus der 
Zielsetzung des Werkes keine Rechtfertigung ergibt. Soweit die Beweise durchgefuhrt werden, 
griinden sie sich immer auf genau formulierte Voraussetzungen in der jeweils fur notwendig 
gehaltenen Allgemeinheit ; der Verfasser scheut sich jedoch nicht, oftmals auch Satze ohne Beweise 
anzufiihren und sich darauf zu beschranken, die Erfiillung der Voraussetzungen fiir die Giltigkeit 
an Beispielen zu illustrieren. Die Fehlerhaftigkeit mancher Figuren, die raumliche Gebilde 
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darstellen, kann dadurch erklart werden, daB die Pflege der konstruktiven Methoden der 
Geometrie im englischen Sprachgebiet etwas vernachlissigt wird. Eine tiberaus umfangreiche 
Sammlung von Ubungsaufgaben erméglicht es dem Leser, sich die routinemabige Beherrschung 
des Stoffes anzueignen. Der reiche Inhalt wird jedoch dem Werk auch auBerhalb Englands die 
verdiente Verbreitung bringen. Das Werk gliedert sich in die Kapitel: Limits-convergence of 
series —- Binomial, exponential and logarithmic series hyperbolic functions; Imaginary and 
complex quantities — Moivre’s Theorem — Algebraic equations; Determinants; The elements 
of plane coordinate geometry; The elements of solid coordinate geometry; Vectors and vector 
products; Differentiation and its applications; The expansion of functions in series — The 
Theorems of Taylor and Maclaurin; Integration; Applications of integration ; Differential equations ; 
Functions of more tham one independent variable — partial differentiation and its applications; 
The integration of linear differential equations in series; Fourier Series — Harmonie Analysis; 
Partial differential equations; Multiple integrals and their applications; Spherical Trigonometry ; 
Differentiation of Vectors; Errors of observation — Method of least squares; Miscellaneous 
Exercises; Answers. R. Inzinger, Wien. 


Handbook of Elliptic Integrals for Engineers and Physicists. Von P.#. Byrd und M. D. Fried- 
man. (Die Grundlehren der mathematischen Wissenschaften in Einzeldarstellungen mit be- 
sonderer Beriicksichtigung der Anwendungsgebiete. Herausgegeben von R. Grammel, E. Hopf, 
H. Hopf, F. Rellich, F. K. Schmidt, B. L. van der Waerden: Band LXVII.) Mit 22 Textabb., 
XIII, 3558. Berlin-G6ttingen-Heidelberg: Springer-Verlag. 1954. DM 36.—, geb. DM 39.60. 


Der Techniker und Physiker ist an einer raschen Auswertung der bei seinen Untersuchungen 
auftretenden Integralen sehr interessiert, weshalb er es als lastig empfindet, die mathematische 
Literatur zur Durchfiihrung oftmals umfangreicher Substitutionen heranziehen zu miissen. Dies 
ist der Grund, weshalb sich Integraltafeln in letzter Zeit einer stets zunehmenden Beliebtheit 
erfreuen. Ein besonderes Bediirfnis lag fiir die Anlage einer umfangreichen Tafel elliptischer 
Integrale vor, das durch das vorliegende Werk in jeder Hinsicht voll befriedigt wurde. Die Ver- 
fasser stellten sich die Aufgabe, ein Handbuch der elliptischen Integrale zusammenzustellen, 
dessen Anwendung keine wesentlichen Kenntnisse aus der Theorie voraussetzt. Hs enthalt rund 
3000 Integrale und Formeln in der Legendreschen und Jacobischen Bezeichnung, die den meisten 
praktischen Beditirfnissen in vollem Umfange gentigen diirften. Daneben werden viele Trans- 
formationen und Rekursionsformeln angegeben, die es gestatten, auftretende Integrale auf die in 
der Sammlung enthaltenen Typen zuriickzufiihren. Eine tibersichtliche Anlage des Stoffes sowie 
ausfithrliche Erklarungen und wechselseitige Verweisungen erleichtern die Handhabung des Werkes. 
Fur numerische Auswertungen stehen Zahlentafeln fiir die Integrale der ersten und zweiten Gattung 
in der Legendreschen Normalform sowie ftir die Jacobische Zetafunktion und die Heumansche 
Lambdafunktion zur Verfiigung. Inhaltsverzeichnis: Definitions and Fundamental Relations; 
Reduction of Algebraic Integrands to Jacobian Elliptic Functions; Reduction of Trigonometric 
Integrands to Jacobian Elliptic Functions; Reduction of Hyperbolic Integrands to Jacobian 
Elliptic Functions; Tables of Integrals of Jacobian Elliptic Functions; Elliptic Integrals of the 
Third Kind; Table of Miscellaneous Elliptic Integrals Involving Trigonometric of Hyperbolic 
Integrands; Elliptic Integrals Resulting from Laplace Transformations; Hyperelliptic Integrals; 
Integrals of the Elliptic Integrals; Derivatives; Miscellaneous Integrals and Formulas; Expansions 
in Series; Appendix; Bibliography; Index. R. Inzinger, Wien. 


Model Analysis of Structures. Von 7’. M. Charlton. (Spon’s Civil Engineering Series. Heraus- 
gegeben von W. Fisher Cassie.) Mit 58 Textabb., XII, 142 S. London: E. & F. N. Spon 
Limited. 1954. Geb. 21s. net. 


Ks ist sehr verdienstlich, daB es der Autor unternommen hat, Verfahren fiir Modellversuche 
zur Bestimmung der Beanspruchung statisch unbestimmter Baukonstruktionen in Buchform 
zusammenzufassen. Behandelt werden darin vornehmlich Fachwerke und Rahmen, bei denen 
ein linearer Zusammenhang zwischen Belastung und Verformung besteht, doch wird auch auf 
den Fall eines nichtlinearen Zusammenhanges eingegangen. Hervorzuheben sind ferner die Aus- 
fuhrungen uber die Wirtschaftlichkeit von Modellversuchen. 

Das sehr klar geschriebene und gut ausgestattete Buch verdient daher die Beachtung jedes 
Statikers. E. Czitary, Wien. 


Technische Strémungslehre. Von B. Eck. Vierte, verbesserte Auflage. Mit 407 Textabb., 
X, 4228. Berlin-G6ttingen-Heidelberg: Springer-Verlag. 1954. Geb. DM 29.40. 


Ein Werk, das in vierter Auflage erscheint, ist bestimmt so bekannt und beliebt, daB es 
nicht weiter empfohlen werden mu. Der Verfasser hat aber die neue Auflage auch noch etwas 
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erweitert und verbessert. In der Gasdynamik werden das Ranquesche Wirbelrohr und die 
A-StoBe besprochen. Auch sonst findet man Erginzungen und Verbesserungen. Auch wurde 
die Anzahl der Abbildungen vermehrt, womit dem Wunsche mancher Leser entsprochen wurde. 
Die Strémungslehre von Eck wird wegen ihrer Reichhaltigkeit und erstklassigen Ausstattung 
wieder vielen Ingenieuren willkommen sein. Rk. Bruniak, Wien. 


Aufgaben aus der Hydromechanik. Von K. Federhofer. 245 Aufgaben nebst Lésungen. Mit 
235 Textabb., V, 221 S. Wien: Springer-Verlag. 1954. S 144.—, geb. S 162.—; sfr. 24.50, 
geb. sfr. 27.70; $5.70, geb. $6.45; DM 24.—, geb. DM 27.—. 


Mit der Fertigstellung dieses Bandes hat der Verfasser seine Sammlung von Ubungsaufgaben 
abgeschlossen, die nunmehr die Mechanik des Punktes und der starren Systeme sowie die Hydro- 
mechanik umfaBt. Die Flugmechanik und Gasdynamik, die inzwischen eigene Spezialgebiete 
geworden sind, werden nicht behandelt. Da diese Aufgabensammlung an Stelle der bekannten 
von Wittenbauer gedacht ist, verdient hervorgehoben zu werden, daB Herr Prof. Federhofer 
Aufgaben, zu deren Lésung besondere Tricks erforderlich sind, soweit wie méglich vermieden 
hat, dafitr aber physikalisch und technisch wichtige Probleme bevorzugt. Dadurch ist die Aufgaben- 
sammlung ‘nicht nur eine wertvolle Erganzung der Vorlesung, sondern auch ein ausgezeichneter 
Behelf zum Studium der Hydromechanik geworden. Da die Lésungen iiberdies noch erklart 
werden, kann eine Durcharbeitung des Buches auBer Studierenden auch praktisch titigen 
Technikern groBen Nutzen bringen. F. Magyar, Wien. 


Einfiihrung in die Atomphysik. Von. W. Finkelnburg. Dritte, umgearbeitete und erweiterte 
Auflage. Mit 266 Textabb., XI, 5438S. Berlin-Gottingen-Heidelberg: Springer-Verlag. 1954. 
Geb. DM 44.—. 


Die ,,Atomphysik** des Autors ist in der Fachwelt bereits ein Begriff geworden, so daB das 
Werk tiberhaupt keiner Besprechung mehr bediirfte. Wenn ein streng wissenschaftliches Buch 
seit 1948 drei deutsche und eine amerikanische Auflage aufweist, so ist durch diese Tatsache 
allein die Nachfrage nach dem Werk bewiesen. Da das Buch dennoch gerade in der vorliegenden 
Zeitschrift ausfiihrlich besprochen wird, verfolgt den Zweck, weite Kreise der 6sterreichischen 
Ingenieure auf diese ftir sie unerschdpfliche Fundgrube moderner physikalischer Erkenntnisse 
hinzuweisen. Gerade der. Techniker neigt haufig zur Ablehnung neuerer physikalischer An- 
schauungen und Vorstellungen, da er entweder ihren praktischen Wert nicht erkennt oder aber 
die Schwierigkeiten beim Einarbeiten in das Gebiet tberschatzt. 

Der besondere Wert des vorliegenden Buches besteht nun darin, da es nicht ein einzelnes 
Spezialgebiet der Atomphysik bis ins letzte Detail behandelt, sondern den Leser in das gesamte 
Gebiet in méglichst einfacher, aber exakter Weise einftihrt. Es vermittelt ihm ein anschauliches 
Bild der Grundprobleme und Ergebnisse aller Gebiete der Atomphysik, wobei Experiment und 
theoretische Deutung in gleicher Weise zu ihrem Recht kommen. Daher mu8 besonders der 
Experimentalphysiker das Buch begriiBen, da der Autor auch die experimentellen Grundlagen 
der Atomphysik behandelt und damit die Auffassung beseitigt, daB die moderne Physik nur 
eine mathematisch-philosophische Angelegenheit sei. 

- Es ist in neuerer Zeit Mode geworden, unter dem Begriff der Atomphysik fast nur die Kern- 
physik, und darunter in erster Linie die theoretische Kernphysik, zu verstehen. Wenn auch 
die groBen Erfolge dieser physikalischen Richtung unleugbar sind, so bendtigt der moderne 
Techniker jeder Fachrichtung eine grundlegende Ausbildung, sowohl auf dem Gebiete der 
experimentellen Kernphysik als auch insbesondere auf dem Gebiete der Physik der Atomhiille. 
Daher mu8 es als besonderer Vorzug des vorliegenden Buches betrachtet werden, daB ein ebenso 
breiter Rahmen wie fiir die Physik der Atomkerne den Atomspektren, dem Atombau, der 
quantenmechanischen Atomtheorie, der Physik der Molektile und der Festkérper-Atomphysik 
zur Verfiigung steht. 

Damit ist das Buch nicht nur den Physikern, besonders natiirlich den studierenden, sondern 
auch allen in der Praxis stehenden Ingenieuren, Metallkundlern, Chemikern, Astronomen, 
Mineralogen, Kristallographen und Biologen in gleicher Weise zu empfehlen. 

DaB die gefallige Ausstattung des Buches sein Studium besonders genuBreich gestaltet, sei 
noch am Schlusse hervorgehoben. F, Regler, Wien. 


Flachentragwerke. Einfiihrung in die Elastostatik der Scheiben, Platten, Schalen und Falt- 
werke. Von K.Girkmann. Dritte, verbesserte und vermehrte Auflage. Mit 308 Textabb., 
XVII, 5588S. Wien: Springer-Verlag. 1954. Geb. S 398.—, sfr. 67.90, $ 15.80, DM 66.40. 


Das vor 8 Jahren erschienene Werk erfreut sich dank der Fiille seines in strenger, aber leicht 
verstandlicher Darstellung gebotenen Inhaltes der ihm gebithrenden weitesten Verbreitung im 
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In- und Auslande, so da®B nach der 1948 erschienenen 2. vermehrten und verbesserten Auflage 
nun bereits eine dritte Auflage notwendig geworden ist. 


Bei Beibehaltung der Gliederung in fiinf Abschnitte (Allgemeine Grundlagen der mathe- 
matischen Theorie der Elastizitat, Scheiben, Platten, Schalen, Faltwerke) erfuhr das Buch gegen- 
iiber der 2. Auflage eine Erweiterung um 56 Seiten. Dadurch war u. a. die eingehendere Be- 
handlung der Kreisscheibe und des beiderseits frei drehbar gelagerten Plattenstreifens mit be- 
liebiger Belastung mit den Sonderfallen des Plattenstreifens bei Angriff von Einzelmomenten 
méglich. In den Abschnitten iiber Scheiben und Platten konnte eine Reihe von Lésungsergebnissen 
nunmehr in geschlossener Form entwickelt werden, ferner wurden manche frither in Doppelreihen 
gewonnene Erbenisse durch einfache Reihen dargestellt, wodurch die numerische Auswertung ~ 
erleichtert wird. ; 


Der Abschnitt iiber Schalen wurde durch Hinzufiigung der von K. Girkmann im Vorjahre 
entwickelten vereinfachten Theorie des Tonnendaches mit Kreiszylinderschale nebst einem voll- 
stiindigen Zahlenbeispiel und durch Wiedergabe der Ergebnisse neuerer experimenteller Unter- 
suchungen uber die Schalenbeulung bereichert. 


Der letzte Abschnitt iiber Faltwerke erhielt eine wertvolle Erweiterung durch Vorfithrung 
der vollstandigen Berechnung eines Faltwerkturmes mit Dachscheibe 


GroBe Sorgfalt ist in der neuen Auflage auch der Vervollstaéndigung des Schrifttumverzeich- 
nisses durch Vermehrung um rund 140 neue Hinweise zugewendet worden. 


Man kann abschlieBend feststellen, daB dieses riihmlichst bekannte, mit zahlreichen guten 
Abbildungen versehene Werk in seiner neuen Bearbeitung den gegenwartigen Stand der Theorie 
der Flachentragwerke mit Ausschlu8 der Probleme groBer Formanderungen und des Durch- 
schlagproblems der Kugel- und Zylinderschale vollauf beriicksichtigt. 


Es eriibrigt sich, zum Lobe dieses vom Verlage schon ausgestatteten Buches etwas zu sagen. 
Der Unterzeichnete ist der Zustimmung vieler Fachgenossen sicher, wenn er den Wunsch aus- 
spricht, es mége diesem Buche iiber Elastostatik der Flachentragwerke in absehbarer Zeit ein 
solches uber Elastokinetik folgen. K. Federhofer, Graz. 


Grundziige der Ausgleichsrechnung nach der Methode der kleinsten Quadrate nebst Anwendungen 
in der Geodisie. Von W. GroBmann. Mit 54 Textabb., VIII, 261 S. Berlin-Gottingen- Heidelberg: 
Springer-Verlag. 1953. Geb. DM 19.80. 


Dieses Buch ist die Erweiterung einer im Jahre 1952 von W. GroBmann erschienenen 
Autographie. Die friihere Gliederung des Stoffes in die fiinf Abschnitte: Grundztige der Fehler- 
lehre, Ausgleichung von direkten Beobachtungen, Ausgleichung nach vermittelnden Beobachtungen, 
Ausgleichung nach bedingten Beobachtungen und Anhang ist beibehalten worden, doch wurde 
der Inhalt des Werkes durch die Aufnahme einiger wichtiger Kapitel weiter ausgebaut. Neben 
kleinen Erganzungen sind neu hinzugekommen die Begriindung des GauBschen Fehlergesetzes 
nach Hagen, das Entwicklungsverfahren von Boltz, eine kurze Darstellung der vermittelnden 
Beobachtungen mit Bedingungsgleichungen, Bedingungsgleichungen mit Unbekannten, die 
Anwendung aquivalenter Fehlergleichungen, die Ausgleichung durch schrittweise Ann&herung, 
die Darstellung von Funktionen durch Potenzreihen und ein Uberblick iiber die mittleren Fehler 
der Genauigkeitsma8e. Dieses Buch halt damit etwa die Mitte zwischen den einfacheren Kurz- 


ausgaben von Hegemann, Werkmeister u. a. und den umfassenderen Werken von Jordan 
und Helmert. 


Die guten Eigenschaften der ersten Veréffentlichung kennzeichnen auch diese Ausgabe mit 
ihrem klaren Aufbau und der guten Auswahl der Beispiele, wozu nun auch noch die schéne Aus- 
stattung des Werkes durch den Verlag mit einem sauberen Druck und einer iibersichtlichen 
Anordnung der Formelsysteme hinzukommt. Die vorliegende Neuherausgabe wird daher von 
den Lehrenden und Studierenden, aber auch von den Praktikern lebhaft begrii8t werden. 


F, Hauer, Wien. 


Warmespannungen infolge stationiirer Temperaturfelder. Von EH. Melan und H. Parkus. Mit 
30 Textabb., V, 114 S. Wien: Springer-Verlag. 1953. Geb. S 110.—, sfr. 18.90, $ 4.40, DM 18.50. 


Mit diesem Biichlein hat das baustatische Schrifttum eine sehr erfreuliche Bereicherung 
erfahren. Das Gebiet der Warmespannungen ist zwar im Laufe von Jahrzehnten in zahlreichen 
Einzelveréffentlichungen bearbeitet worden, die speziellen Problemen gewidmet sind. Das 
Schriftenverzeichnis dieses Buches fiihrt fast 150 soleher Arbeiten an. Doch fehlte bisher eine 
einheitliche zusammenfassende Darstellung dieser fiir die Praxis au erst wichtigen Unter- 
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suchungen. Die beiden Verfasser haben in den letzten Jahren wiederholt wertvolle Beitrage 
zur Berechnung der Warmespannungen in Scheiben, Platten und Rotationsschalen gegeben und 


sind daher fiir die hier erstmals gegebene einheitliche Darstellung der Theorie der Warme- 
spannungen besonders legitimiert. 


Nach einleitenden Bemerkungen iiber die Grundgesetze der Warmeleitung, iiber die thermisch- 
elastischen Grundgleichungen und iiber spannungsfreie Temperaturfelder und tiber die Be- 
dingungen der Spannungsfreiheit im Falle ebener Felder werden die durch Temperaturanderungen 
in einem ebenen Stabsystem entstehenden Spannungen und Verformungen berechnet. Auf den 
noch folgenden 80 Seiten dieser Schrift werden dann zahlreiche Probleme der Warmespannungen 
in Scheiben, Platten, Umdrehungskérpern und in diinnen Rotationsschalen behandelt und auch 
Probleme der Wiarmespannungen in Kérpern mit Einschliissen erértert. 


Vom Leser wird die Kenntnis der wichtigsten Tatsachen der Elastizitatslehre und der Be- 
rechnung statisch unbestimmter Konstruktionen sowie einiges iiber die Lésung_partieller 
Differentialgleichungen und aus der Potentialtheorie vorausgesetzt. In einem folgenden Bande 
beabsichtigen die Verfasser auch die nichtstationdéren Temperaturfelder zu behandeln und einiges 
uber elastisch-plastische Spannungszustinde zu bringen. 


Das Buch zeichnet sich durch eine klare Darstellung aus und kann nicht nur dem ent- 
werfenden Ingenieur, sondern auch den Hdérern héherer Hochschulsemester zum Studium 
warmstens empfohlen werden. K. Federhofer, Graz. 


Lehrbuch des Stahlbetonbaues. Grundlagen und Anwendungen im Hoch- und Briickenbau. Von 
A. Pucher. Zweite, verbesserte und vermehrte Auflage. Mit 321 Textabb., X, 331 S. Wien: 
Springer-Verlag. 1953. Geb. S 192.—, sfr. 32.70, $7.60, DM 32.—. 


Auf die Vorziige dieses sehr modernen Lehrbuches habe ich schon in meiner Buchbesprechung 
zur ersten Auflage hingewiesen und wiederhole, da8 Pucher bestrebt ist, alle Probleme grund- 
sdtzlich méglichst allgemein zu erfassen und unter Einsatz seiner ausgezeichneten Beherrschung 
mathematischer Berechnungsverfahren zu lésen. Dies gilt ebenso fiir die in der vorliegenden 
zweiten Auflage neu hinzugekommenen Abschnitte. 


Im Kapitel tiber das Biegungsproblem des Stahlbetons im plastischen Bereich sind der 
allgemeinen Lésung die aus den Vorschriften der ONORM B 4200/4. Teil folgenden Spezialisie- 
rungen zugefiigt und neue Bemessungstafeln entwickelt worden. So notwendig ein dringender 
Hinweis auf die héhere Verantwortung der planenden und ausfiihrenden Ingenieure bei An- 
wendung dieses neuen Traglastverfahrens zufolge der damit gegebenen hoheren Baustoffausniitzung 
auch sein mag, sollte doch gerade in einem Lehrbuch eine Kritik, wie sie der letzte Satz auf Seite 137 
ausspricht, nicht aufscheinen. 


Neu eingefiigt ist die Behandlung des Knickproblems im plastischen Bereich des Stahlbetons 
mit einer ausfiihrlichen Darstellung der Theorie von Engesser und darauf aufbauend im Abschnitt 
, Briickenbau“* das Stabilitatsproblem schlanker Bogen. Ferner ein Kapitel iiber die Berechnung 
von ,,Bogen mit stiitzliniennaher Ache“‘ und eine Erganzung des Kapitels Torsion fur gleich- 
zeitiges Auftreten von Schub- und Verdrehungsbeanspruchungen. 


Nach Inhalt und Ausstattung mu8B dieses Werk zu den besten Neuerscheinungen der letzten 
Jahre gezihlt werden, das allerdings vom Leser ein tiberdurchschnittliches MaB von Vorkennt- 
nissen voraussetzt. Auch die vorliegende zweite Auflage wird viele Freunde finden. 


F.. Pongratz, Wien. 


Elektrische Triebfahrzeuge. Ein Handbuch fiir die Praxis sowie fiir Studierende. Von K. Sachs. 
In zwei Banden. Erster Band: Mit 847 Textabb., XV, 700 S. — Zweiter Band: Mit 850 Text- 
abb. und 16 Tafeln, XII, 696 S. Frauenfeld: Verlag Huber & Co. AG. 1953. Geb. sfr. 65.—. 


Im Jahre 1928 erschien im Springer-Verlag in Berlin das Buch ,,Elektrische Vollbahn- 
lokomotiven‘‘ von Prof. Sachs und zehn Jahre spater sein Werk itiber ,,Die ortsfesten Anlagen 
elektrischer Bahnen‘‘ im Verlag Orell FiBli in Ziirich. Nun hat der Verfasser ein zweibandiges 
Werk geschaffen, das umfassend die elektrischen Triebfahrzeuge behandelt. 


Der erste Band bringt die theoretischen Grundlagen: Zugkraft, Bremskraft und Leistung 
und den mechanischen Teil: Gesamtaufbau, Laufwerk und Triebwerk. Im zweiten Band findet 
man die elektrische Ausriistung: Triebfahrzeuge fiir Gleichstrom, Einphasen-Wechselstrom, 
Drehstrom und Umformer-Triebfahrzeuge. In eigenen Abschriften werden die Triebfahrzeuge 
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fiir Zahnstangenstrecken besprochen, dann der Trolleybus und die Speicher-Triebfahrzeuge. 
Fiinfzehn Tafeln mit ausgefiihrten Triebfahrzeugen bilden den Abschlu8. 


Hier liegt ein Werk vor, das mit groBer Sachkenntnis und mustergiiltiger Sorgfalt angelegt 
ist. Man sehe sich nur die 1700 Zeichnungen und Abbildungen an, oder die Hinweise auf das 
Schrifttum und die Patente! 


Es soll nicht unerwahnt bleiben, daB der Preis dieses Werkes ungewohnlich niedrig ist. Dies 
wurde durch eine Subventionierung durch schweizerische Industrieunternehmungen, Bahnen, 
Elektrizitatslieferungsgesellschaften usw. erméglicht. 


Das Werk wird durch viele Jahre das Standardwerk tiber die elektrischen Triebfahrzeuge 
bleiben, ein Lehr- und Handbuch, nach dem jeder greifen muB, der eine eingehende Belehrung 
sucht oder auf eine besondere Frage eine erschépfende Antwort wiinscht. Der Verfasser ist zur 
Vollendung seiner Arbeit zu begliickwiinschen. H. Sequenz, Wien. 


Autorenverzeichnis. 


Basch, A., Uber Schwingungen von Systemen 
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(Buchbesprechung.) 8. 308. 


Czitary, E., Model Analysis of Structures. Von 
T.M. Charlton. (Buchbesprechung.) 8.308. 


Desoyer, K., s. Heinrich, G. 


Federhofer, K., Technische Dynamik. Von 
C. B. Biezeno und R. Grammel. Zweite 
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einzelne Maschinenteile. Zweiter Band: 
Dampfturbinen und Brennkraftmaschinen. 
Knicklast der axial gedriickten Kreis- 
zylinderschale bei Vorhandensein eines ent- 
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elastischen Widerstandes. 8. 90. 
Vorlesungen tiber Baustatik. Zweite Auf- 
lage. Erster Band: Statisch bestimmte 
Systeme. Spannungsberechnung. 
Elastische Formaéanderungen. Stabili- 
tatsprobleme. — Seile. Von F. Stiissi. 
(Buchbesprechung.) 8. 229. 
Flachentragwerke. Von K. Girkmann. 
Dritte Auflage. (Buchbesprechung.) 8. 309. 
Warmespannungen infolge stationarer Tem- 
peraturfelder. Von E. Melan und H. Par- 
kus. (Buchbesprechung.) S. 310. 

Frank, W., Zur Berechnung von Potential- 

str6mungsfeldern. 8S. 97. 
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Geiringer, Hilda, Bemerkung zur Theorie der 
Charakteristiken. S. 107. 

Glaser, W., Principles of Electronics. Von 
L. T. Agger. (Buchbesprechung.) S. 76. 

—, Eine neue Begriindung der wellenmechani- 
schen Elektronentheorie. S. 110. 

—, Die Lecher-Leitung in der Kurzzeit-Impuls- 
technik. 8. 121. 


Hauer, F., Grundziige der Ausgleichsrechnung 
nach der Methode der kleinsten Quadrate 
nebst Anwendungen in der Geodasie. Von 
W. Gro8Bmann. (Buchbesprechung.) 8.310. 


Heinrich, G. und K.Desoyer, Die Kipp- 
steifigkeit eimseitig wirkender Strémungs- 
spurlager. 8. 138. 


Inzinger, R., Mathematics for Students of 
Engineering and Applied Science. Von L. B. 
Benny. (Buchbesprechung.) 8. 307. 

—, Handbook of Elliptic Integrals for Engineers 
and Physicists. Von P. F. Byrd und M. D. 
Friedman. (Buchbesprechung.) 8. 308. 


Kammerer, C., Thermodynamische Grund- 
lagen. Von R. Plank. (Buchbesprechung.) 
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—, Str6émung in einer Expansionsdtise mit 
Reibung. S. 293. 

—, H. Mache gestorben. (Persénliches.) 8. 307. 

Kirste, L., Abwickelbare Verformung dinn- 
wandiger Kreiszylinder. S. 149. 

Korst, H., Auflésung eines ebenen Freistrahl- 
randes bei Berticksichtigung der urspriing- 
lichen Grenzschichtstro6mung. 8. 152. 

Kracmar, F., Das magnetische Feld senkrecht 
gekreuzter, stromdurchflossener Leiter. 
S. 158. 

Krzywoblocki, M. Z. E., Bergman’s Linear 
Integral Operator Method in the Theory of 
Compressible Fluid Flow. C. Axially Sym- 
metric Flow and Singularities. 8. 237. 


Mache, H., Uber die Unerreichbarkeit des 
absoluten Nullpunktes und tiber die Normal- 
form einer jeden Zustandgleichung. 8. 161. 

Magyar, F., Hydraulik. Ihre Grundlagen und 
praktische Anwendung. Von J. Kozeny. 
(Buchbesprechung.) S. 228. 

—, Aufgaben aus der Hydromechanik. Von 
K. Federhofer. (Buchbesprechung.) 8.309. 

Melan, E., Warmespannungen infolge eines 
rotierenden Temperaturfeldes. 8. 165. 

Melan, H., Prof. Dr. Franz Magyar 60 Jahre. 
Sa79e 

Miller, W., Beitrag zur Biegungstheorie der 
Mehrpilzplatte. S. 1. 

—, Die Bewegung eines Rotationskérpers in 
der reibungslosen Flissigkeit und das in- 
stabile Moment der Druckkrafte. S. 171. 

—, Zur Bestimmung der Tragheitskoeffizienten 
unsymmetrischer Rotationskérper. S. 263. 

21* 


€ 


314 


Oswatitsch, K. und L. Sjédin, Kegelige 
Uberschallstro6mung in Schallnihe. 8. 284. 


Parkus, H., Das Anlaufen einer Schubdiise mit 
vorgeschaltetem Rohr. S. 185. 

Pongratz, F., Die Bemessung der Stahlbeton- 
bauteile nach Onorm B 4200. Von K. Jager. 
4. Teil. (Buchbesprechung.) S. 77. 

—, Lehrbuch des Stahlbetonbaues. Von A. 
Pucher. Zweite )Auflage. (Buchbespre- 
chung.) 8. 311. ° 

Proll, A., Beitrage zum Schwingenflugproblem. 
8. 189. 

Puwein, M. G., Die mutierte Pendellange. 8.54. 


Regler, F., Einfiihrung in die Atomphysik. Von 
W. Finkelnburg. Dritte Auflage. (Buch- 
besprechung.) S. 309. 

Richardson, E. G. and R. I. Tait, Measure- 
ments on the Damping of Flexural Vibra- 
tions in Rods Immersed in Liquids. S. 200. 

Richter, W., Koordinatentransformationen 
mit Hilfe eines Fluchtliniennomogramms 
und Anwendungen auf die graphische L6- 
sung von Differentialgleichungen. S. 39. 
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S. 207. 
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Signorini, A., Uber eine Erweiterung der 
linearisierten Theorie der Elastizitat. S. 47. 

Sjédin, L., s. Oswatitsch, K. 
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Tremmel, E., Zur Theorie kreisberandeter 
Bogenscheiben. S. 11. 
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Springer-Verlag in Wien I 


Soeben erschienen: 


Die Grundlagen der Akustik. Von Dipl.-Ing. Dr. rer. nat. E. Skudrzyk, B. Sc. (Eng.), 
A. C. G.I., a. o. Professor fiir Niederfrequenztechnik an. der Roemuchen Hochschule in Wien. 
Mit. 450 PesreCu sone. XXII, 1084 Seiten. 1954. 

~Ganzleinen S 880.—, DM 147.—, § 35.—, sfr. 150.— 

Mit Sictein Werk legt der Verfasser ein Grundla genbuch vor, das mit der gebotenen Ausfiihrlichkeit in die vielgestaltigen 


Probleme der modernen Akustik einfiihrt und, dariiber hinausgehend, das groBe Gebiet so erschépfend darstellt, als es dem 
heutigen Stande der Forschung entspricht. Der Verfasser bedient sich in der Darstellung moderner Methoden und bringt 


eine Reihe zum Teil noch unveroffentlichter Ergebnisse eigener Arbeiten. Das Buch wendet sich nicht nur an den Akustiker, 


sondern an den Ingenieur und Physiker im weiteren Sinn. Fragen der reinen Mechanik werden iibergangen, dagegen hat der 
Vertasser besondere Sorgfalt auf die konsequente Durchfiihrung simtlicher Ableitungen und die Auswahl eines tibersichtlichen 
Systems von Bezeichnungen gelegt. Zur Erleichterung des Studiums ist jedem Kapitel eine kurze Zusammenfassung voran- 
gestellt. Die Grundliteratur wird in einem umfassenden Verzeichnis mit iiber 4000 Zitaten angefiihrt. 


Theorie ideal plastischer Korper. v Von W. Prager, Eng. D., Professor of Applied Mechanics, 
Brown University, Providence, R. I., USA., und P. G. Hodge, Jr., Ph. D., Associate Professor 
of Applied Mechanics, Polytechnic Tiistitute of Brooklyn, N. Y., USA. Ins Deutsche tibertragen 
von F, Chmelka, Dr. phil., Dr. techn., Privatdozent an der Technischen Hochschule in 
Wien. Mit 97 Re X, 274 Seiten. 1954. 

Ganzleinen S 198.—, DM 33.—, $ 7.85, ‘str. 33.80 


In der ietzten Zeit gewinnt die Plastizititstheorie immer mehr an Bedeutung. Sie macht Vorgange, wie Walzen, Ziehen 
Pressen usw., der Berechnung zuginglich, sie erméglicht cine der Wirklichkeit entsprechende Bestimmung des Sicherheitsfaktors 
gegen Kollaps - von Tragwerken und anderes mehr. Die vorliegende Ubersetzung des grundlegenden Werkes aus dem Englischen 
ins Deutsche schlie8t somit eine empfindliche Liicke, indem sie dieses wichtige Wissensgebiet nun auch dem deutschsprachigen 
Leser- und Forscherkreis zuginglich macht. Das Buch geht, was Vorkenntnisse anbelangt, von durchaus elementaren Voraus- 
setzungen aus, und bringt nicht nur eine Darstellung der bereits seit langem gelésten Grundaufgaben, sondern umfa8t auch 
die neuesten Erkenntnisse der Theorie. 


Einfihrung. in die Statik. Von Dr. phil., Dr. oun F. Chmelka, Privatdozent an der 
Technischen Hochschule in Wien, und Dr. techn. E. Melan, o. Professor an der Techni- 
schen Hochschule in Wien, wirkl. Mitglied der Osterreichischen Akademie der Wissenschaften. 
8 echste, umgearbeitete underweiterte Auflage. Mit 161 Textabbildungen. VIL, 210 Seiten. 1954. 

Steif geheftet S 84.—, DM 14.—, $ 3.35, sfr. 14.40. 


Wiihrend die letzten beiden Ausgaben des bekannten Hinfithrungsbuches unverinderte Nachdrucke der dritten Auflage waren, 
was durch die Nachkriegsschwierigkeiten bedingt war, ist die jetzt erscheinende sechste Auflage griindlich bearbeitet 
und wesentlich erginzt worden. Die theoretischen Abschnitte sind wesentlich vertieft und verbreitert und zahlreiche neue 
Beispiele hinzugefiigt worden, wodurch die didaktischen Vorziige der Darstellung noch erhéht wurden. Neu ist ein Abschnitt, 
in dem einige statisch bestimmte Tragwerksformen besprochen sind, deren Behandlung erfahrungsgemifZ dem Anfanger 
Schwierigkeiten bereiten. 


Methoden und Probleme der Wettervorhersage. Von Dr. H. Reuter, Privatdozent 
an der Universitat Wien, Observator an der Zentralanstalt fiir Meteorologie und Geodynamik 
in Bn Mit 46 Tome dun ee VIII, 161 Seiten. 1954. 

Ganzleinen S 132.—, DM 22.—, $ 5.25, sfr. 22.60 


“Das Buch will einen Uberblick iiber die wicitigaten Vorhergagemethoden vermitteln, wobei den modernen Methoden und 


Problemen besondere Beachtung geschenkt wurde, da iiber diesen Gegenstand noch keine zusammeniassenden Darstellungen 
vorliegen. Es tritt dabei naturgema8 die theoretische Seite der Problematik der Wetterprognose starker in Erscheinung als 
in fritheren Darstellungen iiber denselben Gegenstand. Es wird jedoch versucht, mit einem mdoglichst geringen Aufwand an 
Mathematik einem groBeren Leserkreis die nicht immer ganz einfachen Probleme zugénglich zu machen. 5 


Beschaftigungspolitik. Von Dr. W. Taucher, 0. Professor an der Universitét Graz. V, 
126 Seiten. 1954. Steif geheftet S 36.—, DM 6.—, $ 1.45, sfr. 6.20 


» ++. Professor Taucher gibt in seinem Buch einen sehr umfassenden Uberblick iiber die vielseitige Problematik der Be- 
schiiftigungspolitik. Er verzichtet darauf, die einzelnen Seiten des Problems bis in alle Hinzelheiten zu analysieren. Er gewinnt 
dadurch die Méglichkeit, in einer auch einem breiteren Publikum zuginglichen Form, die wichtigsten Gesichtspunkte und 
Erkenntnisse besonders zu unterstreichen. .. .‘‘ Dee Industrie 
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Die Wicklungen elektrischer Maschinen. Von Heinrich Sequenz, Dipl.-Ing., 
Dr. techn., Dr.-Ing., Dr. phil., Professor an der Technischen Hochschule in Wien, korresp. 
Mitglied der Osterreichischen Akademie der Wissenschaften. In vier Banden. 


Soeben erschien: 


Dritter Band: Wechselstrom-Sonderwicklungen. Mit 578 Textabbildungen. XIX, 397 Seiten. — 
1954. S 420.—, DM 70.—, $ 16.65, sfr. 71.60; Gln. § 438.—, DM 73.—, $ 17.40, sfr. 74.80 


Prither erschienen: 4 
Erster Band: Wechselstrom-Ankerwicklungen. Mit 408 Textabbildungen. XX, 365 Seiten. 
1950. S 228.—, DM 37.—, $ 9.—, sfr. 39.—; Gin. S 248.—, DM 40.—, $9.90, sfr. 42.50 

Zweiter Band: Wenderwicklungen. Mit 423 Textabbildungen. XVI, 331 Seiten. 1952. 
S 325.—, DM 54.—, $ 12.90, sfr. 56.—; Gin. S 343.—, DM 57.—, $ 13.50, str. 59.— 


In Vorbereitung : 


Vierter Band: Herstellung der Wicklungen. 


,,... Sequenz hates verstanden, den immerhin spréden Stoff so zu meistern, da8 es dem Ingenieur, der sich mit dem Entwurf 
von Wechselstromwicklungen befaBt, méglich ist, das Buch mit Nutzen durchzustudieren. Das gilt auch fiir einzelne Studierende 
des Elektromaschinenbaues, die das Bediirfnis fiihlen, sich etwa gelegentlich einer Studien- oder Diplomarbeit in die Materie 
zu vertiefen...‘* - ree : VDI-Zeitschrift 


Handbuch der mikrochemischen Methoden. Herausgegeben von Friedrich 
Hecht, Professor an der Universitat Wien, und Michael K. Zacherl, Professor an de 
Tierarztlichen Hochschule in Wien. In fiinf Banden. 


Jeder selbstindig erscheinende Bandteil bew. Band ist einzeln kéuflich. Bei Verpflichtung zur Abnahme des 
Gesamtwerkes sowie bei Vorbestellung der einzelnen Teile ermépigt sich der endgiiltige Ladenpreis um 20% 


Soeb. hien: 
et cat Band I / Teil 4 


Praparative Mikromethoden in der organischen Chemie. Von Prof. Dr. Hans Lieb, Graz, und 
Dipl-Ing. Dr. Wolfgang Schéniger, Graz. ; sie ‘ ; 
Mikroskopische Methoden. Von Prof. Dr. Ludwig Kofler} und Dr. Adelheid Kofler, Innsbruck. — 
Mit 275 Textabbildungen VI, 236 Seiten. 1954 = - . 


| Bei Subskription auf das Gesamtwerk: Ganzleinen Si228.—, DM 38.—, $ 9.05, sfr. 38.90 i as 
Endgiiltiger Ladenpreis: Ganzleinen S 285.—, DM 47.50, $ 11.30 sfr. 48.60 — 


Demmniichst erscheint: = 
Band Il / Teil 1 e. 
Radiochemische Methoden der Mikrochemie. Von Priv.-Doz. Dr. Engelbert Broda, Wien, und 
Dr. Thomas Schénfeld, Wien. ; 
Messung radioaktiver Strahlen in der Mikrochemie. Von Prof. Dr. Berta Karlik, Wien, Dr. Traude 
Bernert, Wien, und Priv.-Doz. Dr. Karl Lintner, Wien. = 
Photographische Methoden in der Radiochemie. Von Dr. Hanne Lauda, Wien. 
Mit etwa 90 Textabbildungen. Etwa 330 Seiten. 1955 


| Bei Vorausbestellung bis zum Erscheinen: Ganzleinen S 388.—, DM 64.80, $ 15.45, sfr. 66.40 | 
Endgiiltiger Ladenpreis nach Erscheinen: Ganzleinen § 486.—, DM 81.—, $ 19.30, sfr. 83.— 


he Bei Verpflichtung zur Abnahme des gesamten Handbuches gilt der Vorbestellpreis 
weiter als Subskriptionspreis 


In der Folge werden erscheinen: 


Band I, Teil 2: Biochemische Methoden einschlieBlich medizinischer Verfahren. Von Th. Leipert, 
Wien. Lebensmittelmikrochemische Methoden. Von F. Miinchberg, Wien, und F. Zaribnicky, Wien. 
Band II, Teil 2: Polarographie. Von H. Hohn, Wien, und Otilie Schlager, Wien. 
Band ITI: Anerganisch-analytische Methoden. Von F. Hecht, Wien. ; 
Band IV: Organisch-analytische Methoden. Von M. K. Zacherl, Wien. , 
pid V: Mikromethoden zur Bestimmung physikalisch-chemischer Konstanten. Von Martha Sobotka, 
TAZ. a 
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